Changeset 2088 for trunk/MESOSCALE/LMD_MM_MARS/SRC

Timestamp:

Feb 4, 2019, 3:22:38 PM (6 years ago)

Author:

mlefevre

Message:

MESOSCALE. Argurment of lmd_driver was called twice and a interpolation function was missing

Location:

trunk/MESOSCALE/LMD_MM_MARS/SRC

Files:

• LES/modif_mars/module_first_rk_step_part1.F (modified) (1 diff)
• WRFV2/dyn_em/module_init_utilities.F (modified) (1 diff)

trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_first_rk_step_part1.F

@@ -672 +672 @@
      &       ,LATENT_HF=grid%LATENT_HF  &
      &       ,SWDOWNZ=grid%SWDOWNZ   &
-     &       ,SWDOWNZ=grid%SWDOWNZ   &
      &       ,TAU_DUST=grid%TAU_DUST   &
      &       ,RDUST=grid%RDUST   &

trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_init_utilities.F

@@ -179 +179 @@
  END FUNCTION interp_0_log
 
+ real function interp_0_log2( v_in,  &
+                         p_in, p_out, nz_in  )
+ implicit none
+ integer nz_in, nz_out
+ real*8  v_in(nz_in)
+ real*8  p_in(nz_in)
+ real*8  p_out
+
+ integer kp, k, im, ip
+ logical interp, increasing_z 
+ real    height, w1, w2
+ logical debug
+ parameter ( debug = .false. )
+
+! does vertical coordinate increase or decrease with increasing k?
+! set offset appropriately
+
+ height = p_out
+
+ if(debug) write(6,*) ' height in interp_0 ',height
+
+ if (p_in(nz_in) .gt. p_in(1)) then
+
+    if(debug) write(6,*) ' monotonic increase in z in interp_0 '
+    IF (height > p_in(nz_in)) then
+      if(debug) write(6,*) ' point 1 in interp_0 '
+      w2 = log(p_in(nz_in)/height)/log(p_in(nz_in)/p_in(nz_in-1))
+      w1 = 1.-w2
+      interp_0_log2 = w1*v_in(nz_in) + w2*v_in(nz_in-1)
+    ELSE IF (height < p_in(1)) then
+      if(debug) write(6,*) ' point 2 in interp_0 '
+      w2 = log(p_in(2)/height)/log(p_in(2)/p_in(1))
+      w1 = 1.-w2
+      interp_0_log2 = w1*v_in(2) + w2*v_in(1)
+    ELSE
+      if(debug) write(6,*) ' point 3 in interp_0 '
+      interp = .false.
+      kp = nz_in
+      DO WHILE ( (interp .eqv. .false.) .and. (kp .ge. 2) )
+        IF(   ((p_in(kp)   .ge. height) .and.     &
+               (p_in(kp-1) .le. height))        )   THEN
+          w2 = log(height/p_in(kp))/log(p_in(kp-1)/p_in(kp))
+          w1 = 1.-w2
+          interp_0_log2 = w1*v_in(kp) + w2*v_in(kp-1)
+          if(debug) write(6,*) ' interp data, kp, w1, w2 ',kp, w1, w2
+          if(debug) write(6,*) ' interp data, v_in(kp), v_in(kp-1), interp_0_p', &
+                     v_in(kp), v_in(kp-1), interp_0_log2
+          interp = .true.
+        END IF
+        kp = kp-1
+      ENDDO
+    ENDIF
+
+ else
+
+    if(debug) write(6,*) ' monotonic decrease in z in interp_0 '
+
+    IF (height < p_in(nz_in)) then
+      if(debug) write(6,*) ' point 1 in interp_0 '
+      w2 = log(p_in(nz_in)/height)/log(p_in(nz_in)/p_in(nz_in-1))
+      w1 = 1.-w2
+      interp_0_log2 = w1*v_in(nz_in) + w2*v_in(nz_in-1)
+    ELSE IF (height > p_in(1)) then
+      if(debug) write(6,*) ' point 2 in interp_0 '
+      w2 = log(p_in(2)/height)/log(p_in(2)/p_in(1))
+      w1 = 1.-w2
+      interp_0_log2 = w1*v_in(2) + w2*v_in(1)
+    ELSE
+      if(debug) write(6,*) ' point 3 in interp_0 '
+      interp = .false.
+      kp = nz_in
+      height = p_out
+      DO WHILE ( (interp .eqv. .false.) .and. (kp .ge. 2) )
+        IF(   ((p_in(kp)   .le. height) .and.     &
+               (p_in(kp-1) .ge. height))             )   THEN
+          w2 = log(height/p_in(kp))/log(p_in(kp-1)/p_in(kp))
+          w1 = 1.-w2
+          interp_0_log2 = w1*v_in(kp) + w2*v_in(kp-1)
+          interp = .true.
+        END IF
+        kp = kp-1
+      ENDDO
+    ENDIF
+
+ end if
+
+ return
+ END FUNCTION interp_0_log2
+
 END MODULE module_init_utilities