Changeset 3686 for LMDZ6/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

Timestamp:

May 27, 2020, 2:59:10 PM (4 years ago)

Author:

fhourdin

Message:

Retropedalage sur le 1D.
Annulation de toutes les modifs de la veille.
Frederic

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/1D_read_forc_cases.h

@@ -17 +17 @@
          call read_SCM_cas
          write(*,*) 'Forcing read' 
-         print*,'PS ps_cas',ps_cas
+         print*,'PS ps_cas',ps_cas
 
 !Time interpolation for initial conditions using interpolation routine
          write(*,*) 'AVT 1ere INTERPOLATION: day,day1 = ',daytime,day1   
-        CALL interp_case_time_std(day,day1,annee_ref                                 &
-     &         ,nt_cas,nlev_cas                                                      &
-     &         ,ts_cas,ps_cas,ps_forc_cas,plev_cas,ppforc_cas,t_cas,theta_cas,thl_cas&
-     &         ,qt0_cas,qv0_cas,ql0_cas,qi0_cas                                      &
-     &         ,rt0_cas,rv0_cas,rl0_cas,ri0_cas,rh0_cas                              &
-     &         ,u_cas,v_cas,w_cas,omega_cas,ug_cas,vg_cas                            &
-     &         ,temp_nudg_cas,th_nudg_cas,thl_nudg_cas,qt_nudg_cas,qv_nudg_cas       &
-     &         ,rt_nudg_cas,rv_nudg_cas,u_nudg_cas,v_nudg_cas                        &
-     &         ,uadv_cas,vadv_cas,tadv_cas,thadv_cas,thladv_cas                      &
-     &         ,qtadv_cas,qvadv_cas,rtadv_cas,rvadv_cas                              &
-     &         ,trad_cas,thrad_cas,thlrad_cas                                        &
-     &         ,tke_cas,lat_cas,sens_cas,ustar_cas                                   &
-     &         ,wpthetap_cas,wpqtp_cas,wpqvp_cas,wprtp_cas,wprvp_cas                 &
+        CALL interp_case_time_std(daytime,day1,annee_ref                                       &
+!    &       ,year_ini_cas,day_ju_ini_cas,nt_cas,pdt_cas,nlev_cas                           &
+     &       ,nt_cas,nlev_cas                                                               &
+     &       ,ts_cas,ps_cas,plev_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,ql_cas,qi_cas      &
+     &       ,u_cas,v_cas,ug_cas,vg_cas                                                     &
+     &       ,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas                               &
+     &       ,vitw_cas,omega_cas,du_cas,hu_cas,vu_cas                                       &
+     &       ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas                           &
+     &       ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,lat_cas,sens_cas,ustar_cas       &
+     &       ,uw_cas,vw_cas,q1_cas,q2_cas,tke_cas                                           &
 !
-     &         ,ts_prof_cas,ps_prof_cas,ps_forc_prof_cas,plev_prof_cas,pforc_prof_cas&
-     &         ,t_prof_cas,th_prof_cas,thl_prof_cas                               &
-     &         ,qt_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas                      &
-     &         ,rt_prof_cas,rv_prof_cas,rl_prof_cas,ri_prof_cas,rh_prof_cas          &
-     &         ,u_prof_cas,v_prof_cas,w_prof_cas,omega_prof_cas                      &
-     &         ,ug_prof_cas,vg_prof_cas                                              &
-     &         ,temp_nudg_prof_cas,th_nudg_prof_cas,thl_nudg_prof_cas                &
-     &         ,qt_nudg_prof_cas,qv_nudg_prof_cas                                    &
-     &         ,rt_nudg_prof_cas,rv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas    &
-     &         ,uadv_prof_cas,vadv_prof_cas,tadv_prof_cas,thadv_prof_cas,thladv_prof_cas&
-     &         ,qtadv_prof_cas,qvadv_prof_cas,rtadv_prof_cas,rvadv_prof_cas          &
-     &         ,trad_prof_cas,thrad_prof_cas,thlrad_prof_cas                         &
-     &         ,tke_prof_cas,lat_prof_cas,sens_prof_cas,ustar_prof_cas               &
-     &         ,wpthetap_prof_cas,wpqtp_prof_cas,wpqvp_prof_cas,wprtp_prof_cas,wprvp_prof_cas)
+     &       ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas,thv_prof_cas  &
+     &       ,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas                              &
+     &       ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas                                 &
+     &       ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas           &
+     &       ,vitw_prof_cas,omega_prof_cas                                                  &
+     &       ,du_prof_cas,hu_prof_cas,vu_prof_cas                                           &
+     &       ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas       &
+     &       ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas                            &
+     &       ,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)
 
       do l = 1, nlev_cas
