Changeset 3666 for LMDZ6/trunk/libf/phylmd/physiq_mod.F90

Timestamp:

Apr 20, 2020, 12:13:34 PM (4 years ago)

Author:

lfalletti

Message:

Adding changes for Reprobus

File:

• LMDZ6/trunk/libf/phylmd/physiq_mod.F90 (modified) (6 diffs)

LMDZ6/trunk/libf/phylmd/physiq_mod.F90

@@ -246 +246 @@
 #endif
 #ifdef REPROBUS
-    USE CHEM_REP, ONLY : Init_chem_rep_xjour
+    USE CHEM_REP, ONLY : Init_chem_rep_xjour, &
+         d_q_rep,d_ql_rep,d_qi_rep,ptrop,ttrop, &
+         ztrop, gravit,itroprep, Z1,Z2,fac,B
 #endif
     USE indice_sol_mod
@@ -1950 +1952 @@
 #endif
        ENDIF
+       IF (type_trac == 'repr') THEN
+#ifdef REPROBUS
+          CALL chemini_rep(  &
+               presnivs, &
+               pdtphys, &
+               annee_ref, &
+               day_ref,  &
+               day_ini, &
+               start_time, &
+               itau_phy, &
+               io_lon, &
+               io_lat)
+#endif
+       ENDIF
 
        !$omp single
@@ -2286 +2302 @@
 
           wo(:,:,1)=ozonecm(latitude_deg, paprs,read_climoz,rjour=zzz)
+#ifdef REPROBUS
+          ptrop=dyn_tropopause(t_seri, ztsol, paprs, pplay, rot)/100.
+          DO i = 1, klon
+             Z1=t_seri(i,itroprep(i)+1)
+             Z2=t_seri(i,itroprep(i))
+             fac=(Z1-Z2)/alog(pplay(i,itroprep(i)+1)/pplay(i,itroprep(i)))
+             B=Z2-fac*alog(pplay(i,itroprep(i)))
+             ttrop(i)= fac*alog(ptrop(i))+B
+!        
+             Z1= 1.e-3 * ( pphi(i,itroprep(i)+1)+pphis(i) ) / gravit
+             Z2= 1.e-3 * ( pphi(i,itroprep(i))  +pphis(i) ) / gravit
+             fac=(Z1-Z2)/alog(pplay(i,itroprep(i)+1)/pplay(i,itroprep(i)))
+             B=Z2-fac*alog(pplay(i,itroprep(i)))
+             ztrop(i)=fac*alog(ptrop(i))+B
+          ENDDO
+#endif
        ELSE
           !--- ro3i = elapsed days number since current year 1st january, 0h
@@ -3754 +3786 @@
 #endif 
     ENDIF !type_trac = inca
-
+    IF (type_trac == 'repr') THEN
+#ifdef REPROBUS
+    !CALL chemtime_rep(itap+itau_phy-1, date0, dtime, itap)
+    CALL chemtime_rep(itap+itau_phy-1, date0, phys_tstep, itap) 
+#endif
+    ENDIF
 
     !
@@ -4719 +4756 @@
 
     IF (type_trac=='repr') THEN
+!MM pas d'impact, car on recupere q_seri,tr_seri,t_seri via phys_local_var_mod
+!MM                               dans Reprobus
        sh_in(:,:) = q_seri(:,:)
+#ifdef REPROBUS
+       d_q_rep(:,:) = 0.
+       d_ql_rep(:,:) = 0.
+       d_qi_rep(:,:) = 0.
+#endif
     ELSE
        sh_in(:,:) = qx(:,:,ivap)
@@ -4769 +4813 @@
          d_tr_dyn, &                                 !<<RomP
          tr_seri, init_source)
+#ifdef REPROBUS
+
+
+          print*,'avt add phys rep',abortphy
+
+     CALL add_phys_tend &
+            (du0,dv0,dt0,d_q_rep,d_ql_rep,d_qi_rep,paprs,&
+             'rep',abortphy,flag_inhib_tend,itap,0)
+        IF (abortphy==1) Print*,'ERROR ABORT REP'
+
+          print*,'apr add phys rep',abortphy
+
+#endif
+
 #endif
     ENDIF    ! (iflag_phytrac=1)