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

Remove useless & at beginning of continued line

File:
1 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
Note: See TracChangeset for help on using the changeset viewer.