/[MITgcm]/MITgcm/pkg/thsice/thsice_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_readparms.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.13 by jscott, Mon Apr 9 17:44:13 2007 UTC revision 1.14 by jmc, Sun Apr 29 23:48:44 2007 UTC
# Line 47  C     iUnit     :: Work variable for IO Line 47  C     iUnit     :: Work variable for IO
47  C--   Th-Sea-ICE parameter  C--   Th-Sea-ICE parameter
48        NAMELIST /THSICE_CONST/        NAMELIST /THSICE_CONST/
49       &  rhos, rhoi, rhosw, rhofw,       &  rhos, rhoi, rhosw, rhofw,
50       &  cpice, cpwater,       &  cpIce, cpWater,
51       &  kice, ksnow,       &  kIce, kSnow,
52       &  transcoef, Lfresh, qsnow,       &  bMeltCoef, Lfresh, qsnow,
53       &  albColdSnow, albWarmSnow, tempSnowAlb,       &  albColdSnow, albWarmSnow, tempSnowAlb,
54       &  albOldSnow, hNewSnowAge, snowAgTime,       &  albOldSnow, hNewSnowAge, snowAgTime,
55       &  albIceMax, albIceMin, hAlbIce, hAlbSnow,       &  albIceMax, albIceMin, hAlbIce, hAlbSnow,
56       &  i0, ksolar,       &  i0swFrac, ksolar, dhSnowLin,
57       &  saltice, S_winton, mu_Tf,       &  saltIce, S_winton, mu_Tf,
58       &  Tf0kel,       &  Tf0kel, Terrmax, nitMaxTsf,
59       &  himin, Terrmax, nitMaxTsf, hiMax, hsMax,       &  hIceMin, hiMax, hsMax, iceMaskMax, iceMaskMin,
60       &  iceMaskmax, iceMaskmin, himin0,       &  fracEnMelt, fracEnFreez, hThinIce, hThickIce, hNewIceMax
      &  frac_energy, hihig  
61    
62        NAMELIST /THSICE_PARM01/        NAMELIST /THSICE_PARM01/
63       &     startIceModel, stepFwd_oceMxL,       &     startIceModel, stepFwd_oceMxL,
# Line 89  C--   Default values (constants) Line 88  C--   Default values (constants)
88        rhoi     = 900. _d 0        rhoi     = 900. _d 0
89        rhosw    = rhoConst        rhosw    = rhoConst
90        rhofw    = rhoConstFresh        rhofw    = rhoConstFresh
91        cpice    = 2106. _d 0        cpIce    = 2106. _d 0
92        cpwater  = HeatCapacity_Cp        cpWater  = HeatCapacity_Cp
93        kice     = 2.03 _d 0        kIce     = 2.03 _d 0
94        ksnow    = 0.30 _d 0        kSnow    = 0.30 _d 0
95        transcoef=0.006 _d 0        bMeltCoef=0.006 _d 0
96        Lfresh   = 3.34 _d 5        Lfresh   = 3.34 _d 5
97        qsnow    = Lfresh        qsnow    = Lfresh
98        albColdSnow= 0.85 _d 0        albColdSnow= 0.85 _d 0
# Line 106  C--   Default values (constants) Line 105  C--   Default values (constants)
105        hAlbSnow   = 0.30 _d 0        hAlbSnow   = 0.30 _d 0
106        hNewSnowAge= 2. _d -3        hNewSnowAge= 2. _d -3
107        snowAgTime = 50. _d 0 * 86400. _d 0        snowAgTime = 50. _d 0 * 86400. _d 0
108        i0       = 0.3 _d 0        i0swFrac = 0.3 _d 0
109        ksolar   = 1.5 _d 0        ksolar   = 1.5 _d 0
110        saltice  = 4. _d 0        dhSnowLin= 0. _d 0
111          saltIce  = 4. _d 0
112        S_winton = 1. _d 0        S_winton = 1. _d 0
113        mu_Tf    = 0.054 _d 0        mu_Tf    = 0.054 _d 0
114        Tf0kel   = celsius2K        Tf0kel   = celsius2K
       himin    = 0.01 _d 0  
