[MITgcm-devel] Beaufort experiment on mac os x

[Menemenlis, Dimitris (3248)]

I add to Beaufort thread a message I just received from Torge.

Torge, thank you for notes below.
I am in process of testing/updating the Beaufort experiment
and will get back later today or tomorrow.

Cheers

On Mar 13, 2012, at 1:04 PM, Torge Martin wrote:

> Hi Dimitris,
> 
> I tested the lab_sea set up on my desktop machine and it works fine when I use option "-ieee" with genmake2, otherwise I get a segmentation fault when running.
> 
> Considering this recent experience I tried again running the beaufort set up. However, I ran into errors when compiling which are most likely due to recent changes in the general code. Could you please help to fix these issues? I have found a few solutions myself (see list below).
> 
> Thanks,
> Torge
> 
> 
> =======================================
> COMPILE TIME ERRORS
> 
> ---------------------------------------
> ERROR:
> cpl_mpmice.for:3838.60:
> 
>       Real*8 TICES      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,MULTDIM,nSx,nSy)   
>                                                             1
> Error: Symbol 'multdim' at (1) has no IMPLICIT type
> make: *** [cpl_mpmice.o] Error 1
> 
> SOLUTION:
> in beaufort/code/cpl_mpmice.F
> add 
> # include "SEAICE_SIZE.h" 
> in line 33
> 
> ---------------------------------------
> ERROR:
> do_oceanic_phys.for:3133.32:
> 
>       IF ( allowInteriorFreezing ) THEN                                 
>                                 1
> Error: Symbol 'allowinteriorfreezing' at (1) has no IMPLICIT type
> make: *** [do_oceanic_phys.o] Error 1
> 
> SOLUTION:
> in beaufort/code/do_ocean_phys.F 
> replace
>  C--   Freeze water in the ocean interior and let it rise to the surface
>  C     temporarily exclude from adjoint computations until
>  C     impact has been vetted in forward integrations
>        IF ( allowInteriorFreezing ) THEN
>  #ifdef ALLOW_AUTODIFF_TAMC
>  CADJ STORE theta, salt = comlev1, key = ikey_dynamics,
>  CADJ &     kind = isbyte
>  CADJ STORE recip_hfacc = comlev1, key = ikey_dynamics,
>  CADJ &     kind = isbyte
>  #endif
>          CALL FREEZE_INTERIOR(  myTime, myIter, myThid )
> 
> with
>  #ifdef ALLOW_SALT_PLUME
>        IF ( useSALT_PLUME ) THEN
>            CALL SALT_PLUME_DO_EXCH( myTime, myIter, myThid )
>  #endif /* ALLOW_SALT_PLUME */
> 
> >> I found this by looking into MITgcm/model/src/do_ocean_phys.F
> 
> 
> 
> =======================================
> RUN TIME ERRORS
> 
> ---------------------------------------
> ERROR:
> (PID.TID 0000.0001)  SEAICE_READPARMS: finished reading data.seaice
> S/R SEAICE_READPARMS: "SEAICE_freeze" no longer allowed in file "data.seaice"
> S/R SEAICE_READPARMS: set instead "SEAICE_tempFrz0" and "SEAICE_dTempFrz_dS"
> S/R SEAICE_READPARMS: "SIsalFrac" is no longer allowed in file "data.seaice"
> S/R SEAICE_READPARMS: set "SEAICE_saltFrac" instead
> S/R SEAICE_READPARMS: Error reading parameter file "data.seaice"
> some out of date parameters were found in the namelist
> STOP ABNORMAL END: S/R SEAICE_READPARMS
> 
> SOLUTION:
> change beaufort/input/data.seaice:
> replace SEAICE_freeze with SEAICE_tempFrz0 using same value
> replace SIsalFRAC with SEAICE_saltFrac using same value
> 
> >> I don't know which value to use for SEAICE_dTempFrz_dS; I guess default doesn't match the value chosen for SEAICE_tempFrz0
> 
> ---------------------------------------
> ERROR:
> (PID.TID 0000.0001) SEAICE_CHECK: #define ALLOW_SEAICE
> SEAICE_CHECK: #define ALLOW_SEAICE
> SEAICE_CHECK: turbulent ice-ocn heat flux default changed.
> SEAICE_CHECK:    To recover the old default : set
> SEAICE_CHECK:    SEAICE_mcPheePiston to the first ocn level
> SEAICE_CHECK:    thickness divided by SEAICE_deltaTtherm.
> SEAICE_CHECK: SEAICE_DO_OPEN_WATER_GROWTH / MELT options are obsolete
> they are replaced with run time parameter SEAICE_doOpenWaterGrowth / Melt
> STOP ABNORMAL END: S/R SEAICE_CHECK
> 
> >> solution under debate ???
>