Changeset 1041 for trunk/MESOSCALE

Timestamp:

Sep 19, 2013, 1:16:54 AM (11 years ago)

Author:

aslmd

Message:

MESOSCALE LMDZ.MARS. added outliers in old physics. get rid of negative qsurf in MESOINI mode.

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/surfini.F (modified) (2 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/mars_lmd/libf/phymars/surfini.F

@@ -14 +14 @@
 #include "callkeys.h"
 #include "tracer.h"
+#include "comgeomfi.h"
+#include "comcstfi.h"
+#include "watercap.h"
+
 c
-      INTEGER ngrid,ig,icap
+      INTEGER ngrid,ig,icap,iq
       REAL  piceco2(ngrid),psolaralb(ngrid,2)
       REAL qsurf(ngrid,nqmx) !tracer on surface (kg/m2)
@@ -23 +27 @@
 c
 c=======================================================================
+
+      do ig=1,ngridmx
+
+         !write(*,*) "all qsurf to zero. dirty."
+         do iq=1,nqmx
+         qsurf(ig,iq)=0.  !! on jette les inputs GCM
+                          !! on regle juste watercaptag
+                          !! il faudrait garder les inputs GCM 
+                          !! si elles sont consequentes
+         enddo
+         if ( ( lati(ig)*180./pi .gt. 70. ) .and.
+     .        ( albedodat(ig) .ge. 0.26   ) )  then
+                 write(*,*)"outlier ",ig,albedodat(ig),inertiedat(ig)
+                 watercaptag(ig)  = .true.
+                 dryness(ig)      = 1.
+                 albedodat(ig)    = 0.45 !albedo_h2o_ice  !! pour output 
+                 inertiedat(ig)   = 800.
+                 write(*,*)"new values ",ig,albedodat(ig),inertiedat(ig)
+         else
+                 watercaptag(ig)  = .false.
+                 dryness(ig)      = 1.
+         endif 
+         if ( inertiedat(ig) .ge. 800. ) then
+            write(*,*)"change inertie from ",inertiedat(ig)," to 800."
+            inertiedat(ig)   = 800.
+         endif
+         
+      enddo
+
 
 c