Changeset 1378 for trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars

Timestamp:

Feb 18, 2015, 3:36:22 PM (10 years ago)

Author:

aslmd

Message:

consistent computations to initialize tracers for LES

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F (modified) (10 diffs)

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

@@ -115 +115 @@
  REAL :: lon_input, lat_input, alt_input, tsurf_input
  ! for mode 3
- REAL, DIMENSION(nl_max) :: profdustq,profdustn,dust_p_level
+ REAL, DIMENSION(nl_max) :: profdustq,profdustn
 !!MARS
 
@@ -757 +757 @@
   ENDDO
 
-!!!!! MARS INITIALIZE MODE 3 FOR LES
-!! RESTART=config_flags%restart -> this routine is not used in restart
-
-  IF ( ( config_flags%mars == 3 ) ) then
-
-  write (*,*) 'MARS MODE 3 INITIALIZATION, READING INPUT_DUST (module_initialize_les.F)'
-    ! load a profile of dust (same for all points)
-    call read_dust(profdustq,profdustn,dust_p_level,nl_in)
-    p_level = grid%znu(1)*(pd_surf - grid%p_top) + grid%p_top
-      IF (dust_p_level(1) .lt. p_level) then !input profile needs rescaling to avoid a plateau. This happens when you use different sources to initialize the pressure of the LES and the input_dust. Typicaly: you use the MCD for input_therm and different runs for input_dust, that dont have the same surface pressure ! This trick is ok because we dont want to initialize with a very precise profile of dust, just a realistic one.
-         dust_p_level = dust_p_level * p_level/dust_p_level(1)
-     ENDIF
-     DO k=1,kte-1
+!!!!! MARS 
+
+    ! interpolate water vapor
+    if (      ( config_flags%mars == 1  ) &
+         .OR. ( config_flags%mars == 11 ) &
+         .OR. ( config_flags%mars == 12 ) ) then
+      print *, '**** INTERPOLATE HV **** RANK 2 in SCALAR'
+      DO k=1,kte-1
          p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
-         print *, '**** INTERPOLATE DUSTQ **** RANK 2 in SCALAR'
-         scalar(its:ite,k,jts:jte,2) = interp_0( profdustq, dust_p_level, p_level, nl_in )
-         print *, '**** INTERPOLATE DUSTN **** RANK 3 in SCALAR'
-         scalar(its:ite,k,jts:jte,3) = interp_0( profdustn, dust_p_level, p_level, nl_in ) 
-     ENDDO
-  ENDIF
+         scalar(its:ite,k,jts:jte,2) = interp_0( qv, pd_in, p_level, nl_in )
+         scalar(its:ite,k,jts:jte,3) = 0.
+           !! water ice is set to 0 (was put into water vapor when building prof from MCD)
+      ENDDO
+      print *, "WATER VAPOR",scalar(its,:,jts,2)
+    endif
+
+    ! interpolate qdust
+    if (      ( config_flags%mars == 11 ) &
+         .OR. ( config_flags%mars == 12 ) ) then
+      call read_dust(profdustq,profdustn,nl_in)
+      print *, '**** INTERPOLATE DUSTQ **** RANK 4 in SCALAR'
+      print *, '**** INTERPOLATE DUSTN **** RANK 5 in SCALAR'
+      DO k=1,kte-1
+         p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
+         scalar(its:ite,k,jts:jte,4) = interp_0( profdustq, pd_in, p_level, nl_in )
+         scalar(its:ite,k,jts:jte,5) = interp_0( profdustn, pd_in, p_level, nl_in )
+      ENDDO
+      print *, "DUST Q", scalar(its,:,jts,4)
+      print *, "DUST N", scalar(its,:,jts,5)
+    endif
+
+    if ( config_flags%mars == 12 ) then
+      scalar(its:ite,1:kte-1,jts:jte,6) = 0.
+      scalar(its:ite,1:kte-1,jts:jte,7) = 0.
+    endif
 
 !!!!! MARS 
@@ -866 +881 @@
 ! -----
 
-      if(dry) then
-       do k=1,nl
-         qv_input(k) = 0.
-       enddo
-      endif
+      !if(dry) then
+      ! do k=1,nl
+      !   qv_input(k) = 0.
+      ! enddo
+      !endif
 
       if(debug) write(6,*) ' number of input levels = ',nl
@@ -884 +899 @@
 
       do k=1,nl
+        !!!!!!!!!!!!!! MARS
+        !! from mol/mol to kg/kg
+        qv_input(k) = qv_input(k)*18./mwdry
         qv_input(k) = 0.001*qv_input(k)
       enddo
@@ -889 +907 @@
       p_surf = 100.*p_surf  ! convert to pascals
       qvf = 1. + rvovrd*qv_input(1) 
+!!MARS
+      qvf = 1.
       rho_surf = 1./((r/p1000mb)*th_surf*qvf*((p_surf/p1000mb)**cvpm))
       pi_surf = (p_surf/p1000mb)**(r/cp)
@@ -904 +924 @@
           qvf = 1. + rvovrd*qv_input(1) 
           qvf1 = 1. + qv_input(1)
+!!MARS
+          qvf = 1.
+          qvf1 = 1.
           rho_input(1) = rho_surf
           dz = h_input(1)
@@ -921 +944 @@
             qvf1 = 0.5*(2.+(qv_input(k-1)+qv_input(k)))
             qvf = 1. + rvovrd*qv_input(k)   ! qv is in g/kg here
- 
+!!MARS
+          qvf = 1.
+          qvf1 = 1.
+!!MARS 
             do it=1,10
               pm_input(k) = pm_input(k-1) &
@@ -1035 +1061 @@
 
 !!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-      subroutine read_dust(pdustq,pdustn,dust_p_level,n)
+      subroutine read_dust(pdustq,pdustn,n)
       implicit none
       integer n
-      real pdustq(n),pdustn(n),dust_p_level(n)
+      real pdustq(n),pdustn(n)
       logical end_of_file
 
@@ -1051 +1077 @@
       do while (.not. end_of_file)
 
-        read(11,*,end=102) pdustq(k+1),pdustn(k+1), dust_p_level(k+1)
-        write(*,*) k, pdustq(k+1), pdustn(k+1), dust_p_level(k+1)
+        read(11,*,end=102) pdustq(k+1),pdustn(k+1)
+        write(*,*) k,pdustq(k+1),pdustn(k+1)
         k = k+1
         go to 113
@@ -1064 +1090 @@
 !!cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 
-
 END MODULE module_initialize_ideal