Changeset 1388 for trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2

Timestamp:

Mar 5, 2015, 11:12:22 AM (10 years ago)

Author:

aslmd

Message:

added radiative heating for NOPHYS. corrected a problem with Q2 dimensions.

File:

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

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

@@ -123 +123 @@
 REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(IN ) :: &
      dz8w,p8w,p,exner,t,t8w,rho,u,v,w,z,th
-REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
+!REAL, DIMENSION( ims:ime, kms:kme+1, jms:jme ), INTENT(INOUT ) :: &
+!     M_Q2
+REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(INOUT ) :: &
      M_Q2
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -164 +166 @@
    REAL, DIMENSION( ims:ime, jms:jme ), INTENT(INOUT)  :: PSFC,TSK
    REAL, DIMENSION( ims:ime, kms:kme, jms:jme ), INTENT(OUT)  :: RTHBLTEN,RUBLTEN,RVBLTEN
+   REAL,DIMENSION(:),ALLOCATABLE,SAVE :: dtsw,dtlw
 #endif
 !-------------------------------------------
@@ -855 +858 @@
 ! specific treatment for tracers !
 !--------------------------------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 IF (MARS_MODE .EQ. 0) THEN
@@ -886 +890 @@
 END SELECT
 #endif
+#endif
 
 !!**********************************************************!!
@@ -1043 +1048 @@
 ! Tracer at surface !
 !-------------------!
+#ifndef NOPHYS
 SELECT CASE (MARS_MODE)
     CASE(0)  
@@ -1080 +1086 @@
 
 END SELECT
+#endif
 
 !-------------------!
@@ -1219 +1226 @@
 ! Tracers !  
 !---------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 pq(subs,:,:) = q_prof(:,:)  !! traceurs generiques, seuls noms sont specifiques
@@ -1234 +1242 @@
     STOP    
 END SELECT
+#endif
 #endif
 ENDDO
@@ -1486 +1495 @@
 ENDIF
 
+#ifdef NOPHYS
+!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!
+IF (firstcall .EQV. .true.) THEN
+  ALLOCATE(dtsw(nlayer))
+  ALLOCATE(dtlw(nlayer))
+  dtsw(1:nlayer) = M_Q2(1,1:nlayer,1)
+  dtlw(1:nlayer) = M_Q2(2,1:nlayer,2)
+ENDIF
+!print *, "SW", dtsw(kps:kpe) ! sw
+!print *, "LW", dtlw(kps:kpe) ! lw
+DO i= 1,ngrid
+  pdt(i,kps:kpe) = dtsw(1:nlayer) + dtlw(1:nlayer) 
+ENDDO
+!!!!!!!!!!!!!!!!!!!
+!!!!!!!!!!!!!!!!!!!
+#endif
+
 DO j = jps,jpe
 DO i = ips,ipe
@@ -1548 +1575 @@
 ! Tracers !  
 !---------!
+#ifndef NOPHYS
 #ifdef NEWPHYS
 SCALAR(i,kps:kpe,j,1)=0.
@@ -1591 +1619 @@
 END SELECT
 #endif
+#endif
         !!TODO: check if adding the whole tendency once, and set the
         !!TODO: following tendencies to 0 until physics is called again