115        Terrmax  = 5.0 _d -1        Terrmax  = 5.0 _d -1
116        nitMaxTsf= 20        nitMaxTsf= 20
117          hIceMin    = 1. _d -2
118        hiMax      = 10. _d 0        hiMax      = 10. _d 0
119        hsMax      = 10. _d 0        hsMax      = 10. _d 0
120        iceMaskmax = 1.  _d 0        iceMaskMax =  1. _d 0
121        iceMaskmin =  .1 _d 0        iceMaskMin = 0.1 _d 0
122        himin0     = 0.2 _d 0        fracEnMelt = 0.4 _d 0
123        frac_energy=  .4 _d 0        fracEnFreez=  0. _d 0
124        hihig      = 2.5 _d 0        hThinIce   = 0.2 _d 0
125          hThickIce  = 2.5 _d 0
126          hNewIceMax = UNSET_RL
127    
128  C--   Default values (parameters)  C--   Default values (parameters)
129        stepFwd_oceMxL  = .FALSE.        stepFwd_oceMxL  = .FALSE.
# Line 175  C--   Read parameters from open data fil Line 177  C--   Read parameters from open data fil
177  C--   Close the open data file  C--   Close the open data file
178        CLOSE(iUnit)        CLOSE(iUnit)
179    
180    C-    neutral default:
181          IF ( hNewIceMax .EQ. UNSET_RL ) hNewIceMax = hiMax
182    
183  C     If using the same time step for both icetop temp solver  C     If using the same time step for both icetop temp solver
184  C     and ice thickness/growth, use thSIce_deltaT value  C     and ice thickness/growth, use thSIce_deltaT value
185        IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT        IF ( thSIce_dtTemp .EQ. UNSET_RL ) thSIce_dtTemp=thSIce_deltaT
# Line 241  C---+----1----+----2----+----3----+----4 Line 246  C---+----1----+----2----+----3----+----4
246        iUnit = standardMessageUnit        iUnit = standardMessageUnit
247  c     CALL MDSFINDUNIT( iUnit, mythid )  c     CALL MDSFINDUNIT( iUnit, mythid )
248  c     OPEN(iUnit,file='thsice_check_params',status='unknown')  c     OPEN(iUnit,file='thsice_check_params',status='unknown')
249        WRITE(iUnit,*) 'ThSI: rhos    =',rhos        WRITE(iUnit,*) 'ThSI: rhos      =',rhos
250        WRITE(iUnit,*) 'ThSI: rhoi    =',rhoi        WRITE(iUnit,*) 'ThSI: rhoi      =',rhoi
251        WRITE(iUnit,*) 'ThSI: rhosw   =',rhosw        WRITE(iUnit,*) 'ThSI: rhosw     =',rhosw
252        WRITE(iUnit,*) 'ThSI: rhofw   =',rhofw        WRITE(iUnit,*) 'ThSI: rhofw     =',rhofw
253        WRITE(iUnit,*) 'ThSI: rhoiw   =',rhoiw        WRITE(iUnit,*) 'ThSI: rhoiw     =',rhoiw
254        WRITE(iUnit,*) 'ThSI: cpice   =',cpice        WRITE(iUnit,*) 'ThSI: cpIce     =',cpIce
255        WRITE(iUnit,*) 'ThSI: cpwater =',cpwater        WRITE(iUnit,*) 'ThSI: cpWater   =',cpWater
256        WRITE(iUnit,*) 'ThSI: kice    =',kice        WRITE(iUnit,*) 'ThSI: kIce      =',kIce
257        WRITE(iUnit,*) 'ThSI: ksnow   =',ksnow        WRITE(iUnit,*) 'ThSI: kSnow     =',kSnow
258        WRITE(iUnit,*) 'ThSI: transcoef=',transcoef        WRITE(iUnit,*) 'ThSI: bMeltCoef =',bMeltCoef
259        WRITE(iUnit,*) 'ThSI: Lfresh  =',Lfresh        WRITE(iUnit,*) 'ThSI: Lfresh    =',Lfresh
260        WRITE(iUnit,*) 'ThSI: qsnow   =',qsnow        WRITE(iUnit,*) 'ThSI: qsnow     =',qsnow
261        WRITE(iUnit,*) 'ThSI: albColdSnow=',albColdSnow        WRITE(iUnit,*) 'ThSI: albColdSnow=',albColdSnow
262        WRITE(iUnit,*) 'ThSI: albWarmSnow=',albWarmSnow        WRITE(iUnit,*) 'ThSI: albWarmSnow=',albWarmSnow
263        WRITE(iUnit,*) 'ThSI: tempSnowAlb=',tempSnowAlb        WRITE(iUnit,*) 'ThSI: tempSnowAlb=',tempSnowAlb
# Line 263  c     OPEN(iUnit,file='thsice_check_para Line 268  c     OPEN(iUnit,file='thsice_check_para
268        WRITE(iUnit,*) 'ThSI: albIceMin =',albIceMin        WRITE(iUnit,*) 'ThSI: albIceMin =',albIceMin
269        WRITE(iUnit,*) 'ThSI: hAlbIce   =',hAlbIce        WRITE(iUnit,*) 'ThSI: hAlbIce   =',hAlbIce
270        WRITE(iUnit,*) 'ThSI: hAlbSnow  =',hAlbSnow        WRITE(iUnit,*) 'ThSI: hAlbSnow  =',hAlbSnow
271        WRITE(iUnit,*) 'ThSI: i0      =',i0        WRITE(iUnit,*) 'ThSI: i0swFrac  =',i0swFrac
272        WRITE(iUnit,*) 'ThSI: ksolar  =',ksolar        WRITE(iUnit,*) 'ThSI: ksolar    =',ksolar
273        WRITE(iUnit,*) 'ThSI: saltice =',saltice        WRITE(iUnit,*) 'ThSI: dhSnowLin =',dhSnowLin
274        WRITE(iUnit,*) 'ThSI: S_winton=',S_winton        WRITE(iUnit,*) 'ThSI: saltIce   =',saltIce
275        WRITE(iUnit,*) 'ThSI: mu_Tf   =',mu_Tf        WRITE(iUnit,*) 'ThSI: S_winton  =',S_winton
276        WRITE(iUnit,*) 'ThSI: Tf0kel  =',Tf0kel        WRITE(iUnit,*) 'ThSI: mu_Tf     =',mu_Tf
277        WRITE(iUnit,*) 'ThSI: Tmlt1   =',Tmlt1        WRITE(iUnit,*) 'ThSI: Tf0kel    =',Tf0kel
278        WRITE(iUnit,*) 'ThSI: himin   =',himin        WRITE(iUnit,*) 'ThSI: Tmlt1     =',Tmlt1
279        WRITE(iUnit,*) 'ThSI: Terrmax =',Terrmax        WRITE(iUnit,*) 'ThSI: Terrmax   =',Terrmax
280        WRITE(iUnit,*) 'ThSI: nitMaxTsf=',nitMaxTsf        WRITE(iUnit,*) 'ThSI: nitMaxTsf =',nitMaxTsf
281        WRITE(iUnit,*) 'ThSI: hiMax   =',hiMax        WRITE(iUnit,*) 'ThSI: hIceMin   =',hIceMin
282        WRITE(iUnit,*) 'ThSI: hsMax   =',hsMax        WRITE(iUnit,*) 'ThSI: hiMax     =',hiMax
283        WRITE(iUnit,*) 'ThSI: iceMaskmax=',iceMaskmax        WRITE(iUnit,*) 'ThSI: hsMax     =',hsMax
284        WRITE(iUnit,*) 'ThSI: iceMaskmin=',iceMaskmin        WRITE(iUnit,*) 'ThSI: iceMaskMax =',iceMaskMax
285        WRITE(iUnit,*) 'ThSI: himin0  =',himin0        WRITE(iUnit,*) 'ThSI: iceMaskMin =',iceMaskMin
286        WRITE(iUnit,*) 'ThSI: frac_energy',frac_energy        WRITE(iUnit,*) 'ThSI: fracEnMelt =',fracEnMelt
287        WRITE(iUnit,*) 'ThSI: hihig   =',hihig        WRITE(iUnit,*) 'ThSI: fracEnFreez=',fracEnFreez
288          WRITE(iUnit,*) 'ThSI: hThinIce   =',hThinIce
289          WRITE(iUnit,*) 'ThSI: hThickIce  =',hThickIce
290          WRITE(iUnit,*) 'ThSI: hNewIceMax =',hNewIceMax
291        WRITE(iUnit,*) 'ThSI: stressReduction =',stressReduction        WRITE(iUnit,*) 'ThSI: stressReduction =',stressReduction
292        WRITE(iUnit,*) 'ThSI: thSIceAdvScheme =',thSIceAdvScheme        WRITE(iUnit,*) 'ThSI: thSIceAdvScheme =',thSIceAdvScheme
293        WRITE(iUnit,*) 'ThSI: thSIce_diffK  =',thSIce_diffK        WRITE(iUnit,*) 'ThSI: thSIce_diffK  =',thSIce_diffK
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22