Changeset 835 for trunk/LMDZ.MARS/libf

Timestamp:

Nov 8, 2012, 7:42:44 PM (12 years ago)

Author:

jbmadeleine

Message:

Mars GCM:
Cleaned physiq.F (some debugging lines were introduced at revision
833, sorry about that!);

File:

• trunk/LMDZ.MARS/libf/phymars/physiq.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq.F

@@ -621 +621 @@
 c          ------------------
 
-c-------------------------------------------
-c     MODIF_JBM made on 12W29D5,05:29:23 PM
-c-------------------------------------------
-c     WRITE(*,*) "MODIF_JBM DEBUGGING 12W29D5,05:29:23 PM" 
-c     WRITE(*,*) "BEFORE CALLRADITE"
-c     WRITE(*,*) ">>> rice = ",rice(11,1)
-c-------------------------------------------
-
            CALL callradite(icount,ngrid,nlayer,nq,zday,zls,pq,albedo,
      $     emis,mu0,pplev,pplay,pt,tsurf,fract,dist_sol,igout,
@@ -634 +626 @@
      $     tauref,tau,aerosol,tauscaling,taucloudtes,rdust,rice,
      $     nuice,co2ice)
-
-c-------------------------------------------
-c     MODIF_JBM made on 12W29D5,05:29:23 PM
-c-------------------------------------------
-c     WRITE(*,*) "MODIF_JBM DEBUGGING 12W29D5,05:29:23 PM"
-c     WRITE(*,*) "AFTER CALLRADITE"
-c     WRITE(*,*) ">>> rice = ",rice(11,1)
-c-------------------------------------------
-
 
 c          Outputs for basic check (middle of domain)
@@ -1038 +1021 @@
      &                rsedcloud,rhocloud)
 
-c-------------------------------------------
-c     MODIF_JBM made on 12W29D5,05:29:23 PM
-c-------------------------------------------
-c     WRITE(*,*) "MODIF_JBM DEBUGGING 12W29D5,05:29:23 PM"
-c     WRITE(*,*) "AFTER WATERCLOUD"
-c     WRITE(*,*) ">>> rice = ",rice(11,1)
-c-------------------------------------------
-     
 c Temperature variation due to latent heat release
            if (activice) then
@@ -1157 +1132 @@
      &                tau,tauscaling)
 
-c-------------------------------------------
-c     MODIF_JBM made on 12W29D5,05:29:23 PM
-c-------------------------------------------
-c     WRITE(*,*) "MODIF_JBM DEBUGGING 12W29D5,05:29:23 PM"
-c     WRITE(*,*) "AFTER CALLSEDIM"
-c     WRITE(*,*) ">>> rice = ",rice(11,1)
-c------------------------------------------- 
-     
            DO iq=1, nq
              DO l=1,nlayer
@@ -1419 +1386 @@
       ENDDO
 c ******* TEST ******************************************************
-c MODIF_JBM COMMENTED
-c      ztim1 = 999
-c      DO l=1,nlayer
-c        DO ig=1,ngrid
-c           if (pt(ig,l).lt.ztim1) then
-c               ztim1 = pt(ig,l)
-c               igmin = ig
-c               lmin = l 
-c           end if
-c        ENDDO
-c      ENDDO
-c      if(min(pt(igmin,lmin),zt(igmin,lmin)).lt.70.) then
-c        write(*,*) 'PHYSIQ: stability WARNING :'
-c        write(*,*) 'pt, zt Tmin = ', pt(igmin,lmin), zt(igmin,lmin),
-c     &              'ig l =', igmin, lmin
-c      end if
+      ztim1 = 999
+      DO l=1,nlayer
+        DO ig=1,ngrid
+           if (pt(ig,l).lt.ztim1) then
+               ztim1 = pt(ig,l)
+               igmin = ig
+               lmin = l 
+           end if
+        ENDDO
+      ENDDO
+      if(min(pt(igmin,lmin),zt(igmin,lmin)).lt.70.) then
+        write(*,*) 'PHYSIQ: stability WARNING :'
+        write(*,*) 'pt, zt Tmin = ', pt(igmin,lmin), zt(igmin,lmin),
+     &              'ig l =', igmin, lmin
+      end if
 c *******************************************************************
 
@@ -1901 +1867 @@
 c        WRITEDIAGFI can ALSO be called from any other subroutines
 c        for any variables !!
-         call WRITEDIAGFI(ngrid,"emis","Surface emissivity","w.m-1",2,
-     &                  emis)
-         call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",3,zplay)
-         call WRITEDIAGFI(ngrid,"pplev","Pressure","Pa",3,zplev)
+c        call WRITEDIAGFI(ngrid,"emis","Surface emissivity","w.m-1",2,
+c    &                  emis)
+c        call WRITEDIAGFI(ngrid,"pplay","Pressure","Pa",3,zplay)
+c        call WRITEDIAGFI(ngrid,"pplev","Pressure","Pa",3,zplev)
          call WRITEDIAGFI(ngrid,"tsurf","Surface temperature","K",2,
      &                  tsurf)
@@ -1931 +1897 @@
 c        call WRITEDIAGFI(ngrid,"ssurf","Surface stress","N.m-2",2,
 c    &                  zstress)
-         call WRITEDIAGFI(ngridmx,'sw_htrt','sw heat. rate',
-     &                   'w.m-2',3,zdtsw)
-         call WRITEDIAGFI(ngridmx,'lw_htrt','lw heat. rate',
-     &                   'w.m-2',3,zdtlw)
+c        call WRITEDIAGFI(ngridmx,'sw_htrt','sw heat. rate',
+c    &                   'w.m-2',3,zdtsw)
+c        call WRITEDIAGFI(ngridmx,'lw_htrt','lw heat. rate',
+c    &                   'w.m-2',3,zdtlw)
             if (.not.activice) then
                CALL WRITEDIAGFI(ngridmx,'tauTESap',