/[MITgcm]/MITgcm/pkg/exf/exf_summary.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_summary.F

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

revision 1.15 by jmc, Mon May 14 19:34:57 2007 UTC revision 1.20 by jmc, Fri Jun 26 18:02:15 2009 UTC
# Line 97  C--   Print general parameters: Line 97  C--   Print general parameters:
97        CALL WRITE_0D_L( useExfCheckRange, INDEX_NONE,        CALL WRITE_0D_L( useExfCheckRange, INDEX_NONE,
98       &  'useExfCheckRange =',       &  'useExfCheckRange =',
99       &  ' /* check for fields range */')       &  ' /* check for fields range */')
100        CALL WRITE_0D_R8( exf_monFreq, INDEX_NONE, 'exf_monFreq =',        CALL WRITE_0D_RL( exf_monFreq, INDEX_NONE, 'exf_monFreq =',
101       &  ' /* EXF monitor frequency [ s ] */')       &  ' /* EXF monitor frequency [ s ] */')
102        CALL WRITE_0D_R8( repeatPeriod, INDEX_NONE, 'repeatPeriod =',        CALL WRITE_0D_RL( repeatPeriod, INDEX_NONE, 'repeatPeriod =',
103       &  ' /* period for cycling forcing dataset [ s ] */')       &  ' /* period for cycling forcing dataset [ s ] */')
104        CALL WRITE_0D_R8( climtempfreeze, INDEX_NONE,'climTempFreeze=',        CALL WRITE_0D_RL( climtempfreeze, INDEX_NONE,'climTempFreeze=',
105       &  ' /* Minimum climatological temperature [deg.C] */')       &  ' /* Minimum climatological temperature [deg.C] */')
106        CALL WRITE_0D_R8( windStressMax, INDEX_NONE,'windStressMax =',        CALL WRITE_0D_RL( windStressMax, INDEX_NONE,'windStressMax =',
107       &  ' /* Maximum absolute windstress [ Pa ] */')       &  ' /* Maximum absolute windstress [ Pa ] */')
108        CALL WRITE_0D_L( stressIsOnCgrid,INDEX_NONE,'stressIsOnCgrid =',        CALL WRITE_0D_L( stressIsOnCgrid,INDEX_NONE,'stressIsOnCgrid =',
109       &  ' /* set u,v_stress on Arakawa C-grid */')       &  ' /* set u,v_stress on Arakawa C-grid */')
110        CALL WRITE_0D_R8( cen2kel, INDEX_NONE, 'cen2kel =',        CALL WRITE_0D_RL( cen2kel, INDEX_NONE, 'cen2kel =',
111       &  ' /* conversion of deg. Centigrade to Kelvin [K] */')       &  ' /* conversion of deg. Centigrade to Kelvin [K] */')
112        CALL WRITE_0D_R8( gravity_mks, INDEX_NONE, 'gravity_mks=',        CALL WRITE_0D_RL( gravity_mks, INDEX_NONE, 'gravity_mks=',
113       &  ' /* gravitational acceleration [m/s^2] */')       &  ' /* gravitational acceleration [m/s^2] */')
114        CALL WRITE_0D_R8( atmrho, INDEX_NONE, 'atmrho =',        CALL WRITE_0D_RL( atmrho, INDEX_NONE, 'atmrho =',
115       & '  /* mean atmospheric density [kg/m^3] */')       & '  /* mean atmospheric density [kg/m^3] */')
116        CALL WRITE_0D_R8( atmcp, INDEX_NONE, 'atmcp =',        CALL WRITE_0D_RL( atmcp, INDEX_NONE, 'atmcp =',
117       & '  /* mean atmospheric specific heat [J/kg/K] */')       & '  /* mean atmospheric specific heat [J/kg/K] */')
118        CALL WRITE_0D_R8( flamb, INDEX_NONE, 'flamb =',        CALL WRITE_0D_RL( flamb, INDEX_NONE, 'flamb =',
119       & '  /* latent heat of evaporation [J/kg] */')       & '  /* latent heat of evaporation [J/kg] */')
120        CALL WRITE_0D_R8( flami, INDEX_NONE, 'flami =',        CALL WRITE_0D_RL( flami, INDEX_NONE, 'flami =',
121       & '  /* latent heat of pure-ice melting [J/kg] */')       & '  /* latent heat of pure-ice melting [J/kg] */')
122        CALL WRITE_0D_R8( cvapor_fac, INDEX_NONE, 'cvapor_fac =',        CALL WRITE_0D_RL( cvapor_fac, INDEX_NONE, 'cvapor_fac =',
123       &  ' /* const. for Saturation calculation [?] */')       &  ' /* const. for Saturation calculation [?] */')
124        CALL WRITE_0D_R8( cvapor_exp, INDEX_NONE, 'cvapor_exp =',        CALL WRITE_0D_RL( cvapor_exp, INDEX_NONE, 'cvapor_exp =',
125       &  ' /* const. for Saturation calculation [?] */')       &  ' /* const. for Saturation calculation [?] */')
126        CALL WRITE_0D_R8( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',        CALL WRITE_0D_RL( cvapor_fac_ice, INDEX_NONE, 'cvapor_fac_ice=',
127       &  ' /* const. for Saturation calculation [?] */')       &  ' /* const. for Saturation calculation [?] */')
128        CALL WRITE_0D_R8( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',        CALL WRITE_0D_RL( cvapor_exp_ice, INDEX_NONE, 'cvapor_exp_ice=',
129       &  ' /* const. for Saturation calculation [?] */')       &  ' /* const. for Saturation calculation [?] */')
130        CALL WRITE_0D_R8( humid_fac, INDEX_NONE, 'humid_fac =',        CALL WRITE_0D_RL( humid_fac, INDEX_NONE, 'humid_fac =',
131       &  ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')       &  ' /* humidity coef. in virtual temp. [(kg/kg)^-1] */')
132        CALL WRITE_0D_R8( gamma_blk, INDEX_NONE, 'gamma_blk =',        CALL WRITE_0D_RL( gamma_blk, INDEX_NONE, 'gamma_blk =',
133       &  ' /* adiabatic lapse rate [?] */')       &  ' /* adiabatic lapse rate [?] */')
134        CALL WRITE_0D_R8( saltsat, INDEX_NONE, 'saltsat =',        CALL WRITE_0D_RL( saltsat, INDEX_NONE, 'saltsat =',
135       &  ' /* reduction of Qsat over salty water [-] */')       &  ' /* reduction of Qsat over salty water [-] */')
136        CALL WRITE_0D_R8( cdrag_1, INDEX_NONE, 'cDrag_1 =',        CALL WRITE_0D_L( noNegativeEvap,INDEX_NONE,'noNegativeEvap =',
137         &  ' /* prevent negative Evaporation */')
138          CALL WRITE_0D_RL( sstExtrapol, INDEX_NONE, 'sstExtrapol =',
139         &  ' /* extrapolation coeff from lev. 1 & 2 to surf [-] */')
140          CALL WRITE_0D_RL( cdrag_1, INDEX_NONE, 'cDrag_1 =',
141       &  ' /* coef used in drag calculation [?] */')       &  ' /* coef used in drag calculation [?] */')
142        CALL WRITE_0D_R8( cdrag_2, INDEX_NONE, 'cDrag_2 =',        CALL WRITE_0D_RL( cdrag_2, INDEX_NONE, 'cDrag_2 =',
143       &  ' /* coef used in drag calculation [?] */')       &  ' /* coef used in drag calculation [?] */')
144        CALL WRITE_0D_R8( cdrag_3, INDEX_NONE, 'cDrag_3 =',        CALL WRITE_0D_RL( cdrag_3, INDEX_NONE, 'cDrag_3 =',
145       &  ' /* coef used in drag calculation [?] */')       &  ' /* coef used in drag calculation [?] */')
146        CALL WRITE_0D_R8( cstanton_1, INDEX_NONE, 'cStanton_1 =',        CALL WRITE_0D_RL( cstanton_1, INDEX_NONE, 'cStanton_1 =',
147       &  ' /* coef used in Stanton number calculation [?] */')       &  ' /* coef used in Stanton number calculation [?] */')
148        CALL WRITE_0D_R8( cstanton_2, INDEX_NONE, 'cStanton_2 =',        CALL WRITE_0D_RL( cstanton_2, INDEX_NONE, 'cStanton_2 =',
149       &  ' /* coef used in Stanton number calculation [?] */')       &  ' /* coef used in Stanton number calculation [?] */')
150        CALL WRITE_0D_R8( cdalton, INDEX_NONE, 'cDalton =',        CALL WRITE_0D_RL( cdalton, INDEX_NONE, 'cDalton =',
151       &  ' /* coef used in Dalton number calculation [?] */')       &  ' /* coef used in Dalton number calculation [?] */')
152        CALL WRITE_0D_R8( exf_scal_BulkCdn, INDEX_NONE,        CALL WRITE_0D_RL( exf_scal_BulkCdn, INDEX_NONE,
153       &  'exf_scal_BulkCdn=',       &  'exf_scal_BulkCdn=',
154       &  ' /* Drag coefficient scaling factor [-] */')       &  ' /* Drag coefficient scaling factor [-] */')
155        CALL WRITE_0D_R8( zolmin, INDEX_NONE, 'zolmin =',        CALL WRITE_0D_RL( zolmin, INDEX_NONE, 'zolmin =',
156       &  ' /* minimum stability parameter [?] */')       &  ' /* minimum stability parameter [?] */')
157        CALL WRITE_0D_R8( psim_fac, INDEX_NONE, 'psim_fac =',        CALL WRITE_0D_RL( psim_fac, INDEX_NONE, 'psim_fac =',
158       &  ' /* coef used in turbulent fluxes calculation [-] */')       &  ' /* coef used in turbulent fluxes calculation [-] */')
159        CALL WRITE_0D_R8( zref, INDEX_NONE, 'zref =',        CALL WRITE_0D_RL( zref, INDEX_NONE, 'zref =',
160       & '  /* reference height [ m ] */')       & '  /* reference height [ m ] */')
161        CALL WRITE_0D_R8( hu, INDEX_NONE, 'hu =',        CALL WRITE_0D_RL( hu, INDEX_NONE, 'hu =',
162       & '  /* height of mean wind [ m ] */')       & '  /* height of mean wind [ m ] */')
163        CALL WRITE_0D_R8( ht, INDEX_NONE, 'ht =',        CALL WRITE_0D_RL( ht, INDEX_NONE, 'ht =',
164       & '  /* height of mean temperature [ m ] */')       & '  /* height of mean temperature [ m ] */')
165        CALL WRITE_0D_R8( hq, INDEX_NONE, 'hq =',        CALL WRITE_0D_RL( hq, INDEX_NONE, 'hq =',
166       & '  /* height of mean spec.humidity [ m ] */')       & '  /* height of mean spec.humidity [ m ] */')
167        CALL WRITE_0D_R8( umin, INDEX_NONE, 'uMin =',        CALL WRITE_0D_RL( umin, INDEX_NONE, 'uMin =',
168       &  ' /* minimum wind speed [m/s] */')       &  ' /* minimum wind speed [m/s] */')
169        CALL WRITE_0D_L( useStabilityFct_overIce, INDEX_NONE,        CALL WRITE_0D_L( useStabilityFct_overIce, INDEX_NONE,
170       &  'useStabilityFct_overIce=',       &  'useStabilityFct_overIce=',
171       &  ' /* transfert Coeffs over sea-ice depend on stability */')       &  ' /* transfert Coeffs over sea-ice depend on stability */')
172        CALL WRITE_0D_R8( exf_iceCd, INDEX_NONE, 'exf_iceCd =',        CALL WRITE_0D_RL( exf_iceCd, INDEX_NONE, 'exf_iceCd =',
173       &  ' /* drag coefficient over sea-ice (fixed) [-] */')       &  ' /* drag coefficient over sea-ice (fixed) [-] */')
174        CALL WRITE_0D_R8( exf_iceCe, INDEX_NONE, 'exf_iceCe =',        CALL WRITE_0D_RL( exf_iceCe, INDEX_NONE, 'exf_iceCe =',
175       &  ' /* transfert coeff. over sea-ice, for Evap (fixed) [-] */')       &  ' /* transfert coeff. over sea-ice, for Evap (fixed) [-] */')
176        CALL WRITE_0D_R8( exf_iceCh, INDEX_NONE, 'exf_iceCh =',        CALL WRITE_0D_RL( exf_iceCh, INDEX_NONE, 'exf_iceCh =',
177       &  ' /* transfert coeff. over sea-ice, Sens.Heat.(fixed)[-] */')       &  ' /* transfert coeff. over sea-ice, Sens.Heat.(fixed)[-] */')
178        CALL WRITE_0D_R8( exf_albedo, INDEX_NONE, 'exf_albedo =',        CALL WRITE_0D_RL( exf_albedo, INDEX_NONE, 'exf_albedo =',
179       & '  /* Sea-water albedo [-] */')       & '  /* Sea-water albedo [-] */')
180        CALL WRITE_0D_R8( ocean_emissivity, INDEX_NONE,        CALL WRITE_0D_RL( ocean_emissivity, INDEX_NONE,
181       &  'ocean_emissivity =',       &  'ocean_emissivity =',
182       &  ' /* longwave ocean-surface emissivity [-] */')       &  ' /* longwave ocean-surface emissivity [-] */')
183        CALL WRITE_0D_R8( ice_emissivity, INDEX_NONE,'ice_emissivity =',        CALL WRITE_0D_RL( ice_emissivity, INDEX_NONE,'ice_emissivity =',
184       &  ' /* longwave seaice emissivity [-] */')       &  ' /* longwave seaice emissivity [-] */')
185        CALL WRITE_0D_R8(snow_emissivity, INDEX_NONE,'snow_emissivity =',        CALL WRITE_0D_RL(snow_emissivity, INDEX_NONE,'snow_emissivity =',
186       &  ' /* longwave snow  emissivity [-] */')       &  ' /* longwave snow  emissivity [-] */')
187        WRITE(msgBuf,'(A)') ' '        WRITE(msgBuf,'(A)') ' '
188        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
# Line 566  C--   Runoff. Line 570  C--   Runoff.
570       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
571        il = ilnblnk(runofffile)        il = ilnblnk(runofffile)
572        write(msgbuf,'(a,f12.0)')        write(msgbuf,'(a,f12.0)')
573       &'   Runnoff starts at     ',       &'   Runoff starts at     ',
574       &    runoffstartdate       &    runoffstartdate
575        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
576       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
# Line 590  C--   Runoff. Line 594  C--   Runoff.
594       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
595  #endif /* ALLOW_RUNOFF */  #endif /* ALLOW_RUNOFF */
596    
597  #ifdef DOWNWARD_RADIATION  #ifdef ALLOW_DOWNWARD_RADIATION
598  C--   Downward shortwave.  C--   Downward shortwave.
599        il = ilnblnk(swdownfile)        il = ilnblnk(swdownfile)
600        write(msgbuf,'(a)') ' '        write(msgbuf,'(a)') ' '
# Line 665  C--   Atmospheric pressure. Line 669  C--   Atmospheric pressure.
669        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
670       &                    SQUEEZE_RIGHT , mythid)       &                    SQUEEZE_RIGHT , mythid)
671  #endif  #endif
672    
673    #ifdef ALLOW_ICE_AREAMASK
674    C--   fractional ice-covered area MASK.
675          il = ilnblnk(areamaskfile)
676          write(msgbuf,'(a)') ' '
677          call print_message( msgbuf, standardmessageunit,
678         &                    SQUEEZE_RIGHT , mythid)
679          write(msgbuf,'(a,f12.0)')
680         &'   fractional ice-covered area MASK starts at      ',
681         &    areamaskstartdate
682          call print_message( msgbuf, standardmessageunit,
683         &                    SQUEEZE_RIGHT , mythid)
684          write(msgbuf,'(a,f12.0)')
685         &'   fractional ice-covered area MASK period is      ',
686         &    areamaskperiod
687          call print_message( msgbuf, standardmessageunit,
688         &                    SQUEEZE_RIGHT , mythid)
689          write(msgbuf,'(a)')
690         &'   fractional ice-covered area MASK is read from file:'
691          call print_message( msgbuf, standardmessageunit,
692         &                    SQUEEZE_RIGHT , mythid)
693          write(msgbuf,'(a,a,a)')
694         &'   >>  ',areamaskfile(1:il),'  <<'
695          call print_message( msgbuf, standardmessageunit,
696         &                    SQUEEZE_RIGHT , mythid)
697    #endif
698    
699        write(msgbuf,'(a)') ' '        write(msgbuf,'(a)') ' '
700        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.20
  ViewVC Help
Powered by ViewVC 1.1.22