@@ -56 +49 @@
 ! vertical interpolation using interpolation routine:
 !      write(*,*)'avant interp vert', t_prof
-       CALL interp_case_vertical_std(nlev_cas                                                          & 
-     &         ,plev_prof_cas,t_prof_cas,th_prof_cas,thl_prof_cas                                   &
-     &         ,qt_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas                                        &
-     &         ,rt_prof_cas,rv_prof_cas,rl_prof_cas,ri_prof_cas,rh_prof_cas                            &
-     &         ,u_prof_cas,v_prof_cas,w_prof_cas,omega_prof_cas                                        &
-     &         ,ug_prof_cas,vg_prof_cas                                                                &
-     &         ,temp_nudg_prof_cas,th_nudg_prof_cas,thl_nudg_prof_cas                                  &
-     &         ,qt_nudg_prof_cas,qv_nudg_prof_cas                                                      &
-     &         ,rt_nudg_prof_cas,rv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas                      &
-     &         ,uadv_prof_cas,vadv_prof_cas,tadv_prof_cas,thadv_prof_cas,thladv_prof_cas               &
-     &         ,qtadv_prof_cas,qvadv_prof_cas,rtadv_prof_cas,rvadv_prof_cas                            &
-     &         ,trad_prof_cas,thrad_prof_cas,thlrad_prof_cas                                           &
+      CALL interp2_case_vertical_std(play,nlev_cas,plev_prof_cas                                              &
+     &         ,t_prof_cas,theta_prof_cas,thv_prof_cas,thl_prof_cas                                          &
+     &         ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas                                 &
+     &         ,ug_prof_cas,vg_prof_cas                                                                   &
+     &       ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas                         &
+
+     &         ,vitw_prof_cas,omega_prof_cas                                                              &
+     &         ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas                   &
+     &         ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas    &
+     &         ,dth_prof_cas,hth_prof_cas,vth_prof_cas                                                    &
 !
-     &         ,plev_mod_cas,t_mod_cas,theta_mod_cas,thl_mod_cas                                       &
-     &         ,qt_mod_cas,qv_mod_cas,ql_mod_cas,qi_mod_cas                                            &
-     &         ,rt_mod_cas,rv_mod_cas,rl_mod_cas,ri_mod_cas,rh_mod_cas                                 &
-     &         ,u_mod_cas,v_mod_cas,w_mod_cas,omega_mod_cas                                            &
-     &         ,ug_mod_cas,vg_mod_cas                                                                  &
-     &         ,temp_nudg_mod_cas,th_nudg_mod_cas,thl_nudg_mod_cas                                     &
-     &         ,qt_nudg_mod_cas,qv_nudg_mod_cas                                                        &
-     &         ,rt_nudg_mod_cas,rv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas                          &
-     &         ,uadv_mod_cas,vadv_mod_cas,tadv_mod_cas,thadv_mod_cas,thladv_mod_cas                    &
-     &         ,qtadv_mod_cas,qvadv_mod_cas,rtadv_mod_cas,rvadv_mod_cas                                &
-     &         ,trad_mod_cas,thrad_mod_cas,thlrad_mod_cas)
+     &         ,t_mod_cas,theta_mod_cas,thv_mod_cas,thl_mod_cas,qv_mod_cas,ql_mod_cas,qi_mod_cas          &
+     &         ,u_mod_cas,v_mod_cas,ug_mod_cas,vg_mod_cas                                                 &
+     &         ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas                           &
+     &         ,w_mod_cas,omega_mod_cas                                                                   &
+     &         ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas                         &
+     &         ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas           &
+     &         ,dth_mod_cas,hth_mod_cas,vth_mod_cas,mxcalc)