Changeset 1001 for LMDZ4/trunk/libf/phylmd/physiq.F

Timestamp:

Oct 6, 2008, 11:11:53 AM (16 years ago)

Author:

Laurent Fairhead

Message:

• Modifs sur le parallelisme: masquage dans la physique
• Inclusion strato
• mise en coherence etat0
• le mode offline fonctionne maintenant en parallele,
• les fichiers de la dynamiques sont correctement sortis et peuvent etre reconstruit avec rebuild
• la version parallele de la dynamique peut s'executer sans MPI (sur 1 proc)
• L'OPENMP fonctionne maintenant sans la parallelisation MPI.

YM
LF

File:

• LMDZ4/trunk/libf/phylmd/physiq.F (modified) (8 diffs)

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -94 +94 @@
 #include "clesphys.h"
 #include "control.h"
-#include "logic.h"
+!#include "logic.h"
 #include "temps.h"
 cym#include "comgeomphy.h"
@@ -943 +943 @@
       REAL zustrli(klon), zvstrli(klon)
       REAL zustrph(klon), zvstrph(klon)
+      REAL zustrhi(klon), zvstrhi(klon)
       REAL aam, torsfc
 cIM 141004 END
@@ -1377 +1378 @@
 !           ENDDO
 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
-
-           CALL SUGWD(klon,klev,paprs,pplay)
+           IF (ok_strato) THEN
+             CALL SUGWD_strato(klon,klev,paprs,pplay)
+           ELSE
+             CALL SUGWD(klon,klev,paprs,pplay)
+           ENDIF
+           
            DO i=1,klon
              zuthe(i)=0.
@@ -2932 +2937 @@
 c        igwdim=MAX(1,igwd)
 c
+        IF (ok_strato) THEN
+        
+          CALL drag_noro_strato(klon,klev,dtime,paprs,pplay,
+     e                   zmea,zstd, zsig, zgam, zthe,zpic,zval,
+     e                   igwd,idx,itest,
+     e                   t_seri, u_seri, v_seri,
+     s                   zulow, zvlow, zustrdr, zvstrdr,
+     s                   d_t_oro, d_u_oro, d_v_oro)
+
+       ELSE
         CALL drag_noro(klon,klev,dtime,paprs,pplay,
      e                   zmea,zstd, zsig, zgam, zthe,zpic,zval,
      e                   igwd,idx,itest,
      e                   t_seri, u_seri, v_seri,
-cIM 141004    s                   zulow, zvlow, zustr, zvstr,
      s                   zulow, zvlow, zustrdr, zvstrdr,
      s                   d_t_oro, d_u_oro, d_v_oro)
+       ENDIF
 c
 c  ajout des tendances
@@ -2969 +2984 @@
 c        igwdim=MAX(1,igwd)
 c
-        CALL lift_noro(klon,klev,dtime,paprs,pplay,
+        IF (ok_strato) THEN
+
+          CALL lift_noro_strato(klon,klev,dtime,paprs,pplay,
+     e                   rlat,zmea,zstd,zpic,zgam,zthe,zpic,zval,
+     e                   igwd,idx,itest,
+     e                   t_seri, u_seri, v_seri,
+     s                   zulow, zvlow, zustrli, zvstrli,
+     s                   d_t_lif, d_u_lif, d_v_lif               )
+        
+        ELSE
+          CALL lift_noro(klon,klev,dtime,paprs,pplay,
      e                   rlat,zmea,zstd,zpic,
      e                   itest,
@@ -2975 +3000 @@
      s                   zulow, zvlow, zustrli, zvstrli,
      s                   d_t_lif, d_u_lif, d_v_lif)
-c
+       ENDIF
+c   
 !-----------------------------------------------------------------------------------------
 ! ajout des tendances de la portance de l'orographie
@@ -2982 +3008 @@
 c
       ENDIF ! fin de test sur ok_orolf
+C  HINES GWD PARAMETRIZATION
+
+       IF (ok_hines) then
+
+         CALL hines_gwd(klon,klev,dtime,paprs,pplay,
+     i                  rlat,t_seri,u_seri,v_seri,
+     o                  zustrhi,zvstrhi,
+     o                  d_t_hin, d_u_hin, d_v_hin)
+c
+c  ajout des tendances
+        CALL add_phys_tend(d_u_hin,d_v_hin,d_t_hin,dq0,dql0,'lif')
+
+      ENDIF
+c
+
 c
 cIM cf. FLott BEG
@@ -3368 +3409 @@
          itau_phy = itau_phy + itap
          CALL phyredem ("restartphy.nc")
-         open(97,form="unformatted",file="finbin")
-         write(97) u_seri,v_seri,t_seri,q_seri
-         close(97)
+!         open(97,form="unformatted",file="finbin")
+!         write(97) u_seri,v_seri,t_seri,q_seri
+!         close(97)
       ENDIF