[MITgcm-support] Low salinity bias in PI experiment based on cpl_aim+ocn
Jean-Michel Campin
jmc at mit.edu
Tue Jan 20 22:22:57 EST 2026
Hi Fisher,
I assume "PI" is for pre-industrial. I took a quick look at the 2 figures, but the land/sea
mask does not look like the one I am used to, on cs32 grid (e.g, the red-sea seem closed),
see figure attached. Are you using same bathy and orography as in verificaton/cpl_aim+ocn/ ?
The other thing is that I am not sure what the different deltaT setting you have in
input_ocn/data (here in: input_ocn_data.txt)
> deltaTmom =3600.,
> deltaTtracer=3600.,
> deltaTfreesurf= 3600.,
> deltaTClock =86400.,
implies in a couple set-up (it's a bit unusual). Did you check what this is doing ?
And regarding known limitations:
1) low vertical resolution in the atmosphere.
2) the contrast in OLR (Outgoing Long-Wave radiation) map tends to be weaker
than what we could expect.
3) if I remember well (it was a long time ago), the AMOC was weak in this coupled set-up,
a partial re-routing of some of the run-off helped, but I don't remember getting low
salinity in the Med.
May be others have better advise.
Cheers,
Jean-Michel
On Thu, Jan 15, 2026 at 06:22:48AM +0000, ht wrote:
> > Dear MITgcm Support Community,
> >
> > I am running coupled atmosphere???land???ocean experiments based on the cpl_aim+ocn configuration.
> > My goal is to validate the model setup using a fully coupled PI-type experiment without SST/SSS relaxation or prescribed surface fluxes, before proceeding to paleoclimate simulations. The experiment is conducted with aim_fixed_pCO2 = 285e-6, following a PI control configuration.
> >
> > As the simulation progresses (after about 500 years), some parts of the ocean (e.g., the Southern Ocean and the Mediterranean Sea) become unrealistically fresh, as shown in the attached figure, and North Atlantic Deep Water (NADW) formation appears to cease.
> >
> > I would appreciate guidance on the following points:
> >
> > -
> >
> > Are there known issues or common pitfalls in cpl_aim+ocn that could lead to excessive ocean freshening?
> >
> > -
> >
> > In PI-control fully coupled experiments, is some form of SST/SSS relaxation generally required to maintain stable NADW formation?
> >
> > -
> >
> > When SST/SSS relaxation is turned off, are there specific parameters in data.land and/or data.aimphys that should be adjusted to avoid such behavior?
> >
> > Any suggestions on likely causes, useful diagnostics, or where to focus debugging efforts would be greatly appreciated.
> >
> > Thank you very much for your time and support.
> >
> > Best regards,
> > Fisher
> # Land package parameters:
>
> #-- Land model parameter (namelist LAND_MODEL_PAR):
> # land_calc_grT :: step forward ground Temperature
> # land_calc_grW :: step forward soil moiture
> # land_impl_grT :: solve ground Temperature implicitly
> # land_calc_snow :: step forward snow thickness
> # land_calc_alb :: compute albedo of snow over land
> # land_oldPickup :: restart from an old pickup (=before checkpoint 52l)
> # land_grT_iniFile :: File containing initial ground Temp.
> # land_grW_iniFile :: File containing initial ground Water.
> # land_snow_iniFile :: File containing initial snow thickness.
> # land_deltaT :: land model time-step
> # land_taveFreq :: Frequency^-1 for time-Aver. output (s)
> # land_diagFreq :: Frequency^-1 for diagnostic output (s)
> # land_monFreq :: Frequency^-1 for monitor output (s)
> # land_dzF :: layer thickness
> &LAND_MODEL_PAR
> land_calc_grT =.TRUE.,
> land_calc_grW =.TRUE.,
> land_impl_grT =.TRUE.,
> land_calc_snow=.TRUE.,
> land_calc_alb =.TRUE.,
> land_grT_iniFile='land_grT_ini.cpl.bin',
> land_grW_iniFile='land_grW_ini.mn.bin',
> land_snow_iniFile='land_snw_ini.cpl.bin',
> land_deltaT=450.,
> # land_oldPickup=.TRUE.,
> land_dzF= 0.1, 4.0,
> land_diagFreq=311040000.,
> land_monFreq=0.,
> &
>
> #-- Physical constants (namelist LAND_PHYS_PAR):
> # land_grdLambda :: Thermal conductivity of the ground (W/m/K)
> # land_heatCs :: Heat capacity of dry soil (J/m3/K)
> # land_CpWater :: Heat capacity of water (J/kg/K)
> # land_wTauDiff :: soil moisture diffusion time scale
> # land_waterCap :: field capacity per meter of soil
> # land_fractRunOff:: fraction of water in excess which run-off
> # land_rhoLiqW :: density of liquid water (kg/m3)
> # land_rhoSnow :: density of snow (kg/m3)
> # land_Lfreez :: Latent heat of freezing (J/kg)
> # land_hMaxSnow :: Maximum snow-thickness (m)
> # diffKsnow :: thermal conductivity of snow (W/m/K)
> # timeSnowAge :: snow aging time scale (s)
> # hNewSnowAge :: new snow thickness that refresh the snow-age (by 1/e)
> # albColdSnow :: albedo of cold (=dry) new snow (Tsfc < -10)
> # albWarmSnow :: albedo of warm (=wet) new snow (Tsfc = 0)
> # albOldSnow :: albedo of old snow (snowAge > 35.d)
> # hAlbSnow :: snow thickness for albedo transition: snow/ground
> &LAND_PHYS_PAR
> # land_Lfreez=0.,
> land_hMaxSnow=10.,
> &
> # Coupling package parameters, ATM component:
> # cpl_oldPickup :: restart from an old pickup (= until ckpt 59h)
> # useImportMxlD :: True => use Imported Mix.Layer Detph from coupler
> # useImportSST :: True => use the Imported SST from coupler
> # useImportSSS :: True => use the Imported SSS from coupler
> # useImportVsq :: True => use the Imported Surf. velocity^2
> # useImportThSIce :: True => use the Imported thSIce state var from coupler
> # useImportFlxCO2 :: True => use the Imported air-sea CO2 flux from coupler
> # cpl_atmSendFrq :: Frequency^-1 for sending data to coupler (s)
> # maxNumberPrint :: max number of printed Exp/Imp messages
> &CPL_ATM_PARAM
> cpl_oldPickup =.TRUE.,
> useImportMxlD =.TRUE.,
> useImportSST =.TRUE.,
> useImportSSS =.TRUE.,
> useImportVsq =.TRUE.,
> # useImportThSIce=.TRUE.,
> # useImportFlxCO2=.TRUE.,
> cpl_atmSendFrq = 86400.,
> &
> # ====================
> # | Model parameters |
> # ====================
> #
> # Continuous equation parameters
> &PARM01
> tRef=15*20.,
> sRef=15*35.,
> viscAh =3.E5,
> viscAr =1.E-3,
> diffKhT=0.,
> diffK4T=0.,
> diffKrT=3.E-5,
> diffKhS=0.,
> diffK4S=0.,
> diffKrS=3.E-5,
> gravity=9.81,
> rhoConst=1030.,
> rhoConstFresh=1000.,
> eosType='JMD95Z',
> #allowFreezing=.TRUE.,
> ivdc_kappa=10.,
> implicitDiffusion=.TRUE.,
> implicitFreeSurface=.TRUE.,
> exactConserv=.TRUE.,
> select_rStar=2,
> nonlinFreeSurf=4,
> hFacInf=0.2,
> #hFacInf=0.1,
> hFacSup=2.0,
> useRealFreshWaterFlux=.TRUE.,
> temp_EvPrRn=0.,
> hFacMin=.1,
> hFacMinDr=20.,
> vectorInvariantMomentum=.TRUE.,
> staggerTimeStep=.TRUE.,
> readBinaryPrec=64,
> writeBinaryPrec=64,
> debugLevel=0,
> plotLevel=0,
> useSingleCpuIO=.TRUE.,
> # tempAdvScheme=77,
> # saltAdvScheme=77,
> # multiDimAdvection=.TRUE.,
> &
>
> # Elliptic solver parameters
> &PARM02
> cg2dMaxIters=200,
> cg2dTargetResidual=1.E-9,
> # cg2dTargetResWunit=1.E-14,
> &
>
> # Time stepping parameters
> &PARM03
> # pickupSuff='ckptA',
> #0yr
> # nIter0= 0,
> startTime=26282880000.0,
> endTime=26438400000.0,
> # deltaT= 3600.,
> deltaTmom =3600.,
> deltaTtracer=3600.,
> deltaTfreesurf= 3600.,
> deltaTClock =86400.,
> #must be same as cpl_atmSendFrq?
> abEps = 0.1,
> # ChkptFreq =15552000.,
> pChkptFreq=155520000.,
> # pChkptFreq=86400.,
> # taveFreq=3110400000.,
> # taveFreq =933120000.,
> dumpFreq=0.,
> monitorFreq =0.,
> # monitorFreq =1.,
> forcing_In_AB=.FALSE.,
> # periodicExternalForcing=.TRUE.,
> # externForcingPeriod=2592000.,
> # externForcingCycle=31104000.,
> # 2 years restoring timescale for temperature
> # tauThetaClimRelax= 62208000.,
> # 10 years restoring timescale for salinity
> # tauSaltClimRelax = 311040000.,
> # 2 months restoring timescale for temperature
> # tauThetaClimRelax= 5184000.,
> # 2 years restoring timescale for salinity
> # tauSaltClimRelax = 62208000.,
>
> &
>
> # Gridding parameters
> &PARM04
> usingCurvilinearGrid=.TRUE.,
> horizGridFile='grid_cs32',
> radius_fromHorizGrid=6370.E3,
> # delR= 10., 40., 80., 130., 190.,
> # 250., 250., 250., 250., 250.,
> # 250., 250., 250., 250., 250.,
> #^2950m
> #
> delR= 50., 70., 100., 140., 190.,
> 240., 290., 340., 390., 440.,
> 490., 540., 590., 640., 690.,
> &
>
> # Input datasets
> &PARM05
> bathyFile ='bathy_Hmin50.bin',
> hydrogThetaFile='lev_T_cs_15k.bin',
> hydrogSaltFile ='lev_S_cs_15k.bin',
> # uVelInitFile ='iniOcnU.Ex10.cpl_CTRL_CS32_15k.bin',
> # vVelInitFile ='iniOcnV.Ex10.cpl_CTRL_CS32_15k.bin',
> # pSurfInitFile ='iniOcnEta.Ex10.cpl_CTRL_CS32.bin',
> # zonalWindFile= 'J-OFURO_TauX_steady.CS32.bin',
> # meridWindFile= 'J-OFURO_TauY_steady.CS32.bin',
> # thetaClimFile= 'SST.cos_steady.CS32.bin',
> # saltClimFile= 'SSS.WOA2023_steady.CS32.bin',
> # surfQnetFile ='',
> # EmPmRFile ='',
> &
> # Coupling package parameters, OCN component:
> # useImportHFlx :: True => use the Imported HeatFlux from couler
> # useImportFW :: True => use the Imported Fresh Water flux fr cpl
> # useImportTau :: True => use the Imported Wind-Stress from couler
> # useImportSLP :: True => use the Imported Sea-level Atmos. Pressure
> # useImportRunOff :: True => use the Imported RunOff flux from coupler
> # useImportSIce :: True => use the Imported Sea-Ice mass as ice-loading
> # useImportThSIce :: True => use the Imported thSIce state var from coupler
> # useImportSltPlm :: True => use the Imported Salt-Plume flux from coupler
> # useImportFice :: True => use the Imported Seaice fraction (DIC-only)
> # useImportCO2 :: True => use the Imported atmos. CO2 from coupler
> # useImportWSpd :: True => use the Imported surf. Wind speed from coupler
> # cpl_taveFreq :: Frequency^-1 for time-Aver. output (s)
> &CPL_OCN_PARAM
> useImportHFlx =.TRUE.,
> useImportFW =.TRUE.,
> useImportTau =.TRUE.,
> useImportSLP =.TRUE.,
> useImportRunOff=.TRUE.,
> # useImportSIce =.TRUE.,
> # useImportThSIce=.TRUE.,
> # useImportSltPlm=.TRUE.,
> # useImportFice =.TRUE.,
> # useImportCO2 =.TRUE.,
> useImportWSpd =.TRUE.,
> # cpl_taveFreq=2592000.,
> # cpl_taveFreq=18000.,
> &
> # AIM physics package parameters:
>
> C-- AIM interface parameter (namelist AIM_PARAMS):
> # aim_useFMsurfBC :: select surface B.C. from Franco Molteni
> # aim_useMMsurfFc :: select Monthly Mean surface forcing (e.g., NCEP)
> # aim_surfPotTemp :: surf.Temp input file is in Pot.Temp (aim_useMMsurfFc)
> # aim_energPrecip :: account for energy of precipitation (snow & rain temp)
> # aim_splitSIOsFx :: compute separately Sea-Ice & Ocean surf. Flux
> # aim_LandFile :: file name for Land fraction [0-1]
> # aim_MMsufx :: sufix for all Monthly Mean surface forcing files
> # aim_MMsufxLength :: Length of sufix (Monthly Mean surf. forc. files)
> # aim_albFile :: file name for Albedo input file (F.M. surfBC)
> # aim_vegFile :: file name for vegetation fraction (F.M. surfBC)
> # aim_sstFile :: file name for Sea.Surf.Temp (F.M. surfBC)
> # aim_lstFile :: file name for Land.Surf.Temp (F.M. surfBC)
> # aim_oiceFile :: file name for Sea Ice fraction (F.M. surfBC)
> # aim_snowFile :: file name for Snow depth (F.M. surfBC)
> # aim_swcFile :: file name for Soil Water content (F.M. surfBC)
> # aim_dragStrato :: stratospheric-drag damping time scale (s)
> # aim_clrSkyDiag :: compute clear-sky radiation for diagnostics
> # aim_taveFreq :: Frequency^-1 for time-average output (s)
> # aim_diagFreq :: Frequency^-1 for diagnostic output (s)
> # aim_tendFreq :: Frequency^-1 for tendencies output (s)
> &AIM_PARAMS
> aim_useFMsurfBC=.TRUE.,
> aim_useMMsurfFc=.FALSE.,
> # aim_surfPotTemp=.TRUE.,
> aim_energPrecip=.TRUE.,
> aim_splitSIOsFx=.TRUE.,
> aim_LandFile= 'landFrc.cpl_FM.bin',
> aim_albFile = 'albedo_cs32.bin',
> aim_vegFile = 'vegetFrc_zeros_cs32.bin',
> # aim_vegFile = 'vegetFrc.cpl_FM.bin',
> #
> # SST.Ex10.ocn_sym.bin: Original annual average SST data(shape: (32,192))
> # aim_sstFile = 'SST.Ex10.ocn_sym.bin',
> # aim_lstFile = '',
> # aim_oiceFile= '',
> # aim_snowFile= '',
> # aim_swcFile = '',
> aim_dragStrato= 2592000.,
> aim_clrSkyDiag=.TRUE.,
> aim_diagFreq=311040000.,
> aim_select_pCO2 = 1,
> # aim_fixed_pCO2 = 350.E-6,
> aim_fixed_pCO2 = 285.E-6,
> # atmpCO2init=800.E-6,
> &
>
> #-- Physical constants
> # P0 = reference pressure [Pa=N/m2]
> # GG = gravity accel. [m/s2]
> # RD = gas constant for dry air [J/kg/K]
> # CP = specific heat at constant pressure [J/kg/K]
> # ALHC = latent heat of condensation [J/g]
> # ALHF = latent heat of freezing [J/g]
> # SBC = Stefan-Boltzmann constant
> # rainCP = heat capacity of liquid water [J/g/K]
> # tFreeze = freezing temperature of pure water [K]
>
> #-- Constants for forcing fields (namelist AIM_PAR_FOR):
> # SOLC = Solar constant (area averaged) in W/m^2
> # ALBSEA = Albedo over sea
> # ALBICE = Albedo over sea ice (for ice fraction = 1)
> # ALBSN = Albedo over snow (for snow depth > SDALB)
> # SDALB = Snow depth (mm water) corresponding to maximum albedo
> # SWCAP = Soil wetness at field capacity (volume fraction)
> # SWWIL = Soil wetness at wilting point (volume fraction)
> # hSnowWetness :: snow depth (m) corresponding to maximum wetness
> &AIM_PAR_FOR
> &
>
> #-- Constants for surface fluxes (namelist AIM_PAR_SFL):
> # FWIND0 = ratio of near-sfc wind to lowest-level wind
> # FTEMP0 = weight for near-sfc temperature extrapolation (0-1) :
> # 1 : linear extrapolation from two lowest levels
> # 0 : constant potential temperature ( = lowest level)
> # FHUM0 = weight for near-sfc specific humidity extrapolation (0-1) :
> # 1 : extrap. with constant relative hum. ( = lowest level)
> # 0 : constant specific hum. ( = lowest level)
> # CDL = drag coefficient for momentum over land
> # CDS = drag coefficient for momentum over sea
> # CHL = heat exchange coefficient over land
> # CHS = heat exchange coefficient over sea
> # VGUST = wind speed for sub-grid-scale gusts
> # CTDAY = daily-cycle correction (dTskin/dSSRad)
> # DTHETA = Potential temp. gradient for stability correction
> # dTstab = potential temp. increment for stability function derivative
> # FSTAB = Amplitude of stability correction (fraction)
> # HDRAG = Height scale for orographic correction
> # FHDRAG = Amplitude of orographic correction (fraction)
> &AIM_PAR_SFL
> # FTEMP0 = 0.,
> # FHUM0 = 0.,
> # dTstab = 0.,
> &
>
> #-- Convection constants (namelist AIM_PAR_CNV):
> # PSMIN = minimum (norm.) sfc. pressure for the occurrence of convection
> # TRCNV = time of relaxation (in hours) towards reference state
> # QBL = specific hum. threshold in the boundary layer
> # RHBL = relative hum. threshold in the boundary layer
> # RHIL = rel. hum. threshold in intermed. layers for secondary mass flux
> # ENTMAX = max. entrainment as a fraction of cloud-base mass flux
> # SMF = ratio between secondary and primary mass flux at cloud-base
> &AIM_PAR_CNV
> &
>
> #-- Constants for large-scale condensation (namelist AIM_PAR_LSC):
> # TRLSC = Relaxation time (in hours) for specific humidity
> # RHLSC = Maximum relative humidity threshold (at sigma=1)
> # DRHLSC = Vertical range of relative humidity threshold
> # QSMAX = used to define the maximum latent heat release
> &AIM_PAR_LSC
> &
>
> #-- Radiation constants (namelist AIM_PAR_RAD):
> # RHCL1 = relative hum. corresponding to cloud cover = 0
> # RHCL2 = relative hum. corresponding to cloud cover = 1
> # QACL1 = specific hum. threshold for cloud cover in the upper troposphere
> # QACL2 = specific hum. threshold for cloud cover in the upper troposphere
> # ALBCL = cloud albedo (for cloud cover = 1)
> # EPSSW = fraction of incoming solar radiation absorbed by ozone
> # EPSLW = fraction of surface LW radiation emitted directly to space
> # EMISFC = longwave surface emissivity
> #--: shortwave absorptivities (for dp = 10^5 Pa) :
> # ABSDRY = abs. of dry air (visible band)
> # ABSAER = abs. of aerosols (visible band)
> # ABSWV1 = abs. of water vapour (visible band, for dq = 1 g/kg)
> # ABSWV2 = abs. of water vapour (near IR band, for dq = 1 g/kg)
> # ABSCL1 = abs. of clouds (visible band, constant term)
> # ABSCL2 = abs. of clouds (visible band, for dw = 1 g/kg)
> #--: longwave absorptivities (per dp = 10^5 Pa) :
> # ABLWIN = abs. of air in "window" band
> # ABLCO2 = abs. of air in CO2 band
> # ABLWV1 = abs. of water vapour in H2O band 1 (weak), for dq = 1 g/kg
> # ABLWV2 = abs. of water vapour in H2O band 2 (strong), for dq = 1 g/kg
> # ABLCL1 = abs. of clouds in "window" band, constant term
> # ABLCL2 = abs. of clouds in "window" band, for dw = 1 g/kg
> &AIM_PAR_RAD
> &
>
> #-- Constants for vertical dif. and sh. conv. (namelist AIM_PAR_VDI):
> # TRSHC = relaxation time (in hours) for shallow convection
> # TRVDI = relaxation time (in hours) for moisture diffusion
> # TRVDS = relaxation time (in hours) for super-adiab. conditions
> # RHGRAD = maximum gradient of relative humidity (d_RH/d_sigma)
> # SEGRAD = minimum gradient of dry static energy (d_DSE/d_phi)
> &AIM_PAR_VDI
> &
>
> # ====================
> # | Model parameters |
> # ====================
> #
> # Continuous equation parameters
> &PARM01
> tRef=289.6, 298.1, 314.5, 335.8, 437.4,
> sRef=5*0.,
> no_slip_sides=.FALSE.,
> no_slip_bottom=.FALSE.,
> buoyancyRelation='ATMOSPHERIC',
> eosType='IDEALG',
> atm_Rq=0.6078E-3,
> integr_GeoPot=2,
> selectFindRoSurf=1,
> gravity=9.81,
> rhoConstFresh=1000.,
> implicitFreeSurface=.TRUE.,
> exactConserv=.TRUE.,
> select_rStar=2,
> nonlinFreeSurf=4,
> hFacInf=0.2,
> hFacSup=2.0,
> uniformLin_PhiSurf=.FALSE.,
> hFacMin=0.2,
> saltAdvScheme=3,
> vectorInvariantMomentum=.TRUE.,
> useJamartWetPoints=.TRUE.,
> staggerTimeStep=.TRUE.,
> readBinaryPrec=64,
> writeBinaryPrec=64,
> debugLevel=0,
> plotLevel=0,
> # globalFiles=.TRUE.,
> useSingleCpuIO=.TRUE.,
> # tempAdvScheme=77,
> # saltAdvScheme=77,
> # multiDimAdvection=.TRUE.,
> &
>
> # Elliptic solver parameters
> &PARM02
> cg2dMaxIters=200,
> #cg2dTargetResidual=1.E-12,
> cg2dTargetResWunit=1.E-17,
> &
>
> # Time stepping parameters
> &PARM03
> #200yrs
> # nIter0=1800000,
> startTime=0.0,
> endTime=155520000.0,
> deltaT =450.,
> # deltaTmom =450.,
> # deltaTtracer=450.,
> # deltaTfreesurf=450.,
> # deltaTClock=43200., ????????????
> # must be same as cpl_atmSendFrq?
> abEps=0.1,
> forcing_In_AB=.FALSE.,
> cAdjFreq=0.,
> # ChkptFreq =31104000.,
> pChkptFreq=155520000.,
> # pChkptFreq=86400.,
> # taveFreq=311040000.,
> dumpFreq=0.,
> monitorFreq=0.,
> # monitorFreq=1.,
> # pickupStrictlyMatch=.FALSE.,
> &
>
> # Gridding parameters
> &PARM04
> usingCurvilinearGrid=.TRUE.,
> horizGridFile='grid_cs32',
> radius_fromHorizGrid=6370.E3,
> delR=100.E2, 250.E2, 300.E2, 200.E2, 150.E2,
> &
>
> # Input datasets
> &PARM05
> topoFile='topo.cpl_FM.bin',
> # hydrogThetaFile='iniAtmT.ERA2024_CS32_5k.bin',
> # hydrogSaltFile ='iniAtmS.ERA2024_CS32_5k.bin',
> # uVelInitFile='iniAtmU.Ex10.cpl_CTRL_CS32_5k.bin',
> # vVelInitFile='iniAtmV.Ex10.cpl_CTRL_CS32_5k.bin',
> # pSurfInitFile='iniAtmEta.Ex10.cpl_CTRL_CS32.bin',
>
>
> &
> # ATM+OCN Coupler package parameters:
> # cpl_sequential :: =0/1 : selects Synchronous/Sequential Coupling
> #-- cpl_exchange_[xxx] :: controls exchange of [xxx] fields with OCN & ATM components
> # :: =0 : none ; =1,3 : exch [xxx] fields with OCN comp.
> # :: =2,3 : exch [xxx] fileds with ATM comp.
> # cpl_exchange_RunOff :: controls exchange of RunOff fields
> # cpl_exchange1W_sIce :: controls 1-way exchange of seaice (step fwd in ATM)
> # cpl_exchange2W_sIce :: controls 2-way exchange of ThSIce variables
> # cpl_exchange_SaltPl :: controls exchange of Salt-Plume fields
> # cpl_exchange_DIC :: controls exchange of DIC variables
> #--
> # runOffMapSize :: Nb of connection in the runOff map (nROmap)
> # runOffMapFile :: File name for the runOff map (max_length= 80c)
> #-- IMPORTANT: Unlike other MITgcm parameter files, this namelist is read-in
> # directly without prior removing of commented lines (i.e., starting with "#")
> #-------------
> &COUPLER_PARAMS
> runOffMapSize = 3644,
> runOffMapFile = 'runOff_cs32_3644.bin',
> cpl_sequential = 0,
> cpl_exchange_RunOff = 3,
> cpl_exchange1W_sIce = 3,
> cpl_exchange2W_sIce = 0,
> cpl_exchange_SaltPl = 0,
> cpl_exchange_DIC = 3,
> /
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
A non-text attachment was scrubbed...
Name: fig_sst_01.jpg
Type: image/jpeg
Size: 1090035 bytes
Desc: not available
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20260120/0724637e/attachment-0001.jpg>
More information about the MITgcm-support
mailing list