Changeset 91 for trunk/mars

Timestamp:

Mar 8, 2011, 5:59:28 PM (14 years ago)

Author:

aslmd

Message:

mars: test outliers [dans initracer.F, commente] LMD_MM_MARS: modifications mineures [retrocompatibilite ancienne physique, callphys.def test pour nouvelle physique] PLOT: generalisation de la routine map_uvt pour pouvoir tracer des figures en projection polaire complete

File:

• trunk/mars/libf/phymars/initracer.F (modified) (2 diffs)

trunk/mars/libf/phymars/initracer.F

@@ -46 +46 @@
 #include "chimiedata.h"
 
+#include "surfdat.h"
 
       real qsurf(ngridmx,nqmx)       ! tracer on surface (e.g.  kg.m-2)
       real co2ice(ngridmx)           ! co2 ice mass on surface (e.g.  kg.m-2)
-      integer iq,ig,count
+      integer iq,ig,count, yeyey
       real r0_lift , reff_lift, nueff_lift
 c     Ratio of small over large dust particles (used when both 
@@ -501 +502 @@
 c Perennial H20 north cap defined by watercaptag=true (allows surface to be
 c hollowed by sublimation in vdifc).
+         yeyey = 0
          do ig=1,ngridmx
+!          !!! TESTS TESTS outliers
+!          !!! TESTS TESTS outliers
+!          if ( ( lati(ig)*180./pi      .ge.  75 ) .and.
+!     .         ( lati(ig)*180./pi      .le.  77 ) .and.
+!     .         ( ( ( long(ig)*180./pi .ge. 000. ) .and.
+!     .              ( long(ig)*180./pi .le. 120. ) ) 
+!     .             .or.
+!     .             ( ( long(ig)*180./pi .ge. -130. ) .and.
+!     .             ( long(ig)*180./pi .le. -115. ) ) ) ) then
+!             if (yeyey .eq. 0) then  !!! 1/2 en 64x48 sinon trop large en lat
+!              write(*,*) "outliers ", lati(ig)*180./pi, long(ig)*180./pi
+!              if (ngridmx.ne.1) watercaptag(ig)=.true.
+!              dryness(ig) = 1.
+!              albedodat(ig) = 0.45 !! comme alb_surfice
+!              yeyey = 1
+!             else
+!              yeyey = 0
+!             endif
+!          endif
+!          !!! TESTS TESTS outliers
+!          !!! TESTS TESTS outliers
+!
+!          !!! TESTS TESTS addcap
+!          !!! TESTS TESTS addcap     
+!          if ( ( lati(ig)*180./pi      .ge.  82 ) .and.
+!     .         ( lati(ig)*180./pi      .le.  84 ) .and.
+!     .         ( ( long(ig)*180./pi .gt. -030. ) .and.
+!     .              ( long(ig)*180./pi .lt. 090. ) ) ) then
+!              write(*,*) "capadd ", lati(ig)*180./pi, long(ig)*180./pi
+!              if (ngridmx.ne.1) watercaptag(ig)=.true.
+!              albedodat(ig) = 0.45 !! comme alb_surfice
+!              dryness(ig) = 1.
+!          endif
+!          !!! TESTS TESTS addcap
+!          !!! TESTS TESTS addcap 
+    
            if (lati(ig)*180./pi.gt.84) then
              if (ngridmx.ne.1) watercaptag(ig)=.true.