[MITgcm-devel] Beaufort experiment on mac os x

[Menemenlis, Dimitris (3248)]

Torge, sorry that sea ice code is changing so much.
Hopefully things will stabilize now.

I have updated Beaufort experiment to work with latest
MITgcm code, including the fixes that you mention below.

You can run the experiment using the instructions from
README_beaufort.txt for convenience.
I will keep digging for the segmentation-fault
large-memory footprint problem.

>>>>>>>>>>>

  cvs co -d beaufort MITgcm_contrib/MPMice/beaufort
  cd beaufort
  cvs co MITgcm_code

>> Obtain copies of following directories and put them in beaufort
>>  ftp://ecco2.jpl.nasa.gov/data1/beaufort/run_template
>>  ftp://ecco2.jpl.nasa.gov/data1/data/atmos/blend_forcing/cube78_forcing

 cd MITgcm
 mkdir bin exe
 cd bin
 \cp ../../code/* .
 \mv CPP_EEOPTIONS.h_UNCOUPLED CPP_EEOPTIONS.h
 ../tools/genmake2
 make depend
 make -j 16

 cd ../exe
 cp ../../run_template/* .
 cp ../../input/* .
 cp ../bin/mitgcmuv .
 ./mitgcmuv >& output.txt &

<<<<<<<<

Dimitris Menemenlis

On Mar 13, 2012, at 5:29 PM, Menemenlis, Dimitris (3248) wrote:

> 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 ???
>> 
>