Changeset 4271


Ignore:
Timestamp:
Sep 21, 2022, 3:12:14 PM (20 months ago)
Author:
lguez
Message:

Remove useless & at beginning of continued line

Location:
LMDZ6/trunk/libf/phylmd/dyn1d
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read2.F90

    r4104 r4271  
    179179        print*,'Allocations OK'
    180180        call read_cas2(nid,nlev_cas,nt_cas                                       &
    181      &     ,z_cas,plev_cas,t_cas,q_cas,rh_cas,th_cas,rv_cas,u_cas,v_cas         &
    182      &     ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas,dv_cas,hv_cas,vv_cas    &
    183      &     ,dt_cas,dtrad_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas                 &
    184      &     ,dth_cas,hth_cas,vth_cas,dr_cas,hr_cas,vr_cas,sens_cas,lat_cas,ts_cas&
    185      &     ,ustar_cas,uw_cas,vw_cas,q1_cas,q2_cas)
     181           ,z_cas,plev_cas,t_cas,q_cas,rh_cas,th_cas,rv_cas,u_cas,v_cas         &
     182           ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas,dv_cas,hv_cas,vv_cas    &
     183           ,dt_cas,dtrad_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas                 &
     184           ,dth_cas,hth_cas,vth_cas,dr_cas,hr_cas,vr_cas,sens_cas,lat_cas,ts_cas&
     185           ,ustar_cas,uw_cas,vw_cas,q1_cas,q2_cas)
    186186        print*,'Read cas OK'
    187187
    … …  
    301301        print*,'Allocations OK'
    302302        call read2_cas (nid,nlev_cas,nt_cas,                                                                     &
    303      &     ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                    &
    304      &     ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,ug_cas,vg_cas,du_cas,hu_cas,vu_cas,        &
    305      &     dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
    306      &     dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tke_cas,                      &
    307      &     uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
    308      &     o3_cas,rugos_cas,clay_cas,sand_cas)
     303           ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                    &
     304           ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,ug_cas,vg_cas,du_cas,hu_cas,vu_cas,        &
     305           dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
     306           dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tke_cas,                      &
     307           uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
     308           o3_cas,rugos_cas,clay_cas,sand_cas)
    309309        print*,'Read2 cas OK'
    310310        do ii=1,nlev_cas
    … …  
    457457        print*,'Allocations OK'
    458458        call old_read_SCM (nid,nlev_cas,nt_cas,                                                                     &
    459      &     ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                    &
    460      &     ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,ug_cas,vg_cas,du_cas,hu_cas,vu_cas,        &
    461      &     dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
    462      &     dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tke_cas,                      &
    463      &     uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
    464      &     o3_cas,rugos_cas,clay_cas,sand_cas)
     459           ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                    &
     460           ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,ug_cas,vg_cas,du_cas,hu_cas,vu_cas,        &
     461           dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
     462           dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tke_cas,                      &
     463           uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
     464           o3_cas,rugos_cas,clay_cas,sand_cas)
    465465        print*,'Read2 cas OK'
    466466        do ii=1,nlev_cas
    … …  
    539539!=====================================================================
    540540      subroutine read_cas2(nid,nlevel,ntime                          &
    541      &     ,zz,pp,temp,qv,rh,theta,rv,u,v,ug,vg,w,                   &
    542      &     du,hu,vu,dv,hv,vv,dt,dtrad,ht,vt,dq,hq,vq,                     &
    543      &     dth,hth,vth,dr,hr,vr,sens,flat,ts,ustar,uw,vw,q1,q2)
     541           ,zz,pp,temp,qv,rh,theta,rv,u,v,ug,vg,w,                   &
     542           du,hu,vu,dv,hv,vv,dt,dtrad,ht,vt,dq,hq,vq,                     &
     543           dth,hth,vth,dr,hr,vr,sens,flat,ts,ustar,uw,vw,q1,q2)
    544544
    545545!program reading forcing of the case study
    … …  
    575575      character*5 name_var(1:nbvar3d)
    576576      data name_var/'zz','pp','temp','qv','rh','theta','rv','u','v','ug','vg','w','advu','hu','vu',&
    577      &'advv','hv','vv','advT','hT','vT','advq','hq','vq','advth','hth','vth','advr','hr','vr',&
    578      &'radT','uw','vw','q1','q2','sens','flat','ts','ustar'/
     577      'advv','hv','vv','advT','hT','vT','advq','hq','vq','advth','hth','vth','advr','hr','vr',&
     578      'radT','uw','vw','q1','q2','sens','flat','ts','ustar'/
    579579
    580580       do i=1,nbvar3d
    … …  
    660660!======================================================================
    661661      subroutine read2_cas(nid,nlevel,ntime,                                       &
    662      &     ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,ug,vg,&
    663      &     du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
    664      &     dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tke,uw,vw,q1,q2,       &
    665      &     orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
    666      &     heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
     662           ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,ug,vg,&
     663           du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
     664           dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tke,uw,vw,q1,q2,       &
     665           orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
     666           heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
    667667
    668668!program reading forcing of the case study
    … …  
    700700      character*12 name_var(1:nbvar3d)
    701701      data name_var/'coor_par_a','coor_par_b','height_h','pressure_h',&
    702      &'w','omega','ug','vg','uadv','uadvh','uadvv','vadv','vadvh','vadvv','tadv','tadvh','tadvv',&
    703      &'qadv','qadvh','qadvv','thadv','thadvh','thadvv','thladvh','radv','radvh','radvv','radcool','q1','q2','ustress','vstress', &
     702      'w','omega','ug','vg','uadv','uadvh','uadvv','vadv','vadvh','vadvv','tadv','tadvh','tadvv',&
     703      'qadv','qadvh','qadvv','thadv','thadvh','thadvv','thladvh','radv','radvh','radvv','radcool','q1','q2','ustress','vstress', &
    704704     'rh',&
    705      &'height_f','pressure_f','temp','theta','thv','thl','qv','ql','qi','rv','u','v',&
    706      &'sfc_sens_flx','sfc_lat_flx','ts','ps','ustar','tke',&
    707      &'orog','albedo','emiss','t_skin','q_skin','mom_rough','heat_rough','o3','rugos','clay','sand'/
     705      'height_f','pressure_f','temp','theta','thv','thl','qv','ql','qi','rv','u','v',&
     706      'sfc_sens_flx','sfc_lat_flx','ts','ps','ustar','tke',&
     707      'orog','albedo','emiss','t_skin','q_skin','mom_rough','heat_rough','o3','rugos','clay','sand'/
    708708      do i=1,nbvar3d
    709709        missing_var(i)=0.
    … …  
    846846!======================================================================
    847847      subroutine old_read_SCM(nid,nlevel,ntime,                                       &
    848      &     ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,ug,vg,&
    849      &     du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
    850      &     dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tke,uw,vw,q1,q2,       &
    851      &     orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
    852      &     heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
     848           ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,ug,vg,&
     849           du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
     850           dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tke,uw,vw,q1,q2,       &
     851           orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
     852           heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
    853853
    854854!program reading forcing of the case study
    … …  
    889889      character*13 name_var(1:nbvar3d)
    890890      data name_var/'coor_par_a','coor_par_b','height_h','pressure_h',&
    891      &'temp','qv','ql','qi','u','v','tke','pressure',&
    892      &'w','omega','ug','vg','uadv','uadvh','uadvv','vadv','vadvh','vadvv','tadv','tadvh','tadvv',&
    893      &'qvadv','qvadvh','qvadvv','thadv','thadvh','thadvv','thladvh','radv','radvh','radvv','radcool','q1','q2','ustress','vstress', &
     891      'temp','qv','ql','qi','u','v','tke','pressure',&
     892      'w','omega','ug','vg','uadv','uadvh','uadvv','vadv','vadvh','vadvv','tadv','tadvh','tadvv',&
     893      'qvadv','qvadvh','qvadvv','thadv','thadvh','thadvv','thladvh','radv','radvh','radvv','radcool','q1','q2','ustress','vstress', &
    894894     'rh',&
    895      &'height_f','pressure_forc','tempt','theta','thv','thl','qvt','qlt','qit','rv','ut','vt','tket',&
    896      &'sfc_sens_flx','sfc_lat_flx','ts','ps','ustar',&
    897      &'orog','albedo','emiss','t_skin','q_skin','mom_rough','heat_rough','o3','rugos','clay','sand'/
     895      'height_f','pressure_forc','tempt','theta','thv','thl','qvt','qlt','qit','rv','ut','vt','tket',&
     896      'sfc_sens_flx','sfc_lat_flx','ts','ps','ustar',&
     897      'orog','albedo','emiss','t_skin','q_skin','mom_rough','heat_rough','o3','rugos','clay','sand'/
    898898      do i=1,nbvar3d
    899899        missing_var(i)=0.
    … …  
    10781078        SUBROUTINE interp_case_time2(day,day1,annee_ref                &
    10791079!    &         ,year_cas,day_cas,nt_cas,pdt_forc,nlev_cas      &
    1080      &         ,nt_cas,nlev_cas                                       &
    1081      &         ,ts_cas,ps_cas,plev_cas,t_cas,q_cas,u_cas,v_cas               &
    1082      &         ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas           &
    1083      &         ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas   &
    1084      &         ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas,ustar_cas       &
    1085      &         ,uw_cas,vw_cas,q1_cas,q2_cas                           &
    1086      &         ,ts_prof_cas,plev_prof_cas,t_prof_cas,q_prof_cas       &
    1087      &         ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas         &
    1088      &         ,vitw_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas     &
    1089      &         ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas       &
    1090      &         ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas    &
    1091      &         ,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas    &
    1092      &         ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas)
     1080               ,nt_cas,nlev_cas                                       &
     1081               ,ts_cas,ps_cas,plev_cas,t_cas,q_cas,u_cas,v_cas               &
     1082               ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas           &
     1083               ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas   &
     1084               ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas,ustar_cas       &
     1085               ,uw_cas,vw_cas,q1_cas,q2_cas                           &
     1086               ,ts_prof_cas,plev_prof_cas,t_prof_cas,q_prof_cas       &
     1087               ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas         &
     1088               ,vitw_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas     &
     1089               ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas       &
     1090               ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas    &
     1091               ,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas    &
     1092               ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas)
    10931093         
    10941094
    … …  
    12021202       if (it_cas1 .gt. nt_cas) then
    12031203        write(*,*) 'PB-stop: day, day_ju_ini_cas,it_cas1, it_cas2, timeit: '            &
    1204      &        ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
     1204              ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
    12051205        stop
    12061206       endif
    … …  
    12151215
    12161216       lat_prof_cas = lat_cas(it_cas2)                                       &
    1217      &          -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
     1217                -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
    12181218       sens_prof_cas = sens_cas(it_cas2)                                     &
    1219      &          -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
     1219                -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
    12201220       ts_prof_cas = ts_cas(it_cas2)                                         &
    1221      &          -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
     1221                -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
    12221222       ustar_prof_cas = ustar_cas(it_cas2)                                   &
    1223      &          -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
     1223                -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
    12241224
    12251225       do k=1,nlev_cas
    12261226        plev_prof_cas(k) = plev_cas(k,it_cas2)                               &
    1227      &          -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
     1227                -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
    12281228        t_prof_cas(k) = t_cas(k,it_cas2)                               &
    1229      &          -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
     1229                -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
    12301230        q_prof_cas(k) = q_cas(k,it_cas2)                               &
    1231      &          -frac*(q_cas(k,it_cas2)-q_cas(k,it_cas1))
     1231                -frac*(q_cas(k,it_cas2)-q_cas(k,it_cas1))
    12321232        u_prof_cas(k) = u_cas(k,it_cas2)                               &
    1233      &          -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
     1233                -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
    12341234        v_prof_cas(k) = v_cas(k,it_cas2)                               &
    1235      &          -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
     1235                -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
    12361236        ug_prof_cas(k) = ug_cas(k,it_cas2)                               &
    1237      &          -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
     1237                -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
    12381238        vg_prof_cas(k) = vg_cas(k,it_cas2)                               &
    1239      &          -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
     1239                -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
    12401240        vitw_prof_cas(k) = vitw_cas(k,it_cas2)                               &
    1241      &          -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
     1241                -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
    12421242        du_prof_cas(k) = du_cas(k,it_cas2)                                   &
    1243      &          -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
     1243                -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
    12441244        hu_prof_cas(k) = hu_cas(k,it_cas2)                                   &
    1245      &          -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
     1245                -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
    12461246        vu_prof_cas(k) = vu_cas(k,it_cas2)                                   &
    1247      &          -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
     1247                -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
    12481248        dv_prof_cas(k) = dv_cas(k,it_cas2)                                   &
    1249      &          -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
     1249                -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
    12501250        hv_prof_cas(k) = hv_cas(k,it_cas2)                                   &
    1251      &          -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
     1251                -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
    12521252        vv_prof_cas(k) = vv_cas(k,it_cas2)                                   &
    1253      &          -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
     1253                -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
    12541254        dt_prof_cas(k) = dt_cas(k,it_cas2)                                   &
    1255      &          -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
     1255                -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
    12561256        ht_prof_cas(k) = ht_cas(k,it_cas2)                                   &
    1257      &          -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
     1257                -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
    12581258        vt_prof_cas(k) = vt_cas(k,it_cas2)                                   &
    1259      &          -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
     1259                -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
    12601260        dtrad_prof_cas(k) = dtrad_cas(k,it_cas2)                                   &
    1261      &          -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
     1261                -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
    12621262        dq_prof_cas(k) = dq_cas(k,it_cas2)                                   &
    1263      &          -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
     1263                -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
    12641264        hq_prof_cas(k) = hq_cas(k,it_cas2)                                   &
    1265      &          -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
     1265                -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
    12661266        vq_prof_cas(k) = vq_cas(k,it_cas2)                                   &
    1267      &          -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
     1267                -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
    12681268       uw_prof_cas(k) = uw_cas(k,it_cas2)                                   &
    1269      &          -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
     1269                -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
    12701270       vw_prof_cas(k) = vw_cas(k,it_cas2)                                   &
    1271      &          -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
     1271                -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
    12721272       q1_prof_cas(k) = q1_cas(k,it_cas2)                                   &
    1273      &          -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
     1273                -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
    12741274       q2_prof_cas(k) = q2_cas(k,it_cas2)                                   &
    1275      &          -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
     1275                -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
    12761276        enddo
    12771277
    … …  
    12821282        SUBROUTINE interp2_case_time(day,day1,annee_ref                           &
    12831283!    &         ,year_cas,day_cas,nt_cas,pdt_forc,nlev_cas                         &
    1284      &         ,nt_cas,nlev_cas                                                   &
    1285      &         ,ts_cas,ps_cas,plev_cas,t_cas,theta_cas,thv_cas,thl_cas            &
    1286      &         ,qv_cas,ql_cas,qi_cas,u_cas,v_cas                                  &
    1287      &         ,ug_cas,vg_cas,vitw_cas,omega_cas,du_cas,hu_cas,vu_cas             &
    1288      &         ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas               &
    1289      &         ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas                      &
    1290      &         ,lat_cas,sens_cas,ustar_cas                                        &
    1291      &         ,uw_cas,vw_cas,q1_cas,q2_cas,tke_cas                               &
     1284               ,nt_cas,nlev_cas                                                   &
     1285               ,ts_cas,ps_cas,plev_cas,t_cas,theta_cas,thv_cas,thl_cas            &
     1286               ,qv_cas,ql_cas,qi_cas,u_cas,v_cas                                  &
     1287               ,ug_cas,vg_cas,vitw_cas,omega_cas,du_cas,hu_cas,vu_cas             &
     1288               ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas               &
     1289               ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas                      &
     1290               ,lat_cas,sens_cas,ustar_cas                                        &
     1291               ,uw_cas,vw_cas,q1_cas,q2_cas,tke_cas                               &
    12921292!
    1293      &         ,ts_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas               &
    1294      &         ,thv_prof_cas,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas     &
    1295      &         ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas                     &
    1296      &         ,vitw_prof_cas,omega_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas  &
    1297      &         ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas                   &
    1298      &         ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas                &
    1299      &         ,hq_prof_cas,vq_prof_cas,dth_prof_cas,hth_prof_cas,vth_prof_cas    &
    1300      &         ,lat_prof_cas,sens_prof_cas                                        &
    1301      &         ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas)
     1293               ,ts_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas               &
     1294               ,thv_prof_cas,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas     &
     1295               ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas                     &
     1296               ,vitw_prof_cas,omega_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas  &
     1297               ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas                   &
     1298               ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas                &
     1299               ,hq_prof_cas,vq_prof_cas,dth_prof_cas,hth_prof_cas,vth_prof_cas    &
     1300               ,lat_prof_cas,sens_prof_cas                                        &
     1301               ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tke_prof_cas)
    13021302         
    13031303
    … …  
    14181418       if (it_cas1 .gt. nt_cas) then
    14191419        write(*,*) 'PB-stop: day, day_ju_ini_cas,it_cas1, it_cas2, timeit: '            &
    1420      &        ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
     1420              ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
    14211421        stop
    14221422       endif
    … …  
    14311431
    14321432       lat_prof_cas = lat_cas(it_cas2)                                   &
    1433      &          -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
     1433                -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
    14341434       sens_prof_cas = sens_cas(it_cas2)                                 &
    1435      &          -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
     1435                -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
    14361436       tke_prof_cas = tke_cas(it_cas2)                                   &
    1437      &          -frac*(tke_cas(it_cas2)-tke_cas(it_cas1))
     1437                -frac*(tke_cas(it_cas2)-tke_cas(it_cas1))
    14381438       ts_prof_cas = ts_cas(it_cas2)                                     &
    1439      &          -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
     1439                -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
    14401440       ustar_prof_cas = ustar_cas(it_cas2)                               &
    1441      &          -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
     1441                -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
    14421442
    14431443       do k=1,nlev_cas
    14441444        plev_prof_cas(k) = plev_cas(k,it_cas2)                           &     
    1445      &          -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
     1445                -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
    14461446        t_prof_cas(k) = t_cas(k,it_cas2)                                 &       
    1447      &          -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
     1447                -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
    14481448        print *,'k,frac,plev_cas1,plev_cas2=',k,frac,plev_cas(k,it_cas1),plev_cas(k,it_cas2)
    14491449        theta_prof_cas(k) = theta_cas(k,it_cas2)                         &                     
    1450      &          -frac*(theta_cas(k,it_cas2)-theta_cas(k,it_cas1))
     1450                -frac*(theta_cas(k,it_cas2)-theta_cas(k,it_cas1))
    14511451        thv_prof_cas(k) = thv_cas(k,it_cas2)                             &         
    1452      &          -frac*(thv_cas(k,it_cas2)-thv_cas(k,it_cas1))
     1452                -frac*(thv_cas(k,it_cas2)-thv_cas(k,it_cas1))
    14531453        thl_prof_cas(k) = thl_cas(k,it_cas2)                             &             
    1454      &          -frac*(thl_cas(k,it_cas2)-thl_cas(k,it_cas1))
     1454                -frac*(thl_cas(k,it_cas2)-thl_cas(k,it_cas1))
    14551455        qv_prof_cas(k) = qv_cas(k,it_cas2)                               &
    1456      &          -frac*(qv_cas(k,it_cas2)-qv_cas(k,it_cas1))
     1456                -frac*(qv_cas(k,it_cas2)-qv_cas(k,it_cas1))
    14571457        ql_prof_cas(k) = ql_cas(k,it_cas2)                               &
    1458      &          -frac*(ql_cas(k,it_cas2)-ql_cas(k,it_cas1))
     1458                -frac*(ql_cas(k,it_cas2)-ql_cas(k,it_cas1))
    14591459        qi_prof_cas(k) = qi_cas(k,it_cas2)                               &
    1460      &          -frac*(qi_cas(k,it_cas2)-qi_cas(k,it_cas1))
     1460                -frac*(qi_cas(k,it_cas2)-qi_cas(k,it_cas1))
    14611461        u_prof_cas(k) = u_cas(k,it_cas2)                                 &
    1462      &          -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
     1462                -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
    14631463        v_prof_cas(k) = v_cas(k,it_cas2)                                 &
    1464      &          -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
     1464                -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
    14651465        ug_prof_cas(k) = ug_cas(k,it_cas2)                               &
    1466      &          -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
     1466                -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
    14671467        vg_prof_cas(k) = vg_cas(k,it_cas2)                               &
    1468      &          -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
     1468                -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
    14691469        vitw_prof_cas(k) = vitw_cas(k,it_cas2)                           &
    1470      &          -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
     1470                -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
    14711471        omega_prof_cas(k) = omega_cas(k,it_cas2)                         &
    1472      &          -frac*(omega_cas(k,it_cas2)-omega_cas(k,it_cas1))
     1472                -frac*(omega_cas(k,it_cas2)-omega_cas(k,it_cas1))
    14731473        du_prof_cas(k) = du_cas(k,it_cas2)                               &
    1474      &          -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
     1474                -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
    14751475        hu_prof_cas(k) = hu_cas(k,it_cas2)                               &
    1476      &          -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
     1476                -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
    14771477        vu_prof_cas(k) = vu_cas(k,it_cas2)                               &
    1478      &          -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
     1478                -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
    14791479        dv_prof_cas(k) = dv_cas(k,it_cas2)                               &
    1480      &          -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
     1480                -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
    14811481        hv_prof_cas(k) = hv_cas(k,it_cas2)                               &
    1482      &          -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
     1482                -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
    14831483        vv_prof_cas(k) = vv_cas(k,it_cas2)                               &
    1484      &          -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
     1484                -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
    14851485        dt_prof_cas(k) = dt_cas(k,it_cas2)                               &
    1486      &          -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
     1486                -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
    14871487        ht_prof_cas(k) = ht_cas(k,it_cas2)                               &
    1488      &          -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
     1488                -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
    14891489        vt_prof_cas(k) = vt_cas(k,it_cas2)                               &
    1490      &          -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
     1490                -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
    14911491        dth_prof_cas(k) = dth_cas(k,it_cas2)                             &
    1492      &          -frac*(dth_cas(k,it_cas2)-dth_cas(k,it_cas1))
     1492                -frac*(dth_cas(k,it_cas2)-dth_cas(k,it_cas1))
    14931493        hth_prof_cas(k) = hth_cas(k,it_cas2)                             &
    1494      &          -frac*(hth_cas(k,it_cas2)-hth_cas(k,it_cas1))
     1494                -frac*(hth_cas(k,it_cas2)-hth_cas(k,it_cas1))
    14951495        vth_prof_cas(k) = vth_cas(k,it_cas2)                             &
    1496      &          -frac*(vth_cas(k,it_cas2)-vth_cas(k,it_cas1))
     1496                -frac*(vth_cas(k,it_cas2)-vth_cas(k,it_cas1))
    14971497        dtrad_prof_cas(k) = dtrad_cas(k,it_cas2)                         &
    1498      &          -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
     1498                -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
    14991499        dq_prof_cas(k) = dq_cas(k,it_cas2)                               &
    1500      &          -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
     1500                -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
    15011501        hq_prof_cas(k) = hq_cas(k,it_cas2)                               &
    1502      &          -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
     1502                -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
    15031503        vq_prof_cas(k) = vq_cas(k,it_cas2)                               &
    1504      &          -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
     1504                -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
    15051505       uw_prof_cas(k) = uw_cas(k,it_cas2)                                &
    1506      &          -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
     1506                -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
    15071507       vw_prof_cas(k) = vw_cas(k,it_cas2)                                &
    1508      &          -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
     1508                -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
    15091509       q1_prof_cas(k) = q1_cas(k,it_cas2)                                &
    1510      &          -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
     1510                -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
    15111511       q2_prof_cas(k) = q2_cas(k,it_cas2)                                &
    1512      &          -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
     1512                -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
    15131513        enddo
    15141514

  • LMDZ6/trunk/libf/phylmd/dyn1d/mod_1D_cases_read_std.F90

    r4104 r4271  
    235235        print*,'Allocations OK'
    236236        CALL read_SCM (nid,nlev_cas,nt_cas,                                                                     &
    237      &     ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                   &
    238      &     ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,tke_cas,ug_cas,vg_cas,                            &
    239      &     temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas,                                                     &
    240      &     invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas,                         &
    241      &     du_cas,hu_cas,vu_cas,                                                                                &
    242      &     dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
    243      &     dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tkes_cas,                      &
    244      &     uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
    245      &     o3_cas,rugos_cas,clay_cas,sand_cas)
     237           ap_cas,bp_cas,z_cas,plev_cas,zh_cas,plevh_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,                   &
     238           ql_cas,qi_cas,rh_cas,rv_cas,u_cas,v_cas,vitw_cas,omega_cas,tke_cas,ug_cas,vg_cas,                            &
     239           temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas,                                                     &
     240           invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas,                         &
     241           du_cas,hu_cas,vu_cas,                                                                                &
     242           dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,               &
     243           dr_cas,hr_cas,vr_cas,dtrad_cas,sens_cas,lat_cas,ts_cas,ps_cas,ustar_cas,tkes_cas,                      &
     244           uw_cas,vw_cas,q1_cas,q2_cas,orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,heat_rough, &
     245           o3_cas,rugos_cas,clay_cas,sand_cas)
    246246        print*,'read_SCM cas OK'
    247247        do ii=1,nlev_cas
    … …  
    323323!=====================================================================
    324324      SUBROUTINE read_SCM(nid,nlevel,ntime,                                       &
    325      &     ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,tke,ug,vg,&
    326      &     temp_nudg,qv_nudg,u_nudg,v_nudg,                                        &
    327      &     invtau_temp_nudg,invtau_qv_nudg,invtau_u_nudg,invtau_v_nudg,             &
    328      &     du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
    329      &     dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tkes,uw,vw,q1,q2,       &
    330      &     orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
    331      &     heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
     325           ap,bp,zz,pp,zzh,pph,temp,theta,thv,thl,qv,ql,qi,rh,rv,u,v,vitw,omega,tke,ug,vg,&
     326           temp_nudg,qv_nudg,u_nudg,v_nudg,                                        &
     327           invtau_temp_nudg,invtau_qv_nudg,invtau_u_nudg,invtau_v_nudg,             &
     328           du,hu,vu,dv,hv,vv,dt,ht,vt,dq,hq,vq,                                    &
     329           dth,hth,vth,dr,hr,vr,dtrad,sens,flat,ts,ps,ustar,tkes,uw,vw,q1,q2,       &
     330           orog_cas,albedo_cas,emiss_cas,t_skin_cas,q_skin_cas,mom_rough,          &
     331           heat_rough,o3_cas,rugos_cas,clay_cas,sand_cas)
    332332
    333333!program reading forcing of the case study
    … …  
    393393      data name_var/ &
    394394     ! coordonnees pression (n+1 niveaux) #4
    395      & 'coor_par_a','coor_par_b','zf','pressure_h',& ! #1-#4
     395       'coor_par_a','coor_par_b','zf','pressure_h',& ! #1-#4
    396396     ! coordonnees pression (n niveaux) #8
    397      & 'ta','qv','ql','qi','ua','va','tke','pa',& ! #5-#12
     397       'ta','qv','ql','qi','ua','va','tke','pa',& ! #5-#12
    398398     ! coordonnees pression + temps #46
    399      & 'wa','wap','ug','vg','tnua_adv','tnua_advh','tnua_advv','tnva_adv','tnva_advh','tnva_advv','tnta_adv','tnta_advh','tnta_advv',& !  #13 - #25
    400      & 'tnqv_adv','tnqv_advh','tnqv_advv','thadv','thadvh','thadvv','thladvh',                             & ! #26 - #32
    401      & 'radv','radvh','radvv','tnta_rad','q1','q2','ustress','vstress',                           & ! #33 - #40
    402      & 'rh','ta_nud','qv_nud','ua_nud','va_nud',                                       & ! #41-45
    403      & 'zh_forc','pa_forc','tat','thetat','thetavt','thetalt','qvt','qlt','qit','rv','uat','vat',   & ! #46-57
    404      & 'nudging_constant_ta', 'nudging_constant_qv', 'nudging_constant_ua', 'nudging_constant_va',           & ! # 58-61
     399       'wa','wap','ug','vg','tnua_adv','tnua_advh','tnua_advv','tnva_adv','tnva_advh','tnva_advv','tnta_adv','tnta_advh','tnta_advv',& !  #13 - #25
     400       'tnqv_adv','tnqv_advh','tnqv_advv','thadv','thadvh','thadvv','thladvh',                             & ! #26 - #32
     401       'radv','radvh','radvv','tnta_rad','q1','q2','ustress','vstress',                           & ! #33 - #40
     402       'rh','ta_nud','qv_nud','ua_nud','va_nud',                                       & ! #41-45
     403       'zh_forc','pa_forc','tat','thetat','thetavt','thetalt','qvt','qlt','qit','rv','uat','vat',   & ! #46-57
     404       'nudging_constant_ta', 'nudging_constant_qv', 'nudging_constant_ua', 'nudging_constant_va',           & ! # 58-61
    405405     ! coordonnees temps #12
    406      & 'tkes','hfss','hfls','ts_forc','ps_forc','ustar', &                                        ! 62-67
    407      & 'orog','albedo','emiss','t_skin','q_skin','z0','z0h',       &                    ! 68-74
     406       'tkes','hfss','hfls','ts_forc','ps_forc','ustar', &                                        ! 62-67
     407       'orog','albedo','emiss','t_skin','q_skin','z0','z0h',       &                    ! 68-74
    408408     ! scalaires #4                               
    409      & 'O3','rugos','clay','sand'/                                                      ! 75-78
     409       'O3','rugos','clay','sand'/                                                      ! 75-78
    410410
    411411
    … …  
    654654        SUBROUTINE interp_case_time_std(day,day1,annee_ref                           &
    655655!    &         ,year_cas,day_cas,nt_cas,pdt_forc,nlev_cas                         &
    656      &         ,nt_cas,nlev_cas                                                   &
    657      &         ,ts_cas,ps_cas,plev_cas,t_cas,theta_cas,thv_cas,thl_cas            &
    658      &         ,qv_cas,ql_cas,qi_cas,u_cas,v_cas                                  &
    659      &         ,ug_cas,vg_cas,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas     &
    660      &         ,invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas     &
    661      &         ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas             &
    662      &         ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas               &
    663      &         ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas                      &
    664      &         ,lat_cas,sens_cas,ustar_cas                                        &
    665      &         ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas                               &
     656               ,nt_cas,nlev_cas                                                   &
     657               ,ts_cas,ps_cas,plev_cas,t_cas,theta_cas,thv_cas,thl_cas            &
     658               ,qv_cas,ql_cas,qi_cas,u_cas,v_cas                                  &
     659               ,ug_cas,vg_cas,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas     &
     660               ,invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas     &
     661               ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas             &
     662               ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas               &
     663               ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas                      &
     664               ,lat_cas,sens_cas,ustar_cas                                        &
     665               ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas                               &
    666666!
    667      &         ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas               &
    668      &         ,thv_prof_cas,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas     &
    669      &         ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas                     &
    670      &         ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas     &
    671      &         ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas     &     
    672      &         ,vitw_prof_cas,omega_prof_cas,tke_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas  &
    673      &         ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas                   &
    674      &         ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas                &
    675      &         ,hq_prof_cas,vq_prof_cas,dth_prof_cas,hth_prof_cas,vth_prof_cas    &
    676      &         ,lat_prof_cas,sens_prof_cas                                        &
    677      &         ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas)
     667               ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas               &
     668               ,thv_prof_cas,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas     &
     669               ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas                     &
     670               ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas     &
     671               ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas     &     
     672               ,vitw_prof_cas,omega_prof_cas,tke_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas  &
     673               ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas                   &
     674               ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas                &
     675               ,hq_prof_cas,vq_prof_cas,dth_prof_cas,hth_prof_cas,vth_prof_cas    &
     676               ,lat_prof_cas,sens_prof_cas                                        &
     677               ,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas)
    678678
    679679
    … …  
    811811       if (it_cas1 .gt. nt_cas) then
    812812        write(*,*) 'PB-stop: day, day_ju_ini_cas,it_cas1, it_cas2, timeit: '            &
    813      &        ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
     813              ,day,day_ju_ini_cas,it_cas1,it_cas2,timeit
    814814        stop
    815815       endif
    … …  
    824824
    825825       lat_prof_cas = lat_cas(it_cas2)                                   &
    826      &          -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
     826                -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
    827827       sens_prof_cas = sens_cas(it_cas2)                                 &
    828      &          -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
     828                -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
    829829       tkes_prof_cas = tkes_cas(it_cas2)                                   &
    830      &          -frac*(tkes_cas(it_cas2)-tkes_cas(it_cas1))
     830                -frac*(tkes_cas(it_cas2)-tkes_cas(it_cas1))
    831831       ts_prof_cas = ts_cas(it_cas2)                                     &
    832      &          -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
     832                -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
    833833       ps_prof_cas = ps_cas(it_cas2)                                     &
    834      &          -frac*(ps_cas(it_cas2)-ps_cas(it_cas1))
     834                -frac*(ps_cas(it_cas2)-ps_cas(it_cas1))
    835835       ustar_prof_cas = ustar_cas(it_cas2)                               &
    836      &          -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
     836                -frac*(ustar_cas(it_cas2)-ustar_cas(it_cas1))
    837837
    838838       do k=1,nlev_cas
    839839        plev_prof_cas(k) = plev_cas(k,it_cas2)                           &     
    840      &          -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
     840                -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
    841841        t_prof_cas(k) = t_cas(k,it_cas2)                                 &       
    842      &          -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
     842                -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
    843843        !print *,'k,frac,plev_cas1,plev_cas2=',k,frac,plev_cas(k,it_cas1),plev_cas(k,it_cas2)
    844844        theta_prof_cas(k) = theta_cas(k,it_cas2)                         &                     
    845      &          -frac*(theta_cas(k,it_cas2)-theta_cas(k,it_cas1))
     845                -frac*(theta_cas(k,it_cas2)-theta_cas(k,it_cas1))
    846846        thv_prof_cas(k) = thv_cas(k,it_cas2)                             &         
    847      &          -frac*(thv_cas(k,it_cas2)-thv_cas(k,it_cas1))
     847                -frac*(thv_cas(k,it_cas2)-thv_cas(k,it_cas1))
    848848        thl_prof_cas(k) = thl_cas(k,it_cas2)                             &             
    849      &          -frac*(thl_cas(k,it_cas2)-thl_cas(k,it_cas1))
     849                -frac*(thl_cas(k,it_cas2)-thl_cas(k,it_cas1))
    850850        qv_prof_cas(k) = qv_cas(k,it_cas2)                               &
    851      &          -frac*(qv_cas(k,it_cas2)-qv_cas(k,it_cas1))
     851                -frac*(qv_cas(k,it_cas2)-qv_cas(k,it_cas1))
    852852        ql_prof_cas(k) = ql_cas(k,it_cas2)                               &
    853      &          -frac*(ql_cas(k,it_cas2)-ql_cas(k,it_cas1))
     853                -frac*(ql_cas(k,it_cas2)-ql_cas(k,it_cas1))
    854854        qi_prof_cas(k) = qi_cas(k,it_cas2)                               &
    855      &          -frac*(qi_cas(k,it_cas2)-qi_cas(k,it_cas1))
     855                -frac*(qi_cas(k,it_cas2)-qi_cas(k,it_cas1))
    856856        u_prof_cas(k) = u_cas(k,it_cas2)                                 &
    857      &          -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
     857                -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
    858858        v_prof_cas(k) = v_cas(k,it_cas2)                                 &
    859      &          -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
     859                -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
    860860        ug_prof_cas(k) = ug_cas(k,it_cas2)                               &
    861      &          -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
     861                -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
    862862        vg_prof_cas(k) = vg_cas(k,it_cas2)                               &
    863      &          -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
     863                -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
    864864        temp_nudg_prof_cas(k) = temp_nudg_cas(k,it_cas2)                    &
    865      &          -frac*(temp_nudg_cas(k,it_cas2)-temp_nudg_cas(k,it_cas1))
     865                -frac*(temp_nudg_cas(k,it_cas2)-temp_nudg_cas(k,it_cas1))
    866866        qv_nudg_prof_cas(k) = qv_nudg_cas(k,it_cas2)                        &
    867      &          -frac*(qv_nudg_cas(k,it_cas2)-qv_nudg_cas(k,it_cas1))
     867                -frac*(qv_nudg_cas(k,it_cas2)-qv_nudg_cas(k,it_cas1))
    868868        u_nudg_prof_cas(k) = u_nudg_cas(k,it_cas2)                          &
    869      &          -frac*(u_nudg_cas(k,it_cas2)-u_nudg_cas(k,it_cas1))
     869                -frac*(u_nudg_cas(k,it_cas2)-u_nudg_cas(k,it_cas1))
    870870        v_nudg_prof_cas(k) = v_nudg_cas(k,it_cas2)                          &
    871      &          -frac*(v_nudg_cas(k,it_cas2)-v_nudg_cas(k,it_cas1))
     871                -frac*(v_nudg_cas(k,it_cas2)-v_nudg_cas(k,it_cas1))
    872872        invtau_temp_nudg_prof_cas(k) = invtau_temp_nudg_cas(k,it_cas2)                    &
    873      &          -frac*(invtau_temp_nudg_cas(k,it_cas2)-invtau_temp_nudg_cas(k,it_cas1))
     873                -frac*(invtau_temp_nudg_cas(k,it_cas2)-invtau_temp_nudg_cas(k,it_cas1))
    874874        invtau_qv_nudg_prof_cas(k) = invtau_qv_nudg_cas(k,it_cas2)                        &
    875      &          -frac*(invtau_qv_nudg_cas(k,it_cas2)-invtau_qv_nudg_cas(k,it_cas1))
     875                -frac*(invtau_qv_nudg_cas(k,it_cas2)-invtau_qv_nudg_cas(k,it_cas1))
    876876        invtau_u_nudg_prof_cas(k) = invtau_u_nudg_cas(k,it_cas2)                          &
    877      &          -frac*(invtau_u_nudg_cas(k,it_cas2)-invtau_u_nudg_cas(k,it_cas1))
     877                -frac*(invtau_u_nudg_cas(k,it_cas2)-invtau_u_nudg_cas(k,it_cas1))
    878878        invtau_v_nudg_prof_cas(k) = invtau_v_nudg_cas(k,it_cas2)                          &
    879      &          -frac*(invtau_v_nudg_cas(k,it_cas2)-invtau_v_nudg_cas(k,it_cas1))
     879                -frac*(invtau_v_nudg_cas(k,it_cas2)-invtau_v_nudg_cas(k,it_cas1))
    880880        vitw_prof_cas(k) = vitw_cas(k,it_cas2)                           &
    881      &          -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
     881                -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
    882882        omega_prof_cas(k) = omega_cas(k,it_cas2)                         &
    883      &          -frac*(omega_cas(k,it_cas2)-omega_cas(k,it_cas1))
     883                -frac*(omega_cas(k,it_cas2)-omega_cas(k,it_cas1))
    884884        tke_prof_cas(k) = tke_cas(k,it_cas2)                         &
    885      &          -frac*(tke_cas(k,it_cas2)-tke_cas(k,it_cas1))
     885                -frac*(tke_cas(k,it_cas2)-tke_cas(k,it_cas1))
    886886        du_prof_cas(k) = du_cas(k,it_cas2)                               &
    887      &          -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
     887                -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
    888888        hu_prof_cas(k) = hu_cas(k,it_cas2)                               &
    889      &          -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
     889                -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
    890890        vu_prof_cas(k) = vu_cas(k,it_cas2)                               &
    891      &          -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
     891                -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
    892892        dv_prof_cas(k) = dv_cas(k,it_cas2)                               &
    893      &          -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
     893                -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
    894894        hv_prof_cas(k) = hv_cas(k,it_cas2)                               &
    895      &          -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
     895                -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
    896896        vv_prof_cas(k) = vv_cas(k,it_cas2)                               &
    897      &          -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
     897                -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
    898898        dt_prof_cas(k) = dt_cas(k,it_cas2)                               &
    899      &          -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
     899                -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
    900900        ht_prof_cas(k) = ht_cas(k,it_cas2)                               &
    901      &          -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
     901                -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
    902902        vt_prof_cas(k) = vt_cas(k,it_cas2)                               &
    903      &          -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
     903                -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
    904904        dth_prof_cas(k) = dth_cas(k,it_cas2)                             &
    905      &          -frac*(dth_cas(k,it_cas2)-dth_cas(k,it_cas1))
     905                -frac*(dth_cas(k,it_cas2)-dth_cas(k,it_cas1))
    906906        hth_prof_cas(k) = hth_cas(k,it_cas2)                             &
    907      &          -frac*(hth_cas(k,it_cas2)-hth_cas(k,it_cas1))
     907                -frac*(hth_cas(k,it_cas2)-hth_cas(k,it_cas1))
    908908        vth_prof_cas(k) = vth_cas(k,it_cas2)                             &
    909      &          -frac*(vth_cas(k,it_cas2)-vth_cas(k,it_cas1))
     909                -frac*(vth_cas(k,it_cas2)-vth_cas(k,it_cas1))
    910910        dtrad_prof_cas(k) = dtrad_cas(k,it_cas2)                         &
    911      &          -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
     911                -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
    912912        dq_prof_cas(k) = dq_cas(k,it_cas2)                               &
    913      &          -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
     913                -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
    914914        hq_prof_cas(k) = hq_cas(k,it_cas2)                               &
    915      &          -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
     915                -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
    916916        vq_prof_cas(k) = vq_cas(k,it_cas2)                               &
    917      &          -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
     917                -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
    918918       uw_prof_cas(k) = uw_cas(k,it_cas2)                                &
    919      &          -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
     919                -frac*(uw_cas(k,it_cas2)-uw_cas(k,it_cas1))
    920920       vw_prof_cas(k) = vw_cas(k,it_cas2)                                &
    921      &          -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
     921                -frac*(vw_cas(k,it_cas2)-vw_cas(k,it_cas1))
    922922       q1_prof_cas(k) = q1_cas(k,it_cas2)                                &
    923      &          -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
     923                -frac*(q1_cas(k,it_cas2)-q1_cas(k,it_cas1))
    924924       q2_prof_cas(k) = q2_cas(k,it_cas2)                                &
    925      &          -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
     925                -frac*(q2_cas(k,it_cas2)-q2_cas(k,it_cas1))
    926926        enddo
    927927
    … …  
    932932!=====================================================================
    933933       SUBROUTINE interp2_case_vertical_std(play,plev,nlev_cas,plev_prof_cas                           &
    934      &         ,t_prof_cas,th_prof_cas,thv_prof_cas,thl_prof_cas                                       &
    935      &         ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas                              &
    936      &         ,ug_prof_cas,vg_prof_cas                                                                &
    937      &         ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas                    &
    938      &         ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas &     
    939      &         ,vitw_prof_cas,omega_prof_cas,tke_prof_cas                                              &
    940      &         ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas                &
    941      &         ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas &
    942      &         ,dth_prof_cas,hth_prof_cas,vth_prof_cas                                                 &
     934               ,t_prof_cas,th_prof_cas,thv_prof_cas,thl_prof_cas                                       &
     935               ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas                              &
     936               ,ug_prof_cas,vg_prof_cas                                                                &
     937               ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas                    &
     938               ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas &     
     939               ,vitw_prof_cas,omega_prof_cas,tke_prof_cas                                              &
     940               ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas                &
     941               ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas &
     942               ,dth_prof_cas,hth_prof_cas,vth_prof_cas                                                 &
    943943!
    944      &         ,t_mod_cas,theta_mod_cas,thv_mod_cas,thl_mod_cas                                        &
    945      &         ,qv_mod_cas,ql_mod_cas,qi_mod_cas,u_mod_cas,v_mod_cas                                   &
    946      &         ,ug_mod_cas,vg_mod_cas                                                                  &
    947      &         ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas                        &
    948      &         ,invtau_temp_nudg_mod_cas,invtau_qv_nudg_mod_cas,invtau_u_nudg_mod_cas,invtau_v_nudg_mod_cas                        &     
    949      &         ,w_mod_cas,omega_mod_cas,tke_mod_cas                                                    &
    950      &         ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas                      &
    951      &         ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas        &
    952      &         ,dth_mod_cas,hth_mod_cas,vth_mod_cas,mxcalc)
     944               ,t_mod_cas,theta_mod_cas,thv_mod_cas,thl_mod_cas                                        &
     945               ,qv_mod_cas,ql_mod_cas,qi_mod_cas,u_mod_cas,v_mod_cas                                   &
     946               ,ug_mod_cas,vg_mod_cas                                                                  &
     947               ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas                        &
     948               ,invtau_temp_nudg_mod_cas,invtau_qv_nudg_mod_cas,invtau_u_nudg_mod_cas,invtau_v_nudg_mod_cas                        &     
     949               ,w_mod_cas,omega_mod_cas,tke_mod_cas                                                    &
     950               ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas                      &
     951               ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas        &
     952               ,dth_mod_cas,hth_mod_cas,vth_mod_cas,mxcalc)
    953953 
    954954       implicit none
Note: See TracChangeset for help on using the changeset viewer.