[MITgcm-devel] Beaufort experiment on mac os x

[Art Lazanoff]

Hi all,

I can certainly test drive the code using ifort v12 update 3.


Art

****************************************************************************
*
Art Lazanoff                                     NASA Ames Research Center
Senior Scientific Consultant                       Mail Stop 258-6
Arthur.S.Lazanoff at nasa.gov (650) 604-1687        Moffett Field, CA
94035-1000
****************************************************************************
*


> -----Original Message-----
> From: Menemenlis, Dimitris (3248)
> [mailto:Dimitris.Menemenlis at jpl.nasa.gov]
> Sent: Tuesday, March 13, 2012 5:30 PM
> To: MITgcm-devel at mitgcm.org
> Cc: Lazanoff, Arthur S. (ARC-TN)[Computer Sciences Corporation]; Torge
> Martin
> Subject: Re: [MITgcm-devel] Beaufort experiment on mac os x
> 
> 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 ???
> >