- Timestamp:
- Jan 11, 2021, 11:24:08 PM (4 years ago)
- Location:
- LMDZ6/branches/Ocean_skin
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Ocean_skin
-
LMDZ6/branches/Ocean_skin/libf/phylmd/dyn1d/1D_interp_cases.h
r3605 r3798 1 1 2 2 print*,'FORCING CASE forcing_case2' 3 3 ! print*, & 4 4 ! & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/pdt_cas=', & … … 13 13 & ,u_cas,v_cas,ug_cas,vg_cas & 14 14 & ,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas & 15 & ,vitw_cas,omega_cas, du_cas,hu_cas,vu_cas &15 & ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas & 16 16 & ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas & 17 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 &18 & ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas & 19 19 ! 20 20 & ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas,thv_prof_cas & … … 22 22 & ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas & 23 23 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 24 & ,vitw_prof_cas,omega_prof_cas 24 & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & 25 25 & ,du_prof_cas,hu_prof_cas,vu_prof_cas & 26 26 & ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas & 27 27 & ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & 28 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)30 31 t s_cur= ts_prof_cas29 & ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas) 30 ! EV tg instead of ts_cur 31 tg = ts_prof_cas 32 32 ! psurf=plev_prof_cas(1) 33 33 psurf=ps_prof_cas 34 34 35 35 ! vertical interpolation: 36 CALL interp2_case_vertical_std(play, nlev_cas,plev_prof_cas &36 CALL interp2_case_vertical_std(play,plev,nlev_cas,plev_prof_cas & 37 37 & ,t_prof_cas,theta_prof_cas,thv_prof_cas,thl_prof_cas & 38 38 & ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas & 39 39 & ,ug_prof_cas,vg_prof_cas & 40 40 & ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas & 41 & ,vitw_prof_cas,omega_prof_cas &41 & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & 42 42 & ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas & 43 43 & ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas & … … 47 47 & ,u_mod_cas,v_mod_cas,ug_mod_cas,vg_mod_cas & 48 48 & ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas & 49 & ,w_mod_cas,omega_mod_cas 49 & ,w_mod_cas,omega_mod_cas,tke_mod_cas & 50 50 & ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas & 51 51 & ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas & … … 109 109 do l = 1, llm 110 110 ! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170309) 111 print*, l, llm 112 print*, play(l), temp(l) 111 113 omega(l) = -w_mod_cas(l)*play(l)*rg/(rd*temp(l)) 112 114 enddo
Note: See TracChangeset
for help on using the changeset viewer.