Changeset 1709 for trunk/LMDZ.GENERIC/libf

Timestamp:

May 16, 2017, 9:57:31 AM (8 years ago)

Author:

aslmd

Message:

Added surfalbedo and surfemis keywords to be used with startphy_file = .false. Also made default values in phyetat0_mod unambiguously float

Location:

trunk/LMDZ.GENERIC/libf/phystd

Files:

• callcorrk.F90 (modified) (2 diffs)
• callkeys_mod.F90 (modified) (1 diff)
• phyetat0_mod.F90 (modified) (19 diffs)

trunk/LMDZ.GENERIC/libf/phystd/callcorrk.F90

@@ -20 +20 @@
       USE tracer_h
       use comcstfi_mod, only: pi, mugaz, cpp
-      use callkeys_mod, only: varactive,diurnal,tracer,water,nosurf,varfixed,satval,        &
+      use callkeys_mod, only: varactive,diurnal,tracer,water,varfixed,satval,        &
                               kastprof,strictboundcorrk,specOLR,CLFvarying
 
@@ -494 +494 @@
             albv(nw)=albedo(ig,nw)
          ENDDO
-
-         if (nosurf) then ! Case with no surface.
-            DO nw=1,L_NSPECTV
-               if(albv(nw).gt.0.0) then
-                  print*,'For open lower boundary in callcorrk must'
-                  print*,'have spectral surface band albedos all set to zero!'
-                  call abort
-               endif
-            ENDDO         
-         endif
 
       if ((ngrid.eq.1).and.(global1d)) then ! Fixed zenith angle 'szangle' in 1D simulations w/ globally-averaged sunlight.

trunk/LMDZ.GENERIC/libf/phystd/callkeys_mod.F90

@@ -107 +107 @@
       real,save :: MassPlanet
 !$OMP THREADPRIVATE(flatten,Rmean,J2,MassPlanet)
+      real,save :: surfalbedo
+      real,save :: surfemis
+!$OMP THREADPRIVATE(surfalbedo,surfemis)
 
       logical,save :: iscallphys=.false.!existence of callphys.def

trunk/LMDZ.GENERIC/libf/phystd/phyetat0_mod.F90

@@ -10 +10 @@
                      emis,q2,qsurf,cloudfrac,totcloudfrac,hice, &
                      rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
+
+! to use  'getin_p'
+      use ioipsl_getin_p_mod, only: getin_p
 
   use tabfi_mod, only: tabfi
@@ -18 +21 @@
                      inquire_dimension, inquire_dimension_length
   use slab_ice_h, only: noceanmx
-  use callkeys_mod, only: CLFvarying
+  use callkeys_mod, only: CLFvarying,surfalbedo,surfemis
 
   implicit none
@@ -115 +118 @@
   endif
 else
-  phisfi(:)=0
+  phisfi(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: surface geopotential <phisfi> range:", &
@@ -127 +130 @@
   endif
 else
-  albedodat(:)=0.5 ! would be better to read value from def file...
+  ! If no startfi file, use parameter surfalbedo in def file
+  surfalbedo=0.5
+  call getin_p("surfalbedo",surfalbedo)
+  print*,"surfalbedo",surfalbedo
+  albedodat(:)=surfalbedo
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Bare ground albedo <albedodat> range:", &
@@ -139 +146 @@
   endif
 else
-  zmea(:)=0
+  zmea(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: <ZMEA> range:", &
@@ -151 +158 @@
   endif
 else
-  zstd(:)=0
+  zstd(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: <ZSTD> range:", &
@@ -163 +170 @@
   endif
 else
-  zsig(:)=0
+  zsig(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: <ZSIG> range:", &
@@ -175 +182 @@
   endif
 else
-  zgam(:)=0
+  zgam(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: <ZGAM> range:", &
@@ -187 +194 @@
   endif
 else
-  zthe(:)=0
+  zthe(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: <ZTHE> range:", &
@@ -199 +206 @@
   endif
 else
-  tsurf(:)=0 ! will be updated afterwards in physiq !
+  tsurf(:)=0. ! will be updated afterwards in physiq !
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Surface temperature <tsurf> range:", &
@@ -211 +218 @@
   endif
 else
-  emis(:)=1 ! would be better to read value from def file...
+  ! If no startfi file, use parameter surfemis in def file
+  surfemis=1.0
+  call getin_p("surfemis",surfemis)
+  print*,"surfemis",surfemis
+  emis(:)=surfemis
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Surface emissivity <emis> range:", &
@@ -224 +235 @@
     endif
   else
-    cloudfrac(:,:)=0
+    cloudfrac(:,:)=0.0
   endif ! of if (startphy_file)
   write(*,*) "phyetat0: Cloud fraction <cloudfrac> range:", &
@@ -240 +251 @@
     endif
   else
-    totcloudfrac(:)=0
+    totcloudfrac(:)=0.0
   endif ! of if (startphy_file)
   write(*,*) "phyetat0: Total cloud fraction <totcloudfrac> range:", &
@@ -254 +265 @@
     write(*,*) "phyetat0: Failed loading <hice>"
     !  call abort
-    hice(:)=0
-  endif
-else
-  hice(:)=0
+    hice(:)=0.
+  endif
+else
+  hice(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Height of oceanic ice <hice> range:", &
@@ -279 +290 @@
   endif ! of if (.not.found)
 else
-  rnat(:)=1
+  rnat(:)=1.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Nature of surface <rnat> range:", &
@@ -292 +303 @@
   endif
 else
-  pctsrf_sic(:)=0
+  pctsrf_sic(:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Pourcentage of sea ice cover <pctsrf_sic> range:", &
@@ -335 +346 @@
   endif
 else
-  tsea_ice(1:ngrid)=0
+  tsea_ice(1:ngrid)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: Oceanic ice quantity <sea_ice> range:", &
@@ -348 +359 @@
   endif
 else
-  q2(:,:)=0
+  q2(:,:)=0.
 endif ! of if (startphy_file)
 write(*,*) "phyetat0: PBL wind variance <q2> range:", &
@@ -365 +376 @@
       endif
     else
-      qsurf(:,iq)=0
+      qsurf(:,iq)=0.
     endif ! of if (startphy_file)
     write(*,*) "phyetat0: Surface tracer <",trim(txt),"> range:", &