Changeset 1974 for trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

Timestamp:

Jul 18, 2018, 4:48:34 PM (6 years ago)

Author:

mvals

Message:

Mars GCM:
Integration of the detached dust layer parametrizations (rocket dust storm, slope wind lifting, CW, and dust injection scheme, DB).
Still experimental, default behaviour (rdstorm=.false., dustinjection=0) identical to previous revision.
NB: Updated newstart requires an updated "surface.nc" containing the "hmons" field.
EM+MV

File:

• trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/libf/phymars/phyetat0_mod.F90

@@ -12 +12 @@
   use tracer_mod, only: noms ! tracer names
   use surfdat_h, only: phisfi, albedodat, z0, z0_default,&
-                       zmea, zstd, zsig, zgam, zthe
+                       zmea, zstd, zsig, zgam, zthe, hmons
   use iostart, only: nid_start, open_startphy, close_startphy, &
                      get_field, get_var, inquire_field, &
@@ -19 +19 @@
 
   implicit none
+  
+  include "callkeys.h"
 !======================================================================
 ! Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
@@ -179 +181 @@
 endif
 
+! hmons 
+call get_field("hmons",hmons,found)
+if (.not.found) then
+  write(*,*) "WARNING: phyetat0: Failed loading <hmons>"
+  if (rdstorm) then
+    call abort
+  else
+    write(*,*) "will continue anyway..."
+    write(*,*) "because you may not need it."
+    hmons(:)=0.
+  end if
+else
+  do ig=1,ngrid
+    if (hmons(ig) .eq. -999999.)  hmons(ig)=0.
+  enddo
+  write(*,*) "phyetat0: <hmons> range:", &
+             minval(hmons), maxval(hmons)
+endif
       
 ! Time axis
@@ -277 +297 @@
 if (.not.found) then
   write(*,*) "phyetat0: <tauscaling> not in file"
-  tauscaling(:) = -1
+  tauscaling(:) = 1
 else
   write(*,*) "phyetat0: dust conversion factor <tauscaling> range:", &
@@ -296 +316 @@
 call get_field("wstar",wstar,found,indextime)
 if (.not.found) then
-  write(*,*) "phyetat0: <totcloudfrac> not in file! Set to zero"
+  write(*,*) "phyetat0: <wstar> not in file! Set to zero"
   wstar(:)=0
 else