Changeset 1088 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Oct 31, 2013, 2:49:26 PM (12 years ago)

Author:

tnavarro

Message:

Added a freedust mode, to be used with doubleq (for data assimilation). Dust is not lifted, not rescaled, dust opacity is predicted instead of being forced.

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• aeropacity.F (modified) (3 diffs)
• callkeys.h (modified) (2 diffs)
• inifis.F (modified) (3 diffs)

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

@@ -205 +205 @@
 c     Vertical column optical depth at "odpref" Pa 
 c     ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
-      IF(iaervar.eq.1) THEN 
+      IF(freedust) THEN
+         tauref(:) = 0. ! tauref is computed after, instead of being forced
+
+      ELSE IF(iaervar.eq.1) THEN 
          do ig=1, ngrid
           tauref(ig)=max(tauvis,1.e-9) ! tauvis=cste (set in callphys.def
@@ -431 +434 @@
 c -----------------------------------------------------------------
 
-c     Temporary scaling factor
-      taudusttmp(1:ngrid)=0.
-      DO iaer=1,naerdust
-        DO l=1,nlayer
-          DO ig=1,ngrid
-c           Scaling factor
-            taudusttmp(ig) = taudusttmp(ig) + 
-     &                       aerosol(ig,l,iaerdust(iaer))
-          ENDDO
-        ENDDO
-      ENDDO
-
-c     Saved scaling factor
-      DO ig=1,ngrid
-          tauscaling(ig) = tauref(ig) * 
-     &                     pplev(ig,1) / odpref / taudusttmp(ig)
-c          tauscaling(ig) = 1.e-4
-      ENDDO
+      IF (freedust) THEN
+          tauscaling(:) = 1.
+
+      ELSE
+c       Temporary scaling factor
+        taudusttmp(1:ngrid)=0.
+        DO iaer=1,naerdust
+          DO l=1,nlayer
+            DO ig=1,ngrid
+c             Scaling factor
+              taudusttmp(ig) = taudusttmp(ig) + 
+     &                         aerosol(ig,l,iaerdust(iaer))
+            ENDDO
+          ENDDO
+        ENDDO
+
+c       Saved scaling factor
+        DO ig=1,ngrid
+            tauscaling(ig) = tauref(ig) * 
+     &                       pplev(ig,1) / odpref / taudusttmp(ig)
+c            tauscaling(ig) = 1.e-4
+        ENDDO
+
+      ENDIF
 
 c     Opacity computation
@@ -459 +468 @@
         ENDDO
       ENDDO
+
+      IF (freedust) THEN
+        ! tauref has been initialized to 0 before.
+        DO l=1,nlayer
+          DO ig=1,ngrid
+            tauref(ig) = tauref(ig) + 
+     &                aerosol(ig,l,iaerdust(iaer))
+          ENDDO
+        ENDDO
+      ENDIF
       
 c output for debug

trunk/LMDZ.MARS/libf/phymars/callkeys.h

@@ -11 +11 @@
      &   ,callmolvis,callmoldiff,thermochem,thermoswater,callemis       &
      &   ,callg2d,linear,rayleigh,tracer,active,doubleq,submicron       &
-     &   ,lifting,callddevil,scavenging,sedimentation,activice,water    &
-     &   ,tifeedback,microphys,caps,photochem,calltherm                 &
+     &   ,lifting,freedust,callddevil,scavenging,sedimentation          &
+     &   ,activice,water,tifeedback,microphys,caps,photochem,calltherm  &
      &   ,callrichsl,callslope,tituscap
      
@@ -53 +53 @@
       logical tracer
       integer dustbin
+      logical freedust 
       logical active,doubleq,submicron,lifting,callddevil,scavenging
       logical sedimentation

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

@@ -401 +401 @@
            stop
          endif
+
 ! lifting
          write(*,*)"dust lifted by GCM surface winds ?"
@@ -414 +415 @@
            stop
          endif
+
+! free evolving dust
+! freedust=true just says that there is no lifting and no dust opacity scaling.
+         write(*,*)"dust lifted by GCM surface winds ?"
+         freedust=.false. ! default value
+         call getin("freedust",freedust)
+         write(*,*)" freedust = ",freedust
+         if (freedust.and..not.doubleq) then
+           print*,'freedust should be used with doubleq !'
+           stop
+         endif
+         if (freedust.and.lifting) then
+           print*,'if freedust is used, then lifting should not be used'
+           print*,'lifting forced to false !!'
+           lifting=.false.
+         endif
+
 ! callddevil
          write(*,*)" dust lifted by dust devils ?"
@@ -493 +511 @@
          call getin("microphys",microphys)
          write(*,*)" microphys = ",microphys
+         if (freedust.and.microphys) then
+           print*,'freedust can not be used with microphys !'
+           print*,'(although that could be improved...)'
+           stop
+         endif
 
 ! microphysical parameter contact