Changeset 170 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

Timestamp:

Jun 21, 2011, 11:16:18 AM (13 years ago)

Author:

acolaitis

Message:

17/06/2011 == AC

Added a new package for the LES

• This new mode allows the emission of decaying tracers

both at the surface and at a time-varying location (here, the PBL top).

• Package number is 21, and requires the presence of an "input_zipbl"

file alongside the other input files.

• input_zipbl is expected to contains two columns : the altitude index corresponding

to the height of tracer emission, and the corresponding local time.

• input_zipbl should be generated from re-analysis of previous LES simulations. An example

intput_zipbl has been added in "MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/res/LES/input/"

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F (modified) (8 diffs)

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/module_lmd_driver.F

@@ -162 +162 @@
 ! LOCAL  VARIABLES
 !-------------------------------------------
-   INTEGER ::    i,j,k,its,ite,jts,jte,ij
+   INTEGER ::    i,j,k,its,ite,jts,jte,ij,kk
    INTEGER ::    subs,iii
-   REAL ::    tau_decay
+
+   ! *** for Mars Mode 20 and 21 ***
+   REAL ::    tau_decay, lct, zilctmin
+   INTEGER, SAVE :: zipbl(500) !index of zi in file input_zipbl
+   REAL, SAVE :: zilct(500) !corresponding local time in input_zipbl
+   INTEGER :: lctindex,ziindex
+   LOGICAL :: end_of_file
+   ! *** ----------------------- ***
 
    ! *** for LMD physics
@@ -629 +636 @@
 
 SELECT CASE (MARS_MODE) !! ONLY ALLOW FOR MODES DEFINED IN Registry.EM
-   CASE(4-10,12-19,21:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
+   CASE(4-10,12-19,22:)      !! -- CHANGE THIS if YOU ADDED CASES in REGISTRY.EM
    PRINT *, 'NOT SUPPORTED, to be done'
    STOP
@@ -640 +647 @@
 !package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
 !package   radioac      mars==20                     -              scalar:qtrac1
+!package   radioac2     mars==21                     -              scalar:upward,downward
 !!!!!!!!!!!!!!!!!!! FOR REFERENCE
 
@@ -664 +672 @@
     CASE(20) 
       wtnom(1) = 'qtrac1'
+    CASE(21)
+      wtnom(1) = 'upward'
+      wtnom(2) = 'downward'
 END SELECT
 #endif
-
 
 !!*******************************************!!
@@ -705 +715 @@
 q_prof(:,1:nq) = SCALAR(i,kps:kpe,j,2:nq+1)  !! the names were set above !! one dummy tracer in WRF
   !!! CAS DU CO2
-  DO iii=1,nq 
-   IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
-  ENDDO
-
-    !! Mars mode 20 : add a passive tracer with radioactive-like decay
-    !! INIT HERE
-    IF (firstcall .EQV. .true.) THEN
-    IF (MARS_MODE .EQ. 20) THEN
+DO iii=1,nq 
+ IF ( wtnom(iii) .eq. 'co2' ) q_prof(:,iii) = 0.95
+ENDDO
+
+IF ((MARS_MODE .EQ. 20) .OR. (MARS_MODE .EQ. 21)) THEN
+   IF (firstcall .EQV. .true.) THEN
       q_prof(:,:) = 0.95
-    ENDIF
-    ENDIF
+   ENDIF
+ENDIF
 
 #else
@@ -863 +871 @@
     CASE(20)
     qsurf_val(:)=0.
+    CASE(21)
+    qsurf_val(:)=0.
 #else
     CASE(3:)
@@ -1197 +1207 @@
 ! --is the one calculated during the last call to physics          !
 !------------------------------------------------------------------!
+
+! PBL top index reading for MARS_MODE 21 :
+IF (MARS_MODE .EQ. 21) THEN
+      IF (firstcall .EQV. .true.) THEN
+           open(unit=15,file='input_zipbl',form='formatted',status='old')
+           rewind(15)
+           end_of_file = .false.
+           kk = 0
+           do while (.not. end_of_file)
+             read(15,*,end=101) zipbl(kk+1),zilct(kk+1)
+             write(*,*) kk, zipbl(kk+1),zilct(kk+1)
+             kk = kk+1
+             go to 112
+101          end_of_file = .true.
+112          continue
+           enddo
+           close(unit=15,status = 'keep')
+      ENDIF
+      lct=lct_input + elaps/3700.
+      zilctmin=MINVAL(ABS(zilct(:)-lct))
+      lctindex=0
+      DO kk=1,500
+        IF(ABS(zilct(kk)-lct) .eq. zilctmin) lctindex=kk
+      ENDDO
+      IF(lctindex .eq. 0) print*, 'probleme index'
+        ziindex=zipbl(lctindex)
+ENDIF
+
 DO j = jps,jpe
 DO i = ips,ipe
@@ -1228 +1266 @@
 #ifdef NEWPHYS
 SCALAR(i,kps:kpe,j,1)=0.
-!! Mars mode 20 : add a passive tracer with radioactive-like decay
-IF (MARS_MODE .EQ. 20) THEN
+
+SELECT CASE (MARS_MODE)
+CASE(20)
+   !! Mars mode 20 : add a passive tracer with radioactive-like decay
    IF ( (i == ips) .AND. (j == jps) )   print *, 'RADIOACTIVE-LIKE TRACER WITH SOURCE AT SURFACE LAYER.'
    tau_decay=60.*10. !! why not make it a namelist argument?
    SCALAR(i,kps:kpe,j,2) = SCALAR(i,kps:kpe,j,2)*exp(-dt/tau_decay)
    SCALAR(i,1,j,2) = SCALAR(i,1,j,2) + 1. !! this tracer is emitted in the surface layer
-ELSE
+CASE(21)
+   !! Mars mode 21 : add a passive tracer with radioactive-like decay at surface and pbl top
+   IF ( (i == ips) .AND. (j == jps) )   print *, 'RADIOACTIVE-LIKE TRACER WITH SOURCE AT SURFACE LAYER AND PBL TOP.'
+   tau_decay=60.*10. !! why not make it a namelist argument?
+   SCALAR(i,kps:kpe,j,2) = SCALAR(i,kps:kpe,j,2)*exp(-dt/tau_decay)
+   SCALAR(i,kps:kpe,j,3) = SCALAR(i,kps:kpe,j,3)*exp(-dt/tau_decay)
+   SCALAR(i,1,j,2) = SCALAR(i,1,j,2) + 1. !! this tracer is emitted in the surface layer
+   IF ( (i == ips) .AND. (j == jps) )   print *, 'index of pbl emission: ',ziindex
+   IF (ziindex .NE. 0) THEN
+   SCALAR(i,ziindex,j,3) = SCALAR(i,ziindex,j,3) + 1. !! this tracer is emitted at pbl top
+   ENDIF
+CASE DEFAULT
    SCALAR(i,kps:kpe,j,2:nq+1)=SCALAR(i,kps:kpe,j,2:nq+1)+pdq(subs,kps:kpe,1:nq)*dt  !!! here dt is needed
-ENDIF
+END SELECT
 #else
 SELECT CASE (MARS_MODE)