Changeset 9169 for branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

Timestamp:

2017-12-26T17:32:56+01:00 (6 years ago)

Author:

gm

Message:

dev_merge_2017: all SRC: finalize the removal of useless warning when reading namelist_cfg + remove all nn_closea + nn_msh replaced by a logical

File:

• branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90 (modified) (3 diffs)

branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/DOM/phycst.F90

@@ -113 +113 @@
       ! ------------------------
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '       Constants'
+      IF(lwp) WRITE(numout,*) '   Constants'
 
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          mathematical constant                 rpi = ', rpi
+      IF(lwp) WRITE(numout,*) '      mathematical constant                 rpi = ', rpi
 
       rsiyea = 365.25_wp * rday * 2._wp * rpi / 6.283076_wp
@@ -126 +126 @@
 #endif
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          day                                rday   = ', rday,   ' s'
-      IF(lwp) WRITE(numout,*) '          sideral year                       rsiyea = ', rsiyea, ' s'
-      IF(lwp) WRITE(numout,*) '          sideral day                        rsiday = ', rsiday, ' s'
-      IF(lwp) WRITE(numout,*) '          omega                              omega  = ', omega,  ' s^-1'
-
+      IF(lwp) WRITE(numout,*) '      day                                rday   = ', rday,   ' s'
+      IF(lwp) WRITE(numout,*) '      sideral year                       rsiyea = ', rsiyea, ' s'
+      IF(lwp) WRITE(numout,*) '      sideral day                        rsiday = ', rsiday, ' s'
+      IF(lwp) WRITE(numout,*) '      omega                              omega  = ', omega,  ' s^-1'
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          nb of months per year               raamo = ', raamo, ' months'
-      IF(lwp) WRITE(numout,*) '          nb of hours per day                 rjjhh = ', rjjhh, ' hours'
-      IF(lwp) WRITE(numout,*) '          nb of minutes per hour              rhhmm = ', rhhmm, ' mn'
-      IF(lwp) WRITE(numout,*) '          nb of seconds per minute            rmmss = ', rmmss, ' s'
-
+      IF(lwp) WRITE(numout,*) '      nb of months per year               raamo = ', raamo, ' months'
+      IF(lwp) WRITE(numout,*) '      nb of hours per day                 rjjhh = ', rjjhh, ' hours'
+      IF(lwp) WRITE(numout,*) '      nb of minutes per hour              rhhmm = ', rhhmm, ' mn'
+      IF(lwp) WRITE(numout,*) '      nb of seconds per minute            rmmss = ', rmmss, ' s'
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          earth radius                         ra   = ', ra, ' m'
-      IF(lwp) WRITE(numout,*) '          gravity                              grav = ', grav , ' m/s^2'
-
+      IF(lwp) WRITE(numout,*) '      earth radius                         ra   = ', ra, ' m'
+      IF(lwp) WRITE(numout,*) '      gravity                              grav = ', grav , ' m/s^2'
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          triple point of temperature      rtt      = ', rtt     , ' K'
-      IF(lwp) WRITE(numout,*) '          freezing point of water          rt0      = ', rt0     , ' K'
-      IF(lwp) WRITE(numout,*) '          melting point of snow            rt0_snow = ', rt0_snow, ' K'
-      IF(lwp) WRITE(numout,*) '          melting point of ice             rt0_ice  = ', rt0_ice , ' K'
-
-      IF(lwp) WRITE(numout,*) '          reference density and heat capacity now defined in eosbn2.f90'
+      IF(lwp) WRITE(numout,*) '      triple point of temperature      rtt      = ', rtt     , ' K'
+      IF(lwp) WRITE(numout,*) '      freezing point of water          rt0      = ', rt0     , ' K'
+      IF(lwp) WRITE(numout,*) '      melting point of snow            rt0_snow = ', rt0_snow, ' K'
+      IF(lwp) WRITE(numout,*) '      melting point of ice             rt0_ice  = ', rt0_ice , ' K'
+      IF(lwp) WRITE(numout,*)
+      IF(lwp) WRITE(numout,*) '   reference density and heat capacity now defined in eosbn2.f90'
               
 #if defined key_lim3 || defined key_cice
@@ -163 +160 @@
          WRITE(numout,*)
 #if defined key_cice
-         WRITE(numout,*) '          thermal conductivity of the snow          = ', rcdsn   , ' J/s/m/K'
+         WRITE(numout,*) '      thermal conductivity of the snow          = ', rcdsn   , ' J/s/m/K'
 #endif
