Changeset 480

Timestamp:

Dec 19, 2011, 2:48:57 PM (13 years ago)

Author:

tnavarro

Message:

new option in newstart: ini_h2osurf, water ice surface initialisation to keep seasonal frost

but get rid of polar cap positive or negative values.

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/dyn3d/newstart.F (modified) (2 diffs)

trunk/LMDZ.MARS/README

@@ -1352 +1352 @@
 >> minor correction in aeronomars/calchim.F and aeronomars/surfacearea.F : 
    calls to wstats should always be embeded in if (callstats) condition.
+   
+== 19/12/2011 == TN
+>> new option in newstart: 'ini_h2osurf', water ice surface initialisation to keep seasonal frost 
+   but get rid of polar cap positive or negative values.

trunk/LMDZ.MARS/libf/dyn3d/newstart.F

@@ -487 +487 @@
      $ice '
       write(*,*) 'ini_q-iceH2O : tracers initialisation for chemistry on
-     $ly '
+     $ly ' 
+      write(*,*) 'ini_h2osurf : reinitialize surface water ice '
       write(*,*) 'noglacier : Remove tropical H2O ice if |lat|<45'
       write(*,*) 'watercapn : H20 ice on permanent N polar cap '
@@ -900 +901 @@
          write(*,*) '---------------------------south pole='
      *               ,q(1,jjp1,1,igcm_h2o_vap)
+
+c      ini_h2osurf : reinitialize surface water ice
+c      --------------------------------------------------
+        else if (trim(modif) .eq. 'ini_h2osurf') then
+          write(*,*)'max surface ice left?(e.g. 0.2 kg/m2=200microns)'
+ 207      read(*,*,iostat=ierr) val
+          if(ierr.ne.0) goto 207
+          write(*,*)'also set negative values of surf ice to 0'
+           do ig=1,ngridmx
+              qsurf(ig,igcm_h2o_ice)=min(val,qsurf(ig,igcm_h2o_ice))
+              qsurf(ig,igcm_h2o_ice)=max(0.,qsurf(ig,igcm_h2o_ice))
+           end do
 
 c      noglacier : remove tropical water ice (to initialize high res sim)