Changeset 3780 for LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_interp_cases.h

Timestamp:

Oct 22, 2020, 2:50:18 PM (4 years ago)

Author:

evignon

Message:

Premiere comission Etienne: changements pour le 1D (forcage en Ts au dessus des continents) et inclusion drag arbres dans yamada4_num=6

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_interp_cases.h (modified) (7 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/old_1D_interp_cases.h

@@ -62 +62 @@
      &             ,ts_prof,plev_prof,t_prof,q_prof,u_prof,v_prof,w_prof    &
      &             ,ht_prof,vt_prof,hq_prof,vq_prof)
-
-        if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
+! EV: tg instead of ts_cur
+        if (type_ts_forcing.eq.1) tg = ts_prof ! 
 
 ! vertical interpolation:
@@ -113 +113 @@
 !     print *,'llm l omega_profd',llm,l,omega_profd(l)
 !     enddo
-
-        if (type_ts_forcing.eq.1) ts_cur = tg_prof ! SST used in read_tsurf1d
+! EV tg instead of ts_cur
+        if (type_ts_forcing.eq.1) tg = tg_prof ! SST used 
 
 ! vertical interpolation:
@@ -206 +206 @@
      &             ,ug_gabls4,vg_gabls4,ht_gabls4,hq_gabls4,tg_gabls4                            &
      &             ,ug_profg,vg_profg,ht_profg,hq_profg,tg_profg)
-
-        if (type_ts_forcing.eq.1) ts_cur = tg_prof ! SST used in read_tsurf1d
+!EV tg instead of ts_cur
+        if (type_ts_forcing.eq.1) tg = tg_prof ! SST used
 
 ! vertical interpolation:
@@ -499 +499 @@
      &             ,nlev_sandu                                              &
      &             ,ts_sandu,ts_prof)
-
-        if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
+! EV tg instead of ts_cur
+        if (type_ts_forcing.eq.1) tg = ts_prof ! SST used in read_tsurf1d
 
 ! vertical interpolation:
@@ -582 +582 @@
      &             ,ufa_astex,vfa_astex,div_prof,ts_prof,ug_prof,vg_prof    &
      &             ,ufa_prof,vfa_prof)
-
-        if (type_ts_forcing.eq.1) ts_cur = ts_prof ! SST used in read_tsurf1d
-
+! EV tg instead of ts_cur
+        if (type_ts_forcing.eq.1) tg = ts_prof ! SST used 
 ! vertical interpolation:
       CALL interp_astex_vertical(play,nlev_astex,plev_profa                 &
@@ -675 +674 @@
      &       ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas,lat_prof_cas               &
      &       ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas)
-
-             ts_cur = ts_prof_cas 
+! EV tg instead of ts_cur
+
+             tg = ts_prof_cas 
              psurf=plev_prof_cas(1)
 
@@ -850 +850 @@
      &       ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas                           &
      &       ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas)
-
-             ts_cur = ts_prof_cas 
+! EV tg instead of ts_cur
+
+             tg = ts_prof_cas 
 !            psurf=plev_prof_cas(1)
              psurf=ps_prof_cas