-         WRITE(numout,*) '          thermal conductivity of pure ice          = ', rcdic   , ' J/s/m/K'
-         WRITE(numout,*) '          fresh ice specific heat                   = ', cpic    , ' J/kg/K'
-         WRITE(numout,*) '          latent heat of fusion of fresh ice / snow = ', lfus    , ' J/kg'
+         WRITE(numout,*) '      thermal conductivity of pure ice          = ', rcdic   , ' J/s/m/K'
+         WRITE(numout,*) '      fresh ice specific heat                   = ', cpic    , ' J/kg/K'
+         WRITE(numout,*) '      latent heat of fusion of fresh ice / snow = ', lfus    , ' J/kg'
 #if defined key_lim3 || defined key_cice
-         WRITE(numout,*) '          latent heat of subl.  of fresh ice / snow = ', lsub    , ' J/kg'
+         WRITE(numout,*) '      latent heat of subl.  of fresh ice / snow = ', lsub    , ' J/kg'
 #else
-         WRITE(numout,*) '          density times specific heat for snow      = ', rcpsn   , ' J/m^3/K' 
-         WRITE(numout,*) '          density times specific heat for ice       = ', rcpic   , ' J/m^3/K'
-         WRITE(numout,*) '          volumetric latent heat fusion of sea ice  = ', xlic    , ' J/m' 
-         WRITE(numout,*) '          latent heat of sublimation of snow        = ', xsn     , ' J/kg' 
+         WRITE(numout,*) '      density times specific heat for snow      = ', rcpsn   , ' J/m^3/K' 
+         WRITE(numout,*) '      density times specific heat for ice       = ', rcpic   , ' J/m^3/K'
+         WRITE(numout,*) '      volumetric latent heat fusion of sea ice  = ', xlic    , ' J/m' 
+         WRITE(numout,*) '      latent heat of sublimation of snow        = ', xsn     , ' J/kg' 
 #endif
-         WRITE(numout,*) '          volumetric latent heat fusion of snow     = ', xlsn    , ' J/m^3' 
-         WRITE(numout,*) '          density of sea ice                        = ', rhoic   , ' kg/m^3'
-         WRITE(numout,*) '          density of snow                           = ', rhosn   , ' kg/m^3'
-         WRITE(numout,*) '          density of freshwater (in melt ponds)     = ', rhofw   , ' kg/m^3'
-         WRITE(numout,*) '          emissivity of snow or ice                 = ', emic  
-         WRITE(numout,*) '          salinity of ice                           = ', sice    , ' psu'
-         WRITE(numout,*) '          salinity of sea                           = ', soce    , ' psu'
-         WRITE(numout,*) '          latent heat of evaporation (water)        = ', cevap   , ' J/m^3' 
-         WRITE(numout,*) '          correction factor for solar radiation     = ', srgamma 
-         WRITE(numout,*) '          von Karman constant                       = ', vkarmn 
-         WRITE(numout,*) '          Stefan-Boltzmann constant                 = ', stefan  , ' J/s/m^2/K^4'
+         WRITE(numout,*) '      volumetric latent heat fusion of snow     = ', xlsn    , ' J/m^3' 
+         WRITE(numout,*) '      density of sea ice                        = ', rhoic   , ' kg/m^3'
+         WRITE(numout,*) '      density of snow                           = ', rhosn   , ' kg/m^3'
+         WRITE(numout,*) '      density of freshwater (in melt ponds)     = ', rhofw   , ' kg/m^3'
+         WRITE(numout,*) '      emissivity of snow or ice                 = ', emic  
+         WRITE(numout,*) '      salinity of ice                           = ', sice    , ' psu'
+         WRITE(numout,*) '      salinity of sea                           = ', soce    , ' psu'
+         WRITE(numout,*) '      latent heat of evaporation (water)        = ', cevap   , ' J/m^3' 
+         WRITE(numout,*) '      correction factor for solar radiation     = ', srgamma 
+         WRITE(numout,*) '      von Karman constant                       = ', vkarmn 
+         WRITE(numout,*) '      Stefan-Boltzmann constant                 = ', stefan  , ' J/s/m^2/K^4'
          WRITE(numout,*)
-         WRITE(numout,*) '          conversion: degre ==> radian          rad = ', rad
+         WRITE(numout,*) '      conversion: degre ==> radian          rad = ', rad
          WRITE(numout,*)
-         WRITE(numout,*) '          smallest real computer value       rsmall = ', rsmall
+         WRITE(numout,*) '      smallest real computer value       rsmall = ', rsmall
       ENDIF