Changeset 2333 for LMDZ5/trunk/libf/phylmd/physiq.F90

Timestamp:

Jul 22, 2015, 4:14:51 PM (9 years ago)

Author:

lguez

Message:

New parameterization of gravity wave drag due to front/jet systems, by

1. de la Camara and F. Lott. The new Camara-Lott parameterization

replaces the Hines parameterization so it is activated if not ok_hines
and ok_gwd_rando.

Also changed distribution of phase speeds in FLOTT_GWD_rando, from
uniform to Gaussian. Bug fix in sugwd_strato. Bug fix in the arguments
of the call to add_phys_tend for methane oxydation.

For the new Camara-Lott parameterization, we need to compute relative
vorticity in calfis and pass it as a new argument "rot" to
physiq. Interpolation of relative vorticity to the physics grid is not
optimal for now: it is not weighted by cell areas.

Alvaro de la Camara, Fran\c{}cois Lott

File:

• LMDZ5/trunk/libf/phylmd/physiq.F90 (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -5 +5 @@
      debut,lafin,jD_cur, jH_cur,pdtphys, &
      paprs,pplay,pphi,pphis,presnivs, &
-     u,v,t,qx, &
+     u,v,rot,t,qx, &
      flxmass_w, &
      d_u, d_v, d_t, d_qx, d_ps &
@@ -66 +66 @@
   USE CFMIP_point_locations
   use FLOTT_GWD_rando_m, only: FLOTT_GWD_rando
+  use ACAMA_GWD_rando_m, only: ACAMA_GWD_rando
 
   IMPLICIT none
@@ -230 +231 @@
   REAL u(klon,klev)
   REAL v(klon,klev)
+
+  REAL, intent(in):: rot(klon, klev)
+  ! relative vorticity, in s-1, needed for frontal waves
+
   REAL t(klon,klev),thetal(klon,klev)
   ! thetal: ligne suivante a decommenter si vous avez les fichiers     MPL 20130625
@@ -696 +701 @@
   REAL zustrli(klon), zvstrli(klon)
   REAL zustrph(klon), zvstrph(klon)
-  REAL zustrhi(klon), zvstrhi(klon)
   REAL aam, torsfc
   !IM 141004 END
@@ -3835 +3839 @@
              d_t_lif, d_u_lif, d_v_lif)
      ENDIF
-     !   
-     !-----------------------------------------------------------------------------------------
+
      ! ajout des tendances de la portance de l'orographie
-     CALL add_phys_tend(d_u_lif,d_v_lif,d_t_lif,dq0,dql0,dqi0,paprs,'lif',abortphy)
-     !-----------------------------------------------------------------------------------------
-     !
+     CALL add_phys_tend(d_u_lif, d_v_lif, d_t_lif, dq0, dql0, dqi0, paprs, &
+          'lif', abortphy)
   ENDIF ! fin de test sur ok_orolf
-  !  HINES GWD PARAMETRIZATION
 
   IF (ok_hines) then
-
-     CALL hines_gwd(klon,klev,dtime,paprs,pplay, &
-          rlat,t_seri,u_seri,v_seri, &
-          zustrhi,zvstrhi, &
-          d_t_hin, d_u_hin, d_v_hin)
-     !
-     !  ajout des tendances
-     CALL add_phys_tend(d_u_hin,d_v_hin,d_t_hin,dq0,dql0,dqi0,paprs,'hin',abortphy)
-
+     !  HINES GWD PARAMETRIZATION
+     east_gwstress=0.
+     west_gwstress=0.
+     du_gwd_hines=0.
+     dv_gwd_hines=0.
+     CALL hines_gwd(klon, klev, dtime, paprs, pplay, rlat, t_seri, u_seri, &
+          v_seri, zustr_gwd_hines, zvstr_gwd_hines, d_t_hin, du_gwd_hines, &
+          dv_gwd_hines)
+     zustr_gwd_hines=0.
+     zvstr_gwd_hines=0.
+     DO k = 1, klev
+        zustr_gwd_hines(:)=zustr_gwd_hines(:)+ du_gwd_hines(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+        zvstr_gwd_hines(:)=zvstr_gwd_hines(:)+ dv_gwd_hines(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+     ENDDO
+
+     d_t_hin(:, :)=0.
+     CALL add_phys_tend(du_gwd_hines, dv_gwd_hines, d_t_hin, dq0, dql0, dqi0, &
+          paprs, 'hin', abortphy)
+  ENDIF
+
+  IF (.not. ok_hines .and. ok_gwd_rando) then
+     CALL acama_GWD_rando(DTIME, pplay, rlat, t_seri, u_seri, v_seri, rot, &
+          zustr_gwd_front, zvstr_gwd_front, du_gwd_front, dv_gwd_front, &
+          east_gwstress, west_gwstress)
+     zustr_gwd_front=0.
+     zvstr_gwd_front=0.
+     DO k = 1, klev
+        zustr_gwd_front(:)=zustr_gwd_front(:)+ du_gwd_front(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+        zvstr_gwd_front(:)=zvstr_gwd_front(:)+ dv_gwd_front(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+     ENDDO
+
+     CALL add_phys_tend(du_gwd_front, dv_gwd_front, dt0, dq0, dql0, dqi0, &
+          paprs, 'front_gwd_rando', abortphy)
+
+     !  LIGNE TEMPORAIRE POUR TESTER ROT (ECRITE AUSSI AU DESSOUS!!!):
+     east_gwstress=rot
   ENDIF
 
@@ -3859 +3891 @@
      call FLOTT_GWD_rando(DTIME, pplay, t_seri, u_seri, v_seri, &
           rain_fall + snow_fall, zustr_gwd_rando, zvstr_gwd_rando, &
-          du_gwd_rando, dv_gwd_rando)
-     CALL add_phys_tend(du_gwd_rando, dv_gwd_rando, dt0, dq0, dql0,dqi0,paprs, &
-          'flott_gwd_rando',abortphy)
+          du_gwd_rando, dv_gwd_rando, east_gwstress, west_gwstress)
+     CALL add_phys_tend(du_gwd_rando, dv_gwd_rando, dt0, dq0, dql0, dqi0, &
+          paprs, 'flott_gwd_rando', abortphy)
+     zustr_gwd_rando=0.
+     zvstr_gwd_rando=0.
+     DO k = 1, klev
+        zustr_gwd_rando(:)=zustr_gwd_rando(:)+ du_gwd_rando(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+        zvstr_gwd_rando(:)=zvstr_gwd_rando(:)+ dv_gwd_rando(:, k)/dtime &
+             * (paprs(:, k)-paprs(:, k+1))/rg
+     ENDDO
+
+     !  LIGNE TEMPORAIRE POUR TESTER ROT (ECRITE AUSSI AU DESSUS!!!):
+     east_gwstress=rot
   end if
 
@@ -3915 +3958 @@
      CALL METHOX(1,klon,klon,klev,q_seri,d_q_ch4,pplay)
   ! ajout de la tendance d'humidite due au methane
-     CALL add_phys_tend(du0,dv0,dt0,d_q_ch4*dtime,dql0,'q_ch4',abortphy)
+     CALL add_phys_tend(du0, dv0, dt0, d_q_ch4*dtime, dql0, dqi0, paprs, &
+          'q_ch4', abortphy)
   END IF
   !