Changeset 1848

Timestamp:

Dec 10, 2017, 10:07:25 PM (7 years ago)

Author:

aslmd

Message:

idealized mesoscale, implemented the Wee et al. 2012 bias correction

File:

• trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_quarter_ss.F (modified) (8 diffs)

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

@@ -101 +101 @@
 !****Mars
    INTEGER, PARAMETER :: nl_max = 1000
-   REAL, DIMENSION(nl_max) :: zk, p_in, theta, rho, u, v, qv, pd_in
+   REAL, DIMENSION(nl_max) :: zk, p_in, theta, tk, rho, u, v, qv, pd_in
    INTEGER :: nl_in
 
@@ -354 +354 @@
   write(6,*) ' getting dry sounding for base state '
 
-  CALL get_sounding( zk, p_in, pd_in, theta, rho, u, v, qv, dry_sounding, nl_max, nl_in, &
+  CALL get_sounding( zk, p_in, pd_in, theta, tk, rho, u, v, qv, dry_sounding, nl_max, nl_in, &
                         mulu, mulv, addu, addv )
   ENDIF
@@ -361 +361 @@
   CALL wrf_dm_bcast_real( pd_in , nl_max )
   CALL wrf_dm_bcast_real( theta , nl_max )
+  CALL wrf_dm_bcast_real( tk , nl_max )
   CALL wrf_dm_bcast_real( rho , nl_max )
   CALL wrf_dm_bcast_real( u , nl_max )
@@ -538 +539 @@
       p_level = grid%em_znu(k)*(p_surf - grid%p_top) + grid%p_top
       grid%em_pb(i,k,j) = p_level
-      grid%em_t_init(i,k,j) = interp_0( theta, p_in, p_level, nl_in ) - t0
+
+! OLD METHOD
+!      grid%em_t_init(i,k,j) = interp_0( theta, p_in, p_level, nl_in ) - t0
+! NEW METHOD: Wee et al. 2012
+! interpolate temperature. then convert to potential temperature.
+      grid%em_t_init(i,k,j) = interp_0( tk, p_in, p_level, nl_in )
+      !! l un ou l autre pareil
+      grid%em_t_init(i,k,j) = - t0 + (grid%em_t_init(i,k,j) * ((p1000mb/p_level)**rcp))
+      !grid%em_t_init(i,k,j) = - t0 + (grid%em_t_init(i,k,j) * ((610./p_level)**(1.0/3.9)))
       grid%em_alb(i,k,j) = (r_d/p1000mb)*(grid%em_t_init(i,k,j)+t0)*(grid%em_pb(i,k,j)/p1000mb)**cvpm
     ENDDO
@@ -574 +583 @@
 !  dry_sounding = .false.
   dry_sounding = .true.
-  CALL get_sounding( zk, p_in, pd_in, theta, rho, u, v, qv, dry_sounding, nl_max, nl_in, &
+  CALL get_sounding( zk, p_in, pd_in, theta, tk, rho, u, v, qv, dry_sounding, nl_max, nl_in, &
                               mulu, mulv, addu, addv )
 
@@ -922 +931 @@
 !---------------------------------------------------------------------------
 
-      subroutine get_sounding( zk, p, p_dry, theta, rho, &
+      subroutine get_sounding( zk, p, p_dry, theta, tk, rho, &
                                u, v, qv, dry, nl_max, nl_in, &
                                mulu, mulv, addu, addv )
@@ -928 +937 @@
 
       integer nl_max, nl_in
-      real zk(nl_max), p(nl_max), theta(nl_max), rho(nl_max), &
+      real zk(nl_max), p(nl_max), theta(nl_max), tk(nl_max), rho(nl_max), &
            u(nl_max), v(nl_max), qv(nl_max), p_dry(nl_max)
       logical dry
@@ -1095 +1104 @@
           p_dry(k) = p_input(k)
           theta(k) = th_input(k)
+          tk(k) = t_input(k)
           rho(k) = rho_input(k)
           u(k) = mulu*u_input(k) + addu