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

Timestamp:

Oct 31, 2013, 2:49:26 PM (11 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.

File:

• trunk/LMDZ.MARS/libf/phymars/inifis.F (modified) (3 diffs)

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