- Timestamp:
- Jul 24, 2024, 4:23:34 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/dyn1d/old_1D_interp_cases.h
r5116 r5117 5 5 ! Forcing_LES case: constant dq_dyn 6 6 !--------------------------------------------------------------------- 7 if(forcing_LES) THEN7 IF (forcing_LES) THEN 8 8 DO l = 1,llm 9 9 d_q_adv(l,1) = dq_dyn(l,1) … … 14 14 ! Interpolation forcing in time and onto model levels 15 15 !--------------------------------------------------------------------- 16 if(forcing_GCSSold) THEN16 IF (forcing_GCSSold) THEN 17 17 CALL get_uvd(it,timestep,fich_gcssold_ctl,fich_gcssold_dat, & 18 18 & ht_gcssold,hq_gcssold,hw_gcssold, & … … 21 21 & imp_fcg_gcssold,ts_fcg_gcssold, & 22 22 & Tp_fcg_gcssold,Turb_fcg_gcssold) 23 if(prt_level.ge.1) THEN23 IF (prt_level.ge.1) THEN 24 24 print *,' get_uvd -> hqturb_gcssold ',it,hqturb_gcssold 25 25 endif … … 46 46 ! Interpolation Toga forcing 47 47 !--------------------------------------------------------------------- 48 if(forcing_toga) THEN49 if(prt_level.ge.1) THEN48 IF (forcing_toga) THEN 49 IF (prt_level.ge.1) THEN 50 50 PRINT*, & 51 51 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_toga=', & … … 61 61 & ,ht_prof,vt_prof,hq_prof,vq_prof) 62 62 ! EV: tg instead of ts_cur 63 if (type_ts_forcing.eq.1) tg = ts_prof !63 IF (type_ts_forcing.EQ.1) tg = ts_prof ! 64 64 65 65 ! vertical interpolation: … … 91 91 ! Interpolation DICE forcing 92 92 !--------------------------------------------------------------------- 93 if(forcing_dice) THEN94 if(prt_level.ge.1) THEN93 IF (forcing_dice) THEN 94 IF (prt_level.ge.1) THEN 95 95 PRINT*,'#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_dice=',& 96 96 & day,day1,(day-day1)*86400.,(day-day1)*86400/dt_dice … … 111 111 ! enddo 112 112 ! EV tg instead of ts_cur 113 if (type_ts_forcing.eq.1) tg = tg_prof ! SST used113 IF (type_ts_forcing.EQ.1) tg = tg_prof ! SST used 114 114 115 115 ! vertical interpolation: … … 191 191 ! Interpolation gabls4 forcing 192 192 !--------------------------------------------------------------------- 193 if(forcing_gabls4 ) THEN194 if(prt_level.ge.1) THEN193 IF (forcing_gabls4 ) THEN 194 IF (prt_level.ge.1) THEN 195 195 PRINT*,'#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_gabls4=',& 196 196 & day,day1,(day-day1)*86400.,(day-day1)*86400/dt_gabls4 … … 203 203 & ,ug_profg,vg_profg,ht_profg,hq_profg,tg_profg) 204 204 !EV tg instead of ts_cur 205 if (type_ts_forcing.eq.1) tg = tg_prof ! SST used205 IF (type_ts_forcing.EQ.1) tg = tg_prof ! SST used 206 206 207 207 ! vertical interpolation: … … 232 232 ! Interpolation forcing TWPice 233 233 !--------------------------------------------------------------------- 234 if(forcing_twpice) THEN234 IF (forcing_twpice) THEN 235 235 PRINT*, & 236 236 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_twpi=', & … … 277 277 !wind nudging above 500m with a 2h time scale 278 278 do l=1,llm 279 if(nudge_wind) THEN279 IF (nudge_wind) THEN 280 280 ! if (phi(l).gt.5000.) THEN 281 if(phi(l).gt.0.) THEN281 IF (phi(l).gt.0.) THEN 282 282 u(l)=u(l)+timestep*(u_mod(l)-u(l))/(2.*3600.) 283 283 v(l)=v(l)+timestep*(v_mod(l)-v(l))/(2.*3600.) … … 290 290 291 291 !CR:nudging of q and theta with a 6h time scale above 15km 292 if(nudge_thermo) THEN292 IF (nudge_thermo) THEN 293 293 do l=1,llm 294 294 zz(l)=phi(l)/9.8 295 if ((zz(l).le.16000.).and.(zz(l).gt.15000.)) THEN295 IF ((zz(l).le.16000.).AND.(zz(l).gt.15000.)) THEN 296 296 zfact=(zz(l)-15000.)/1000. 297 297 q(l,1)=q(l,1)+timestep*(q_mod(l)-q(l,1))/(6.*3600.)*zfact 298 298 temp(l)=temp(l)+timestep*(t_mod(l)-temp(l))/(6.*3600.)*zfact 299 else if(zz(l).gt.16000.) THEN299 ELSE IF (zz(l).gt.16000.) THEN 300 300 q(l,1)=q(l,1)+timestep*(q_mod(l)-q(l,1))/(6.*3600.) 301 301 temp(l)=temp(l)+timestep*(t_mod(l)-temp(l))/(6.*3600.) … … 309 309 alpha = rd*temp(l)*(1.+(rv/rd-1.)*q(l,1))/play(l) 310 310 !calcul de l'advection totale 311 if(cptadvw) THEN311 IF (cptadvw) THEN 312 312 d_t_adv(l) = alpha*omega(l)/rcpd+ht_mod(l)-d_t_dyn_z(l) 313 313 ! PRINT*,'temp vert adv',l,ht_mod(l),vt_mod(l),-d_t_dyn_z(l) … … 328 328 !--------------------------------------------------------------------- 329 329 330 if(forcing_amma) THEN330 IF (forcing_amma) THEN 331 331 PRINT*, & 332 332 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_amma=', & … … 410 410 ! Interpolation forcing Rico 411 411 !--------------------------------------------------------------------- 412 if(forcing_rico) THEN412 IF (forcing_rico) THEN 413 413 ! CALL lstendH(llm,omega,dt_dyn,dq_dyn,du_dyn, dv_dyn,q,temp,u,v,play) 414 414 CALL lstendH(llm,nqtot,omega,dt_dyn,dq_dyn,q,temp,u,v,play) … … 424 424 ! Interpolation forcing Arm_cu 425 425 !--------------------------------------------------------------------- 426 if(forcing_armcu) THEN426 IF (forcing_armcu) THEN 427 427 PRINT*, & 428 428 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_armcu=', & … … 479 479 ! Interpolation forcing in time and onto model levels 480 480 !--------------------------------------------------------------------- 481 if(forcing_sandu) THEN481 IF (forcing_sandu) THEN 482 482 PRINT*, & 483 483 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_sandu=', & … … 492 492 & ,ts_sandu,ts_prof) 493 493 ! EV tg instead of ts_cur 494 if (type_ts_forcing.eq.1) tg = ts_prof ! SST used in read_tsurf1d494 IF (type_ts_forcing.EQ.1) tg = ts_prof ! SST used in read_tsurf1d 495 495 496 496 ! vertical interpolation: … … 560 560 ! Interpolation forcing in time and onto model levels 561 561 !--------------------------------------------------------------------- 562 if(forcing_astex) THEN562 IF (forcing_astex) THEN 563 563 PRINT*, & 564 564 & '#### ITAP,day,day1,(day-day1)*86400,(day-day1)*86400/dt_astex=', & … … 574 574 & ,ufa_prof,vfa_prof) 575 575 ! EV tg instead of ts_cur 576 if (type_ts_forcing.eq.1) tg = ts_prof ! SST used576 IF (type_ts_forcing.EQ.1) tg = ts_prof ! SST used 577 577 ! vertical interpolation: 578 578 CALL interp_astex_vertical(play,nlev_astex,plev_profa & … … 642 642 ! Interpolation forcing standard case 643 643 !--------------------------------------------------------------------- 644 if(forcing_case) THEN644 IF (forcing_case) THEN 645 645 PRINT*,'FORCING CASE forcing_case' 646 646 … … 714 714 715 715 !wind nudging 716 if(nudge_u.gt.0.) THEN716 IF (nudge_u.gt.0.) THEN 717 717 do l=1,llm 718 718 u(l)=u(l)+timestep*(u_mod_cas(l)-u(l))/(nudge_u) … … 724 724 endif 725 725 726 if(nudge_v.gt.0.) THEN726 IF (nudge_v.gt.0.) THEN 727 727 do l=1,llm 728 728 v(l)=v(l)+timestep*(v_mod_cas(l)-v(l))/(nudge_v) … … 734 734 endif 735 735 736 if(nudge_w.gt.0.) THEN736 IF (nudge_w.gt.0.) THEN 737 737 do l=1,llm 738 738 w(l)=w(l)+timestep*(w_mod_cas(l)-w(l))/(nudge_w) … … 745 745 746 746 !nudging of q and temp 747 if(nudge_t.gt.0.) THEN747 IF (nudge_t.gt.0.) THEN 748 748 do l=1,llm 749 749 temp(l)=temp(l)+timestep*(t_mod_cas(l)-temp(l))/(nudge_t) 750 750 enddo 751 751 endif 752 if(nudge_q.gt.0.) THEN752 IF (nudge_q.gt.0.) THEN 753 753 do l=1,llm 754 754 q(l,1)=q(l,1)+timestep*(q_mod_cas(l)-q(l,1))/(nudge_q) … … 762 762 763 763 !calcul advection 764 if ((tend_u.eq.1).and.(tend_w.eq.0)) THEN764 IF ((tend_u.EQ.1).AND.(tend_w.EQ.0)) THEN 765 765 d_u_adv(l)=du_mod_cas(l) 766 else if ((tend_u.eq.1).and.(tend_w.eq.1)) THEN766 ELSE IF ((tend_u.EQ.1).AND.(tend_w.EQ.1)) THEN 767 767 d_u_adv(l)=hu_mod_cas(l)-d_u_dyn_z(l) 768 768 endif 769 769 770 if ((tend_v.eq.1).and.(tend_w.eq.0)) THEN770 IF ((tend_v.EQ.1).AND.(tend_w.EQ.0)) THEN 771 771 d_v_adv(l)=dv_mod_cas(l) 772 else if ((tend_v.eq.1).and.(tend_w.eq.1)) THEN772 ELSE IF ((tend_v.EQ.1).AND.(tend_w.EQ.1)) THEN 773 773 d_v_adv(l)=hv_mod_cas(l)-d_v_dyn_z(l) 774 774 endif 775 775 776 if ((tend_t.eq.1).and.(tend_w.eq.0)) THEN776 IF ((tend_t.EQ.1).AND.(tend_w.EQ.0)) THEN 777 777 ! d_t_adv(l)=alpha*omega(l)/rcpd+dt_mod_cas(l) 778 778 d_t_adv(l)=alpha*omega(l)/rcpd-dt_mod_cas(l) 779 else if ((tend_t.eq.1).and.(tend_w.eq.1)) THEN779 ELSE IF ((tend_t.EQ.1).AND.(tend_w.EQ.1)) THEN 780 780 ! d_t_adv(l)=alpha*omega(l)/rcpd+ht_mod_cas(l)-d_t_dyn_z(l) 781 781 d_t_adv(l)=alpha*omega(l)/rcpd-ht_mod_cas(l)-d_t_dyn_z(l) 782 782 endif 783 783 784 if ((tend_q.eq.1).and.(tend_w.eq.0)) THEN784 IF ((tend_q.EQ.1).AND.(tend_w.EQ.0)) THEN 785 785 ! d_q_adv(l,1)=dq_mod_cas(l) 786 786 d_q_adv(l,1)=-1*dq_mod_cas(l) 787 else if ((tend_q.eq.1).and.(tend_w.eq.1)) THEN787 ELSE IF ((tend_q.EQ.1).AND.(tend_w.EQ.1)) THEN 788 788 ! d_q_adv(l,1)=hq_mod_cas(l)-d_q_dyn_z(l) 789 789 d_q_adv(l,1)=-1*hq_mod_cas(l)-d_q_dyn_z(l) 790 790 endif 791 791 792 if (tend_rayo.eq.1) THEN792 IF (tend_rayo.EQ.1) THEN 793 793 dt_cooling(l) = dtrad_mod_cas(l) 794 794 ! print *,'dt_cooling=',dt_cooling(l) … … 814 814 ! Interpolation forcing standard case 815 815 !--------------------------------------------------------------------- 816 if(forcing_case2 .OR. forcing_SCM) THEN816 IF (forcing_case2 .OR. forcing_SCM) THEN 817 817 PRINT*,'FORCING CASE forcing_case2' 818 818 PRINT*, & … … 904 904 905 905 !geostrophic wind 906 if (forc_geo.eq.1) THEN906 IF (forc_geo.EQ.1) THEN 907 907 do l=1,llm 908 908 ug(l) = ug_mod_cas(l) … … 911 911 endif 912 912 !wind nudging 913 if(nudging_u.gt.0.) THEN913 IF (nudging_u.gt.0.) THEN 914 914 do l=1,llm 915 915 u(l)=u(l)+timestep*(u_mod_cas(l)-u(l))/(nudge_u) … … 921 921 endif 922 922 923 if(nudging_v.gt.0.) THEN923 IF (nudging_v.gt.0.) THEN 924 924 do l=1,llm 925 925 v(l)=v(l)+timestep*(v_mod_cas(l)-v(l))/(nudge_v) … … 931 931 endif 932 932 933 if(nudging_w.gt.0.) THEN933 IF (nudging_w.gt.0.) THEN 934 934 do l=1,llm 935 935 w(l)=w(l)+timestep*(w_mod_cas(l)-w(l))/(nudge_w) … … 942 942 943 943 !nudging of q and temp 944 if(nudging_t.gt.0.) THEN944 IF (nudging_t.gt.0.) THEN 945 945 do l=1,llm 946 946 temp(l)=temp(l)+timestep*(t_mod_cas(l)-temp(l))/(nudge_t) 947 947 enddo 948 948 endif 949 if(nudging_qv.gt.0.) THEN949 IF (nudging_qv.gt.0.) THEN 950 950 do l=1,llm 951 951 q(l,1)=q(l,1)+timestep*(q_mod_cas(l)-q(l,1))/(nudge_q) … … 960 960 961 961 !calcul advections 962 if ((forc_u.eq.1).and.(forc_w.eq.0)) THEN962 IF ((forc_u.EQ.1).AND.(forc_w.EQ.0)) THEN 963 963 d_u_adv(l)=du_mod_cas(l) 964 else if ((forc_u.eq.1).and.(forc_w.eq.1)) THEN964 ELSE IF ((forc_u.EQ.1).AND.(forc_w.EQ.1)) THEN 965 965 d_u_adv(l)=hu_mod_cas(l)-d_u_dyn_z(l) 966 966 endif 967 967 968 if ((forc_v.eq.1).and.(forc_w.eq.0)) THEN968 IF ((forc_v.EQ.1).AND.(forc_w.EQ.0)) THEN 969 969 d_v_adv(l)=dv_mod_cas(l) 970 else if ((forc_v.eq.1).and.(forc_w.eq.1)) THEN970 ELSE IF ((forc_v.EQ.1).AND.(forc_w.EQ.1)) THEN 971 971 d_v_adv(l)=hv_mod_cas(l)-d_v_dyn_z(l) 972 972 endif … … 974 974 ! Puisque dth a ete converti en dt, on traite de la meme facon 975 975 ! les flags tadv et thadv 976 if ((tadv.eq.1.or.thadv.eq.1) .and. (forc_w.eq.0)) THEN976 IF ((tadv.EQ.1.OR.thadv.EQ.1) .AND. (forc_w.EQ.0)) THEN 977 977 ! d_t_adv(l)=alpha*omega(l)/rcpd-dt_mod_cas(l) 978 978 d_t_adv(l)=alpha*omega(l)/rcpd+dt_mod_cas(l) 979 else if ((tadv.eq.1.or.thadv.eq.1) .and. (forc_w.eq.1)) THEN979 ELSE IF ((tadv.EQ.1.OR.thadv.EQ.1) .AND. (forc_w.EQ.1)) THEN 980 980 ! d_t_adv(l)=alpha*omega(l)/rcpd-ht_mod_cas(l)-d_t_dyn_z(l) 981 981 d_t_adv(l)=alpha*omega(l)/rcpd+ht_mod_cas(l)-d_t_dyn_z(l) 982 982 endif 983 983 984 ! if ((thadv. eq.1) .and. (forc_w.eq.0)) THEN984 ! if ((thadv.EQ.1) .AND. (forc_w.EQ.0)) THEN 985 985 ! d_t_adv(l)=alpha*omega(l)/rcpd-dth_mod_cas(l) 986 986 ! d_t_adv(l)=alpha*omega(l)/rcpd+dth_mod_cas(l) 987 ! else if ((thadv.eq.1) .and. (forc_w.eq.1)) THEN987 ! ELSE IF ((thadv.EQ.1) .AND. (forc_w.EQ.1)) THEN 988 988 ! d_t_adv(l)=alpha*omega(l)/rcpd-hth_mod_cas(l)-d_t_dyn_z(l) 989 989 ! d_t_adv(l)=alpha*omega(l)/rcpd+hth_mod_cas(l)-d_t_dyn_z(l) 990 990 ! endif 991 991 992 if ((qadv.eq.1) .and. (forc_w.eq.0)) THEN992 IF ((qadv.EQ.1) .AND. (forc_w.EQ.0)) THEN 993 993 d_q_adv(l,1)=dq_mod_cas(l) 994 994 ! d_q_adv(l,1)=-1*dq_mod_cas(l) 995 else if ((qadv.eq.1) .and. (forc_w.eq.1)) THEN995 ELSE IF ((qadv.EQ.1) .AND. (forc_w.EQ.1)) THEN 996 996 d_q_adv(l,1)=hq_mod_cas(l)-d_q_dyn_z(l) 997 997 ! d_q_adv(l,1)=-1*hq_mod_cas(l)-d_q_dyn_z(l) 998 998 endif 999 999 1000 if (trad.eq.1) THEN1000 IF (trad.EQ.1) THEN 1001 1001 tend_rayo=1 1002 1002 dt_cooling(l) = dtrad_mod_cas(l)
Note: See TracChangeset
for help on using the changeset viewer.
