Changeset 3680 for LMDZ6/trunk/libf/phylmd/dyn1d/1D_interp_cases.h
- Timestamp:
- May 26, 2020, 8:10:45 AM (4 years ago)
- File:
-
- 1 edited
-
LMDZ6/trunk/libf/phylmd/dyn1d/1D_interp_cases.h (modified) (6 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/dyn1d/1D_interp_cases.h
r3594 r3680 7 7 8 8 ! time interpolation: 9 CALL interp_case_time_std(daytime,day1,annee_ref & 10 ! & ,year_ini_cas,day_ju_ini_cas,nt_cas,pdt_cas,nlev_cas & 11 & ,nt_cas,nlev_cas & 12 & ,ts_cas,ps_cas,plev_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,ql_cas,qi_cas & 13 & ,u_cas,v_cas,ug_cas,vg_cas & 14 & ,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas & 15 & ,vitw_cas,omega_cas,du_cas,hu_cas,vu_cas & 16 & ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas & 17 & ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,lat_cas,sens_cas,ustar_cas & 18 & ,uw_cas,vw_cas,q1_cas,q2_cas,tke_cas & 9 CALL interp_case_time_std(day,day1,annee_ref & 10 & ,nt_cas,nlev_cas & 11 & ,ts_cas,ps_cas,ps_forc_cas,plev_cas,ppforc_cas,t_cas,theta_cas,thl_cas& 12 & ,qt0_cas,qv0_cas,ql0_cas,qi0_cas & 13 & ,rt0_cas,rv0_cas,rl0_cas,ri0_cas,rh0_cas & 14 & ,u_cas,v_cas,w_cas,omega_cas,ug_cas,vg_cas & 15 & ,temp_nudg_cas,th_nudg_cas,thl_nudg_cas,qt_nudg_cas,qv_nudg_cas & 16 & ,rt_nudg_cas,rv_nudg_cas,u_nudg_cas,v_nudg_cas & 17 & ,uadv_cas,vadv_cas,tadv_cas,thadv_cas,thladv_cas & 18 & ,qtadv_cas,qvadv_cas,rtadv_cas,rvadv_cas & 19 & ,trad_cas,thrad_cas,thlrad_cas & 20 & ,tke_cas,lat_cas,sens_cas,ustar_cas & 21 & ,wpthetap_cas,wpqtp_cas,wpqvp_cas,wprtp_cas,wprvp_cas & 19 22 ! 20 & ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas,thv_prof_cas & 21 & ,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas & 22 & ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas & 23 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 24 & ,vitw_prof_cas,omega_prof_cas & 25 & ,du_prof_cas,hu_prof_cas,vu_prof_cas & 26 & ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas & 27 & ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & 28 & ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas & 29 & ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas) 23 & ,ts_prof_cas,ps_prof_cas,ps_forc_prof_cas,plev_prof_cas,pforc_prof_cas& 24 & ,t_prof_cas,th_prof_cas,thl_prof_cas & 25 & ,qt_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas & 26 & ,rt_prof_cas,rv_prof_cas,rl_prof_cas,ri_prof_cas,rh_prof_cas & 27 & ,u_prof_cas,v_prof_cas,w_prof_cas,omega_prof_cas & 28 & ,ug_prof_cas,vg_prof_cas & 29 & ,temp_nudg_prof_cas,th_nudg_prof_cas,thl_nudg_prof_cas & 30 & ,qt_nudg_prof_cas,qv_nudg_prof_cas & 31 & ,rt_nudg_prof_cas,rv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 32 & ,uadv_prof_cas,vadv_prof_cas,tadv_prof_cas,thadv_prof_cas,thladv_prof_cas& 33 & ,qtadv_prof_cas,qvadv_prof_cas,rtadv_prof_cas,rvadv_prof_cas & 34 & ,trad_prof_cas,thrad_prof_cas,thlrad_prof_cas & 35 & ,tke_prof_cas,lat_prof_cas,sens_prof_cas,ustar_prof_cas & 36 & ,wpthetap_prof_cas,wpqtp_prof_cas,wpqvp_prof_cas,wprtp_prof_cas,wprvp_prof_cas) 30 37 31 38 ts_cur = ts_prof_cas … … 34 41 35 42 ! vertical interpolation: 36 CALL interp2_case_vertical_std(play,nlev_cas,plev_prof_cas & 37 & ,t_prof_cas,theta_prof_cas,thv_prof_cas,thl_prof_cas & 38 & ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas & 39 & ,ug_prof_cas,vg_prof_cas & 40 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 41 & ,vitw_prof_cas,omega_prof_cas & 42 & ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas & 43 & ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & 44 & ,dth_prof_cas,hth_prof_cas,vth_prof_cas & 43 CALL interp_case_vertical_std(nlev_cas & 44 & ,plev_prof_cas,t_prof_cas,th_prof_cas,thl_prof_cas & 45 & ,qt_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas & 46 & ,rt_prof_cas,rv_prof_cas,rl_prof_cas,ri_prof_cas,rh_prof_cas & 47 & ,u_prof_cas,v_prof_cas,w_prof_cas,omega_prof_cas & 48 & ,ug_prof_cas,vg_prof_cas & 49 & ,temp_nudg_prof_cas,th_nudg_prof_cas,thl_nudg_prof_cas & 50 & ,qt_nudg_prof_cas,qv_nudg_prof_cas & 51 & ,rt_nudg_prof_cas,rv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 52 & ,uadv_prof_cas,vadv_prof_cas,tadv_prof_cas,thadv_prof_cas,thladv_prof_cas & 53 & ,qtadv_prof_cas,qvadv_prof_cas,rtadv_prof_cas,rvadv_prof_cas & 54 & ,trad_prof_cas,thrad_prof_cas,thlrad_prof_cas & 45 55 ! 46 & ,t_mod_cas,theta_mod_cas,thv_mod_cas,thl_mod_cas,qv_mod_cas,ql_mod_cas,qi_mod_cas & 47 & ,u_mod_cas,v_mod_cas,ug_mod_cas,vg_mod_cas & 48 & ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas & 49 & ,w_mod_cas,omega_mod_cas & 50 & ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas & 51 & ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas & 52 & ,dth_mod_cas,hth_mod_cas,vth_mod_cas,mxcalc) 56 & ,plev_mod_cas,t_mod_cas,theta_mod_cas,thl_mod_cas & 57 & ,qt_mod_cas,qv_mod_cas,ql_mod_cas,qi_mod_cas & 58 & ,rt_mod_cas,rv_mod_cas,rl_mod_cas,ri_mod_cas,rh_mod_cas & 59 & ,u_mod_cas,v_mod_cas,w_mod_cas,omega_mod_cas & 60 & ,ug_mod_cas,vg_mod_cas & 61 & ,temp_nudg_mod_cas,th_nudg_mod_cas,thl_nudg_mod_cas & 62 & ,qt_nudg_mod_cas,qv_nudg_mod_cas & 63 & ,rt_nudg_mod_cas,rv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas & 64 & ,uadv_mod_cas,vadv_mod_cas,tadv_mod_cas,thadv_mod_cas,thladv_mod_cas & 65 & ,qtadv_mod_cas,qvadv_mod_cas,rtadv_mod_cas,rvadv_mod_cas & 66 & ,trad_mod_cas,thrad_mod_cas,thlrad_mod_cas) 53 67 54 68 … … 79 93 else 80 94 DO l=2,llm-1 81 95 IF (omega(l)>0.) THEN 82 96 d_t_z(l)=(temp(l+1)-temp(l))/(play(l+1)-play(l)) 83 97 d_th_z(l)=(teta(l+1)-teta(l))/(play(l+1)-play(l)) … … 85 99 d_u_z(l)=(u(l+1)-u(l))/(play(l+1)-play(l)) 86 100 d_v_z(l)=(v(l+1)-v(l))/(play(l+1)-play(l)) 87 101 ELSE 88 102 d_t_z(l)=(temp(l-1)-temp(l))/(play(l-1)-play(l)) 89 103 d_th_z(l)=(teta(l-1)-teta(l))/(play(l-1)-play(l)) … … 91 105 d_u_z(l)=(u(l-1)-u(l))/(play(l-1)-play(l)) 92 106 d_v_z(l)=(v(l-1)-v(l))/(play(l-1)-play(l)) 93 94 107 ENDIF 108 ENDDO 95 109 endif 96 110 d_t_z(1)=d_t_z(2) … … 157 171 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! 158 172 159 if ( trad.eq.1) then173 if (rad_temp.eq.1) then 160 174 tend_rayo=1 161 175 dt_cooling(l) = dtrad_mod_cas(l)
Note: See TracChangeset
for help on using the changeset viewer.
