subroutine qu2ll2 USE mod_extra USE mod_pvrbls USE mod_vrbls USE mod_pvrbls_ll implicit none include 'param_o.h' !!JNT include 'extra.comm' !! real,dimension(0:im+1,0:jm+1,nm,lsm)::fsl,usl,vsl,ull,vll,qsl,tsl,q2sl !! real,dimension(igm,jgm,lsm)::fsl_ll,hgt_ll,usp,vsp,qsl_ll,tsl_ll,q2sl_ll real, allocatable, dimension(:,:,:,:)::fsl,usl,vsl,ull,vll,qsl,tsl,q2sl real, allocatable, dimension(:,:,:)::fsl_ll,hgt_ll,usp,vsp,qsl_ll,tsl_ll,q2sl_ll real,dimension(jgm,lsm)::usp_avg real,dimension(jgm,lsm)::tsl_ll_avg !pvrbls_ll real,dimension(igm,jgm)::psurf_ll,tsfc_ll,pslp_ll,sfcshx_ll,sfclhx_ll,vegfrc_ll,z0_ll,albedo_ll !pvrbls_ll real, allocatable, dimension(:,:)::psurf_ll,tsfc_ll,pslp_ll,sfcshx_ll,sfclhx_ll,vegfrc_ll,z0_ll,albedo_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::tsfc,sfclhx,sfcshx,vegfrc,z0,albedo !! real,dimension(0:im+1,0:jm+1,nm) :: tsfc real, allocatable, dimension(:,:,:) :: tsfc !pvrbls real,dimension(0:im+1,0:jm+1,nm)::RSWIN,RSWOUT,RSWTOA,ASWIN,ASWOUT,ASWTOA,rswnet !pvrbls real,dimension(0:im+1,0:jm+1,nm)::RLWIN,RLWOUT,RLWTOA,ALWIN,ALWOUT,ALWTOA,rlwnet !!pvrbls_ll real,dimension(igm,jgm)::RSWIN_ll,RSWOUT_ll,RSWTOA_ll,ASWIN_ll,ASWOUT_ll,ASWTOA_ll,rswnet_ll !!pvrbls_ll real,dimension(igm,jgm)::RLWIN_ll,RLWOUT_ll,RLWTOA_ll,ALWIN_ll,ALWOUT_ll,ALWTOA_ll,rlwnet_ll !pvrbls_ll real, allocatable, dimension(:,:)::RSWIN_ll,RSWOUT_ll,RSWTOA_ll,ASWIN_ll,ASWOUT_ll,ASWTOA_ll,rswnet_ll !pvrbls_ll real, allocatable, dimension(:,:)::RLWIN_ll,RLWOUT_ll,RLWTOA_ll,ALWIN_ll,ALWOUT_ll,ALWTOA_ll,rlwnet_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm,lsm)::RSWTT,RLWTT !!pvrbls_ll real,dimension(igm,jgm,lsm)::RSWTT_ll,RLWTT_ll !pvrbls_ll real, allocatable, dimension(:,:,:)::RSWTT_ll,RLWTT_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::CUPPT,cfracl,cfracm,cfrach !!pvrbls_ll real,dimension(igm,jgm)::CUPPT_ll,cfracl_ll,cfracm_ll,cfrach_ll !pvrbls_ll real, allocatable, dimension(:,:)::CUPPT_ll,cfracl_ll,cfracm_ll,cfrach_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::acprec,cuprec,aprec,acsnow,prec !!pvrbls_ll real,dimension(igm,jgm)::acprec_ll,cuprec_ll,aprec_ll,acsnow_ll,prec_ll !pvrbls_ll real, allocatable, dimension(:,:)::acprec_ll,cuprec_ll,aprec_ll,acsnow_ll,prec_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::SMSTAV,BGROFF,SSROFF !! real,dimension(igm,jgm)::SMSTAV_ll,BGROFF_ll,SSROFF_ll real, allocatable, dimension(:,:)::SMSTAV_ll,BGROFF_ll,SSROFF_ll !pvrbls integer,parameter::nsoil=4 !pvrbls real,dimension(0:im+1,0:jm+1,nm,nsoil)::stc,smc !!pvrbls_ll real,dimension(igm,jgm,nsoil)::stc_ll,smc_ll !pvrbls_ll real, allocatable, dimension(:,:,:)::stc_ll,smc_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::tshltr,qshltr,albase,cnvtop,cnvbot,htop,hbot !!pvrbls_ll real,dimension(igm,jgm)::tshltr_ll,qshltr_ll,albase_ll,cnvtop_ll,cnvbot_ll,htop_ll,hbot_ll !pvrbls_ll real, allocatable, dimension(:,:)::tshltr_ll,qshltr_ll,albase_ll,cnvtop_ll,cnvbot_ll,htop_ll,hbot_ll !pvrbls integer,dimension(0:im+1,0:jm+1,nm)::IVGTYP,ISLTYP !!pvrbls_ll integer,dimension(igm,jgm)::IVGTYP_ll,ISLTYP_ll !pvrbls_ll integer, allocatable, dimension(:,:)::IVGTYP_ll,ISLTYP_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::acfrcv,acfrst,acsnom,cmc, & !pvrbls czen,czmean,epsr,glat,glon,grnflx,mxsnal,potevp,potflx !!pvrbls_ll real,dimension(igm,jgm)::acfrcv_ll,acfrst_ll,acsnom_ll,cmc_ll, & !pvrbls_ll real, allocatable, dimension(:,:)::acfrcv_ll,acfrst_ll,acsnom_ll,cmc_ll, & !pvrbls_ll czen_ll,czmean_ll,epsr_ll,glat_ll,glon_ll,grnflx_ll,mxsnal_ll,potevp_ll,potflx_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm,lm)::omgalf,q2,tcucn,train,div,rtop,htm,vtm, w, dwdt, z !!pvrbls_ll real,dimension(igm,jgm,lsm)::omgalf_ll,q2_ll,tcucn_ll,train_ll,div_ll,rtop_ll,htm_ll,vtm_ll, w_ll, dwdt_ll, z_ll !pvrbls_ll real, allocatable, dimension(:,:,:)::omgalf_ll,q2_ll,tcucn_ll,train_ll,div_ll,rtop_ll,htm_ll,vtm_ll, w_ll, dwdt_ll, z_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::qwbs,twbs,sm,fis,q10,th10,qz0,thz0,radot,sigt4, & !pvrbls sfcevp,subshx,snopcx,si,sno,sice,smstot,soiltb,tg,sr,akhs,akms,ths,plm !!pvrbls_ll real,dimension(igm,jgm)::qwbs_ll,twbs_ll,sm_ll,fis_ll,q10_ll,th10_ll,qz0_ll,thz0_ll, & !pvrbls_ll real, allocatable, dimension(:,:)::qwbs_ll,twbs_ll,sm_ll,fis_ll,q10_ll,th10_ll,qz0_ll,thz0_ll, & !pvrbls_ll radot_ll,sigt4_ll,sfcevp_ll,subshx_ll,snopcx_ll,si_ll,sno_ll,sice_ll,smstot_ll,soiltb_ll, & !pvrbls_ll tg_ll,sr_ll,akhs_ll,akms_ll,ths_ll,plm_ll !pvrbls real,dimension(0:im+1,0:jm+1,nm)::u10,v10,uz0,vz0,ustar,u10ll,v10ll,uz0ll,vz0ll !!pvrbls_ll real,dimension(0:im+1,0:jm+1,nm):: u10ll,v10ll,uz0ll,vz0ll !!pvrbls_ll real,dimension(igm,jgm)::u10_ll,v10_ll,uz0_ll,vz0_ll,ustar_ll !pvrbls_ll real, allocatable, dimension(:,:,:):: u10ll,v10ll,uz0ll,vz0ll !pvrbls_ll real, allocatable, dimension(:,:)::u10_ll,v10_ll,uz0_ll,vz0_ll,ustar_ll integer::l,i,j,n,k !+++++++++++++++++++ ! open(1,file='../data_out/tmpfld.dat',form='unformatted') ALLOCATE(fsl(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(usl(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(vsl(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(ull(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(vll(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(qsl(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(tsl(0:im+1,0:jm+1,nm,lsm)) ALLOCATE(q2sl(0:im+1,0:jm+1,nm,lsm)) ull = 0.0 vll = 0.0 ALLOCATE(tsfc(0:im+1,0:jm+1,nm)) open(1,file='tmpfld.dat',form='unformatted') !!! open(1,file='./ETAPOST/tmpfld.dat',form='unformatted') do l=1,lsm read(1)fsl(:,:,:,l) read(1)usl(:,:,:,l) read(1)vsl(:,:,:,l) read(1)qsl(:,:,:,l) !DRAGAN 23.09.'10 read(1)tsl(:,:,:,l) !!!avg2017 read(1)q2sl(:,:,:,l) !non-hydrostatic !!! read(1)w(:,:,:,l) !!! read(1)dwdt(:,:,:,l) !!! read(1)z(:,:,:,l) ! read(1)tshltr ! read(1)qshltr ! read(1)plm enddo close(1) !! do j=1,jm !! do i=1,im ! print *,i,j,q2sl(i,j,3,1) !! enddo !! enddo !++++++++++++++++++ open(1,file='tmpsfc.dat',form='unformatted') !!! open(1,file='./ETAPOST/tmpsfc.dat',form='unformatted') read(1)tsfc !jnt read(1)tshltr !jnt read(1)qshltr !jnt read(1)plm !jnt !jnt read(1)sfclhx !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !jnt read(1)sfcshx !!!!!!!!!!!! !jnt read(1)qwbs !jnt read(1)twbs !jnt !jnt read(1)grnflx !jnt read(1)vegfrc !!!!!!!!!!!!! !jnt read(1)IVGTYP !jnt read(1)ISLTYP !jnt read(1)glat !jnt read(1)glon !jnt read(1)sm !jnt! read(1)htm !jnt! read(1)vtm !jnt read(1)sice !jnt read(1)z0 !jnt read(1)omgalf !jnt read(1)albedo !jnt read(1)albase !jnt read(1)mxsnal !jnt !jnt read(1)czen !jnt read(1)czmean !jnt read(1)fis !jnt read(1)q10 !jnt read(1)th10 !jnt read(1)qz0 !jnt read(1)thz0 !jnt !jnt read(1)u10 !jnt read(1)v10 !jnt read(1)uz0 !jnt read(1)vz0 !jnt read(1)ustar !jnt !jnt read(1)ths !jnt read(1)q2 !jnt read(1)tcucn !jnt read(1)train !jnt read(1)div !jnt read(1)rtop !jnt !jnt !jnt close(1) !jnt!+++++++++++++++++++ !jnt open(2,file='./ETAPOST/sw_rdtn.dat',form='unformatted') !jnt read(2)RSWIN !jnt read(2)RSWOUT !jnt read(2)RSWTOA !!!!!!!!!!!!!!!!DRAGAN !jnt read(2)ASWIN !jnt read(2)ASWOUT !jnt !jnt read(2)ASWTOA !jnt read(2)rswnet !jnt read(2)RSWTT !jnt !jnt close(2) !jnt!+++++++++++++++++ !jnt open(3,file='./ETAPOST/lw_rdtn.dat',form='unformatted') !jnt read(3)RLWIN !jnt read(3)RLWOUT !jnt read(3)RLWTOA !!!!!!!!!!!!!!!!DRAGAN !jnt read(3)ALWIN !jnt read(3)ALWOUT !jnt !jnt read(3)ALWTOA !jnt read(3)rlwnet !jnt read(3)RLWTT !jnt read(3)epsr !jnt read(3)radot !jnt read(3)sigt4 !jnt close(3) !jnt !jnt !jnt!++++++++++++++++ !jnt !jnt !jnt open(4,file='./ETAPOST/clouds.dat',form='unformatted') !jnt read(4)CUPPT !jnt read(4)cfracl !jnt read(4)cfracm !!!!!!!!!!!!!!!!DRAGAN !jnt read(4)cfrach !jnt read(4)acfrcv !jnt read(4)acfrst !jnt !jnt read(4)cnvtop !jnt read(4)cnvbot !jnt! read(4)htop !jnt! read(4)hbot !jnt close(4) !jnt !jnt!++++++++++++++++ !jnt !jnt open(5,file='./ETAPOST/precANDsoil.dat',form='unformatted') !jnt read(5)acprec !jnt read(5)cuprec !jnt read(5)aprec !!!!!!!!!!!!!!!!DRAGAN !jnt read(5)acsnow !jnt read(5)acsnom !jnt !jnt read(5)prec !jnt !jnt read(5)stc(:,:,:,1) !jnt read(5)stc(:,:,:,4) !jnt read(5)smc(:,:,:,1) !jnt read(5)smc(:,:,:,4) !jnt !jnt !jnt read(5)SMSTAV !jnt read(5)BGROFF !jnt read(5)SSROFF !jnt read(5)cmc !jnt read(5)potevp !jnt read(5)potflx !jnt read(5)sfcevp !jnt read(5)subshx !jnt read(5)snopcx !jnt !jnt read(5)si !jnt read(5)sno !jnt !jnt read(5)smstot !jnt read(5)soiltb !!!! !jnt! read(5)tg !!!! !jnt read(5)sr !jnt read(5)akhs !jnt read(5)akms !jnt close(5) !+++++++++++++++++ ALLOCATE(fsl_ll(igm,jgm,lsm)) ALLOCATE(hgt_ll(igm,jgm,lsm)) ALLOCATE(usp(igm,jgm,lsm)) ALLOCATE(vsp(igm,jgm,lsm)) ALLOCATE(qsl_ll(igm,jgm,lsm)) ALLOCATE(tsl_ll(igm,jgm,lsm)) ALLOCATE(q2sl_ll(igm,jgm,lsm)) fsl_ll = 0.0 hgt_ll = 0.0 usp = 0.0 vsp = 0.0 qsl_ll = 0.0 tsl_ll = 0.0 q2sl_ll = 0.0 ALLOCATE(SMSTAV_ll(igm,jgm)) ALLOCATE(BGROFF_ll(igm,jgm)) ALLOCATE(SSROFF_ll(igm,jgm)) SMSTAV_ll = 0.0 BGROFF_ll = 0.0 SSROFF_ll = 0.0 print *,'tsl,30,14,4,',tsl(30,14,4,17),tsl(30,15,4,17),tsl(31,14,4,17), & tsl(31,15,4,17) print *,'tsl,30,14,4,',tsl(30,14,4,16),tsl(30,15,4,16),tsl(31,14,4,16), & tsl(31,15,4,16) !! do n=1,nm !! do i=1,im !! do j=1,jm ! print *,i,j,n,usl(i,j,n,1),vsl(i,j,n,1) !! enddo !! enddo !! enddo if(nm.eq.14)then ! print *,pslp(34,57,3),pslp(34,58,3),pslp(35,57,3),pslp(35,58,3) call octa2llH(fsl, fsl_ll, lsm) call octa2llH(qsl, qsl_ll, lsm) call octa2llH(tsl, tsl_ll, lsm) call octa2llH(q2sl, q2sl_ll, lsm) !!! call octa2llH(tshltr, tshltr_ll, 1) !!! call octa2llH(qshltr, qshltr_ll, 1) ! call octa2llH(plm, plm_ll, 1) call octa2llH(pint(:,:,:,lp1),psurf_ll,1) call octa2llH(pslp,pslp_ll,1) call octa2llH(tsfc,tsfc_ll,1) !++++++++++++++++++++++++++++++++++++++ call octa2llH(tshltr, tshltr_ll, 1) call octa2llH(qshltr, qshltr_ll, 1) call octa2llH(plm, plm_ll, 1) call octa2llH(sfclhx,sfclhx_ll,1) call octa2llH(sfcshx,sfcshx_ll,1) call octa2llH(qwbs,qwbs_ll,1) call octa2llH(twbs,twbs_ll,1) call octa2llH(grnflx,grnflx_ll,1) call octa2llH(vegfrc,vegfrc_ll,1) call octa2llH(IVGTYP,IVGTYP_ll,1) call octa2llH(ISLTYP,ISLTYP_ll,1) call octa2llH(glat,glat_ll,1) call octa2llH(glon,glon_ll,1) call octa2llH(sm,sm_ll,1) ! call octa2llH(htm,htm_ll,lsm) ! call octa2llH(vtm,vtm_ll,lsm) call octa2llH(sice,sice_ll,1) call octa2llH(z0,z0_ll,1) call octa2llH(omgalf,omgalf_ll,lsm) call octa2llH(albedo,albedo_ll,1) call octa2llH(albase,albase_ll,1) call octa2llH(mxsnal,mxsnal_ll,1) call octa2llH(czen,czen_ll,1) call octa2llH(czmean,czmean_ll,1) call octa2llH(fis,fis_ll,1) call octa2llH(q10,q10_ll,1) call octa2llH(th10,th10_ll,1) call octa2llH(qz0,qz0_ll,1) call octa2llH(thz0,thz0_ll,1) call octa2llH(ustar,ustar_ll,1) call octa2llH(ths,ths_ll,1) call octa2llH(q2,q2_ll,lsm) call octa2llH(tcucn,tcucn_ll,lsm) call octa2llH(train,train_ll,lsm) call octa2llH(div,div_ll,lsm) call octa2llH(rtop,rtop_ll,lsm) !++++++++++++++++++++++++++++++++++++ call octa2llH(RSWIN,RSWIN_ll,1) call octa2llH(RSWOUT,RSWOUT_ll,1) call octa2llH(RSWTOA,RSWTOA_ll,1) call octa2llH(ASWIN,ASWIN_ll,1) call octa2llH(ASWOUT,ASWOUT_ll,1) call octa2llH(ASWTOA,ASWTOA_ll,1) call octa2llH(rswnet,rswnet_ll,1) call octa2llH(RSWTT,RSWTT_ll,lsm) !+++++++++++++++++++++++++++++++++++++++ call octa2llH(RLWIN,RLWIN_ll,1) call octa2llH(RLWOUT,RLWOUT_ll,1) call octa2llH(RLWTOA,RLWTOA_ll,1) call octa2llH(ALWIN,ALWIN_ll,1) call octa2llH(ALWOUT,ALWOUT_ll,1) call octa2llH(ALWTOA,ALWTOA_ll,1) call octa2llH(rlwnet,rlwnet_ll,1) call octa2llH(RLWTT,RLWTT_ll,lsm) call octa2llH(epsr,epsr_ll,1) call octa2llH(radot,radot_ll,1) call octa2llH(sigt4,sigt4_ll,1) !++++++++++++++++++++++++++++++++++++++ call octa2llH(CUPPT,CUPPT_ll,1) call octa2llH(cfracl,cfracl_ll,1) call octa2llH(cfracm,cfracm_ll,1) call octa2llH(cfrach,cfrach_ll,1) call octa2llH(acfrcv,acfrcv_ll,1) call octa2llH(acfrst,acfrst_ll,1) !++++++++++++++++++++++++++++++++++++++ call octa2llH(cnvtop,cnvtop_ll,1) call octa2llH(cnvbot,cnvbot_ll,1) ! call octa2llH(htop,htop_ll,1) ! call octa2llH(hbot,hbot_ll,1) call octa2llH(acprec,acprec_ll,1) call octa2llH(cuprec,cuprec_ll,1) call octa2llH(aprec,aprec_ll,1) call octa2llH(acsnow,acsnow_ll,1) call octa2llH(acsnom,acsnom_ll,1) call octa2llH(prec,prec_ll,1) !++++++++++++++++++++++++++++++++++++++ call octa2llH(stc,stc_ll,1) call octa2llH(stc,stc_ll,4) call octa2llH(smc,smc_ll,1) call octa2llH(smc,smc_ll,4) call octa2llH(SMSTAV,SMSTAV_ll,1) call octa2llH(BGROFF,BGROFF_ll,1) call octa2llH(SSROFF,SSROFF_ll,1) call octa2llH(cmc,cmc_ll,1) call octa2llH(potevp,potevp_ll,1) call octa2llH(potflx,potflx_ll,1) call octa2llH(sfcevp,sfcevp_ll,1) call octa2llH(subshx,subshx_ll,1) call octa2llH(snopcx,snopcx_ll,1) call octa2llH(si,si_ll,1) call octa2llH(sno,sno_ll,1) call octa2llH(smstot,smstot_ll,1) call octa2llH(soiltb,soiltb_ll,1) ! call octa2llH(tg,tg_ll,1) call octa2llH(sr,sr_ll,1) call octa2llH(akhs,akhs_ll,1) call octa2llH(akms,akms_ll,1) !++++++++++++++++++++++++++++++++++++++ call llwinds_oc(usl, vsl, ull, vll,lsm) call bocovs_oc(ull,im,jm,lsm) call bocovs_oc(vll,im,jm,lsm) call avrv_oc(ull, usl,lsm) call avrv_oc(vll, vsl,lsm) usl = usl * .25 vsl = vsl * .25 do k=1,nm if(k.eq.1.or.k.eq.2.or.k.eq.7.or.k.eq.8.or.k.eq.9 & .or.k.eq.12)then usl(1 ,1 ,k,:)=usl(1 ,1 ,k,:)*4./3. vsl(1 ,1 ,k,:)=vsl(1 ,1 ,k,:)*4./3. endif if(k.eq.1.or.k.eq.4.or.k.eq.6.or.k.eq.7.or.k.eq.11 & .or.k.eq.14)then usl(im,1 ,k,:)=usl(im,1 ,k,:)*4./3. vsl(im,1 ,k,:)=vsl(im,1 ,k,:)*4./3. endif if(k.eq.2.or.k.eq.5.or.k.eq.6.or.k.eq.9.or.k.eq.13 & .or.k.eq.14)then usl(1 ,jm,k,:)=usl(1 ,jm,k,:)*4./3. vsl(1 ,jm,k,:)=vsl(1 ,jm,k,:)*4./3. endif if(k.eq.4.or.k.eq.5.or.k.eq.8.or.k.eq.11.or.k.eq.12 & .or.k.eq.13)then usl(im,jm,k,:)=usl(im,jm,k,:)*4./3. vsl(im,jm,k,:)=vsl(im,jm,k,:)*4./3. endif enddo call octa2llH(usl, usp, lsm) call octa2llH(vsl, vsp, lsm) !+++++++++++++++++++++++++++++ ! call llwinds_oc(u10, v10, u10ll, v10ll,1) ! call bocovs_oc(u10ll,im,jm,1) ! call bocovs_oc(v10ll,im,jm,1) ! call avrv_oc(u10ll, u10,1) ! call avrv_oc(v10ll, v10,1) ! u10 = u10 * .25 ! v10 = v10 * .25 ! do k=1,nm ! if(k.eq.1.or.k.eq.2.or.k.eq.7.or.k.eq.8.or.k.eq.9 & ! .or.k.eq.12)then ! u10(1 ,1 ,k)=u10(1 ,1 ,k)*4./3. ! v10(1 ,1 ,k)=v10(1 ,1 ,k)*4./3. ! endif ! ! if(k.eq.1.or.k.eq.4.or.k.eq.6.or.k.eq.7.or.k.eq.11 & ! .or.k.eq.14)then ! u10(im,1 ,k)=u10(im,1 ,k)*4./3. ! v10(im,1 ,k)=v10(im,1 ,k)*4./3. ! endif ! ! if(k.eq.2.or.k.eq.5.or.k.eq.6.or.k.eq.9.or.k.eq.13 & ! .or.k.eq.14)then ! u10(1 ,jm,k)=u10(1 ,jm,k)*4./3. ! v10(1 ,jm,k)=v10(1 ,jm,k)*4./3. ! endif ! ! if(k.eq.4.or.k.eq.5.or.k.eq.8.or.k.eq.11.or.k.eq.12 & ! .or.k.eq.13)then ! u10(im,jm,k)=u10(im,jm,k)*4./3. ! v10(im,jm,k)=v10(im,jm,k)*4./3. ! endif ! enddo call octa2llH(u10, u10_ll, 1) call octa2llH(v10, v10_ll, 1) !++++++++++++++++++++++++++++ ! call llwinds_oc(uz0, vz0, uz0ll, vz0ll,1) ! call bocovs_oc(uz0ll,im,jm,1) ! call bocovs_oc(vz0ll,im,jm,1) ! call avrv_oc(uz0ll, uz0,1) ! call avrv_oc(v10ll, vz0,1) ! uz0 = uz0 * .25 ! vz0 = vz0 * .25 ! do k=1,nm ! if(k.eq.1.or.k.eq.2.or.k.eq.7.or.k.eq.8.or.k.eq.9 & ! .or.k.eq.12)then ! uz0(1 ,1 ,k)=uz0(1 ,1 ,k)*4./3. ! vz0(1 ,1 ,k)=vz0(1 ,1 ,k)*4./3. ! endif ! if(k.eq.1.or.k.eq.4.or.k.eq.6.or.k.eq.7.or.k.eq.11 & ! .or.k.eq.14)then ! uz0(im,1 ,k)=uz0(im,1 ,k)*4./3. ! vz0(im,1 ,k)=vz0(im,1 ,k)*4./3. ! endif ! if(k.eq.2.or.k.eq.5.or.k.eq.6.or.k.eq.9.or.k.eq.13 & ! .or.k.eq.14)then ! uz0(1 ,jm,k)=uz0(1 ,jm,k)*4./3. ! vz0(1 ,jm,k)=vz0(1 ,jm,k)*4./3. ! endif ! if(k.eq.4.or.k.eq.5.or.k.eq.8.or.k.eq.11.or.k.eq.12 & ! .or.k.eq.13)then ! uz0(im,jm,k)=uz0(im,jm,k)*4./3. ! vz0(im,jm,k)=vz0(im,jm,k)*4./3. ! endif ! enddo call octa2llH(uz0, uz0_ll, 1) call octa2llH(vz0, vz0_ll, 1) !+++++ else if(nm.eq.6)then print *,'calling cube2llH - fsl ...' call cube2llH(fsl, fsl_ll, lsm) print *,'calling cube2llH - qsl ...' call cube2llH(qsl, qsl_ll, lsm) print *,'calling cube2llH - tsl ...' call cube2llH(tsl, tsl_ll, lsm) !!!avg2017 print *,'calling cube2llH - q2sl ...' !!!avg2017 call cube2llH(q2sl, q2sl_ll, lsm) !!! call cube2llH(tshltr, tshltr_ll, 1) !!! call cube2llH(qshltr, qshltr_ll, 1) ! call cube2llH(plm, plm_ll, 1) !!!avg2017 print *,'calling cube2llH - psurf_ll ...' !!!avg2017 call cube2llH(pint(:,:,:,lp1),psurf_ll,1) print *,'calling cube2llH - pslp_ll ...' call cube2llH(pslp,pslp_ll,1) print *,'calling cube2llH - tsfc_ll ...' call cube2llH(tsfc,tsfc_ll,1) !non-hydrostatic !!! print *,'calling cube2llH - w_ll ...' !!! call cube2llH(w, w_ll, lsm) !!! print *,'calling cube2llH - dwdt_ll ...' !!! call cube2llH(dwdt, dwdt_ll, lsm) !!! print *,'calling cube2llH - z_ll ...' !!! call cube2llH(z, z_ll, lsm) !++++++++++++++++++++++++++++++++++++++++++++++++++++++ print *,'calling cube2llH - tshltr_ll ...' call cube2llH(tshltr, tshltr_ll, 1) print *,'calling cube2llH - qshltr_ll ...' call cube2llH(qshltr, qshltr_ll, 1) print *,'calling cube2llH - plm_ll ...' call cube2llH(plm, plm_ll, 1) print *,'calling cube2llH - sfclhx_ll ...' call cube2llH(sfclhx,sfclhx_ll,1) print *,'calling cube2llH - sfcshx_ll ...' call cube2llH(sfcshx,sfcshx_ll,1) !!!avg2017 print *,'calling cube2llH - qwbs_ll ...' !!!avg2017 call cube2llH(qwbs,qwbs_ll,1) !!!avg2017 print *,'calling cube2llH - twbsb ...' !!!avg2017 call cube2llH(twbs,twbs_ll,1) !!!avg2017 call cube2llH(grnflx,grnflx_ll,1) call cube2llH(vegfrc,vegfrc_ll,1) !!!avg2017 call cube2llH(IVGTYP,IVGTYP_ll,1) !!!avg2017 call cube2llH(ISLTYP,ISLTYP_ll,1) !!!avg2017 print *,'calling cube2llH - glat_ll ...' !!!avg2017 call cube2llH(glat,glat_ll,1) !!!avg2017 call cube2llH(glon,glon_ll,1) call cube2llH(sm,sm_ll,1) ! call cube2llH(htm,htm_ll,lsm) ! call cube2llH(vtm,vtm_ll,lsm) call cube2llH(sice,sice_ll,1) call cube2llH(z0,z0_ll,1) print *,'calling cube2llH - omgalf_ll...' call cube2llH(omgalf,omgalf_ll,lsm) call cube2llH(albedo,albedo_ll,1) call cube2llH(albase,albase_ll,1) !!!avg2017 call cube2llH(mxsnal,mxsnal_ll,1) !!!avg2017 call cube2llH(czen,czen_ll,1) !!!avg2017 call cube2llH(czmean,czmean_ll,1) !!!avg2017 print *,'calling cube2llH - fis ...' !!!avg2017 call cube2llH(fis,fis_ll,1) !!!avg2017 call cube2llH(q10,q10_ll,1) !!!avg2017 call cube2llH(th10,th10_ll,1) !!!avg2017 call cube2llH(qz0,qz0_ll,1) !!!avg2017 call cube2llH(thz0,thz0_ll,1) !!!avg2017 call cube2llH(ustar,ustar_ll,1) !!!avg2017 print *,'calling cube2llH - ths_ll ...' !!!avg2017 call cube2llH(ths,ths_ll,1) !!!avg2017 call cube2llH(q2,q2_ll,lsm) !!!avg2017 call cube2llH(tcucn,tcucn_ll,lsm) !!!avg2017 call cube2llH(train,train_ll,lsm) !!!avg2017 call cube2llH(div,div_ll,lsm) !!!avg2017 call cube2llH(rtop,rtop_ll,lsm) !++++++++++++++++++++++++++++++++++++++++++++++++++++++ !!!avg2017 print *,'calling cube2llH - RSWIN_ll ...' !!!avg2017 call cube2llH(RSWIN,RSWIN_ll,1) !!!avg2017 print *,'calling cube2llH - 01 ...' !!!avg2017 call cube2llH(RSWOUT,RSWOUT_ll,1) !!!avg2017 print *,'calling cube2llH - 02 ...' !!!avg2017 call cube2llH(RSWTOA,RSWTOA_ll,1) print *,'calling cube2llH - 03 ...' call cube2llH(ASWIN,ASWIN_ll,1) print *,'calling cube2llH - 04 ...' call cube2llH(ASWOUT,ASWOUT_ll,1) !!!avg2017 print *,'calling cube2llH - 05 ...' !!!avg2017 call cube2llH(ASWTOA,ASWTOA_ll,1) !!!avg2017 print *,'calling cube2llH - 06 ...' !!!avg2017 call cube2llH(rswnet,rswnet_ll,1) !!!avg2017 print *,'calling cube2llH - 07 ...' !!!avg2017 call cube2llH(RSWTT,RSWTT_ll,lsm)!!!!!!!!!!!!!!!!!??? !+++++++++++++++++++++++++++++++++++++++ !!!avg2017 print *,'calling cube2llH - RLWIN_l ...' !!!avg2017 call cube2llH(RLWIN,RLWIN_ll,1) !!!avg2017 call cube2llH(RLWOUT,RLWOUT_ll,1) !!!avg2017 call cube2llH(RLWTOA,RLWTOA_ll,1) print *,'calling cube2llH - ALWIN_llf ...' call cube2llH(ALWIN,ALWIN_ll,1) call cube2llH(ALWOUT,ALWOUT_ll,1) call cube2llH(ALWTOA,ALWTOA_ll,1) !!!avg2017 call cube2llH(rlwnet,rlwnet_ll,1) !!!avg2017 call cube2llH(RLWTT,RLWTT_ll,lsm) !!!!!!!!!!??? !!!avg2017 call cube2llH(epsr,epsr_ll,1) !!!avg2017 print *,'calling cube2llH - radot_ll ...' !!!avg2017 call cube2llH(radot,radot_ll,1) !!!avg2017 call cube2llH(sigt4,sigt4_ll,1) !++++++++++++++++++++++++++++++++++++++ !!!avg2017 call cube2llH(CUPPT,CUPPT_ll,1) !!!avg2017 call cube2llH(cfracl,cfracl_ll,1) !!!avg2017 call cube2llH(cfracm,cfracm_ll,1) !!!avg2017 call cube2llH(cfrach,cfrach_ll,1) !!!avg2017 call cube2llH(acfrcv,acfrcv_ll,1) !!!avg2017 call cube2llH(acfrst,acfrst_ll,1) !!!avg2017 call cube2llH(cnvtop,cnvtop_ll,1) !!!avg2017 call cube2llH(cnvbot,cnvbot_ll,1) call cube2llH(htop,htop_ll,1) call cube2llH(hbot,hbot_ll,1) !++++++++++++++++++++++++++++++++++++++ print *,'calling cube2llH - acprec_l ...' call cube2llH(acprec,acprec_ll,1) call cube2llH(cuprec,cuprec_ll,1) call cube2llH(aprec,aprec_ll,1) call cube2llH(acsnow,acsnow_ll,1) !!!avg2017 call cube2llH(acsnom,acsnom_ll,1) call cube2llH(prec,prec_ll,1) !++++++++++++++++++++++++++++++++++++++ print *,'calling cube2llH - stc_ll ...' call cube2llH(stc,stc_ll,1) print *,'calling cube2llH - 01 ...' call cube2llH(stc,stc_ll,4) print *,'calling cube2llH - 02 ...' call cube2llH(smc,smc_ll,1) print *,'calling cube2llH - 03 ...' call cube2llH(smc,smc_ll,4) print *,'calling cube2llH - 04 ...' call cube2llH(SMSTAV,SMSTAV_ll,1) print *,'calling cube2llH - 05 ...' call cube2llH(BGROFF,BGROFF_ll,1) print *,'calling cube2llH - 06 ...' call cube2llH(SSROFF,SSROFF_ll,1) !!!avg2017 print *,'calling cube2llH - 07 ...' !!!avg2017 call cube2llH(cmc,cmc_ll,1) print *,'calling cube2llH - 08 ...' call cube2llH(potevp,potevp_ll,1) print *,'calling cube2llH - 09 ...' call cube2llH(potflx,potflx_ll,1) print *,'calling cube2llH - 10 ...' call cube2llH(sfcevp,sfcevp_ll,1) print *,'calling cube2llH - 11 ...' call cube2llH(subshx,subshx_ll,1) print *,'calling cube2llH - 12 ...' call cube2llH(snopcx,snopcx_ll,1) print *,'calling cube2llH - 13 ...' call cube2llH(si,si_ll,1) !!!avg2017 print *,'calling cube2llH - 14 ...' !!!avg2017 call cube2llH(sno,sno_ll,1) !!!avg2017 print *,'calling cube2llH - smstot_ll ...' !!!avg2017 call cube2llH(smstot,smstot_ll,1) !!!avg2017 call cube2llH(soiltb,soiltb_ll,1) ! call cube2llH(tg,tg_ll,1) !!!avg2017 call cube2llH(sr,sr_ll,1) !!!avg2017 call cube2llH(akhs,akhs_ll,1) !!!avg2017 call cube2llH(akms,akms_ll,1) !++++++++++++++++++++++++++++++++++++++ print *,'calling llwinds - ...' call llwinds(usl, vsl, ull, vll,lsm) call bocovs(ull,im,jm,lsm) call bocovs(vll,im,jm,lsm) call avrv(ull, usl,lsm) call avrv(vll, vsl,lsm) usl = usl * .25 vsl = vsl * .25 do k=1,nm usl(1 ,1 ,k,:)=usl(1 ,1 ,k,:)*4./3. vsl(1 ,1 ,k,:)=vsl(1 ,1 ,k,:)*4./3. usl(im,1 ,k,:)=usl(im,1 ,k,:)*4./3. vsl(im,1 ,k,:)=vsl(im,1 ,k,:)*4./3. usl(1 ,jm,k,:)=usl(1 ,jm,k,:)*4./3. vsl(1 ,jm,k,:)=vsl(1 ,jm,k,:)*4./3. usl(im,jm,k,:)=usl(im,jm,k,:)*4./3. vsl(im,jm,k,:)=vsl(im,jm,k,:)*4./3. enddo print *,'calling cube2llH - usp,vsp ...' call cube2llH(usl, usp, lsm) call cube2llH(vsl, vsp, lsm) !++++++++++++++++++++++++++++++++++++++ ! call llwinds(u10, v10, u10ll, v10ll,1) ! call bocovs(u10ll,im,jm,1) ! call bocovs(v10ll,im,jm,1) ! call avrv(u10ll, u10,1) ! call avrv(v10ll, v10,1) ! u10 = u10 * .25 ! v10 = v10 * .25 ! do k=1,nm ! u10(1 ,1 ,k)=u10(1 ,1 ,k)*4./3. ! v10(1 ,1 ,k)=v10(1 ,1 ,k)*4./3. ! u10(im,1 ,k)=u10(im,1 ,k)*4./3. ! v10(im,1 ,k)=v10(im,1 ,k)*4./3. ! u10(1 ,jm,k)=u10(1 ,jm,k)*4./3. ! v10(1 ,jm,k)=v10(1 ,jm,k)*4./3. ! u10(im,jm,k)=u10(im,jm,k)*4./3. ! v10(im,jm,k)=v10(im,jm,k)*4./3. ! enddo ! call cube2llH(u10, u10_ll, 1) call cube2llH(v10, v10_ll, 1) !+++++++++++++++++++++++++++++++++++++ ! call llwinds(uz0, vz0, uz0ll, vz0ll,1) ! call bocovs(uz0ll,im,jm,1) ! call bocovs(vz0ll,im,jm,1) ! call avrv(uz0ll, uz0,1) ! call avrv(vz0ll, vz0,1) ! uz0 = uz0 * .25 ! vz0 = vz0 * .25 ! do k=1,nm ! uz0(1 ,1 ,k)=uz0(1 ,1 ,k)*4./3. ! vz0(1 ,1 ,k)=vz0(1 ,1 ,k)*4./3. ! uz0(im,1 ,k)=uz0(im,1 ,k)*4./3. ! vz0(im,1 ,k)=vz0(im,1 ,k)*4./3. ! uz0(1 ,jm,k)=uz0(1 ,jm,k)*4./3. ! vz0(1 ,jm,k)=vz0(1 ,jm,k)*4./3. ! uz0(im,jm,k)=uz0(im,jm,k)*4./3. ! vz0(im,jm,k)=vz0(im,jm,k)*4./3. ! enddo ! !!!avg2017 print *,'calling cube2llH - uz0 ...' !!!avg2017 call cube2llH(uz0, uz0_ll, 1) !!!avg2017 print *,'calling cube2llH - vz0 ...' !!!avg2017 call cube2llH(vz0, vz0_ll, 1) !+++++++++++++++++++++++++++++++++++++ endif print *,' Endif nm.eq.6 ' print *,' update hgt_ll ' hgt_ll=fsl_ll/9.80616 ! print *,usp(82,44,17) print *,' update usp_avg, tsl_ll_avg ' do l=1,lsm do j=1,jgm usp_avg(j,l)=0. tsl_ll_avg(j,l)=0. do i=1,igm ! do i=14,250 usp_avg(j,l)=usp_avg(j,l)+usp(i,j,l) tsl_ll_avg(j,l)=tsl_ll_avg(j,l)+tsl_ll(i,j,l) if(l.eq.17.and.j.eq.44)then ! print *,i,j,l,tsl_ll(i,j,l),usp_avg(j,l),usp(i,j,l) endif enddo ! do i=335,360 !igm ! usp_avg(j,l)=usp_avg(j,l)+usp(i,j,l) ! tsl_ll_avg(j,l)=tsl_ll_avg(j,l)+tsl_ll(i,j,l) ! if(l.eq.17.and.j.eq.44)then ! print *,i,j,l,tsl_ll(i,j,l),usp_avg(j,l),usp(i,j,l) ! endif ! enddo usp_avg(j,l)=usp_avg(j,l)/igm ! usp_avg(j,l)=usp_avg(j,l)/237 tsl_ll_avg(j,l)=tsl_ll_avg(j,l)/igm ! tsl_ll_avg(j,l)=tsl_ll_avg(j,l)/237 enddo enddo !! do l=lsm-1,lsm !! do j=90,91 ! print *,j,l,tsl_ll_avg(j,l) !! enddo !! enddo != = = = = = = = = = = = = = = = = = = = = = = = = = = !VARIABLES WHICH GO TO THE GLOBRUN***.DAT FILES !DRAGAN 21.09.10. != = = = = = = = = = = = = = = = = = = = = = = = = = = print *, 'writeout to unit 9 .... ' write(9)tsfc_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)psurf_ll(:,:) !DRAGAN 21.09. write(9)pslp_ll(:,:) !DRAGAN 21.09. !!! write(9)tshltr_ll(:,:) !DRAGAN 21.09. !!! write(9)qshltr_ll(:,:) !DRAGAN 21.09. ! write(9)plm_ll(:,:) !++++++++++++++++++++++++++++++++++++++++++++ write(9)tshltr_ll(:,:) !DRAGAN 21.09. write(9)qshltr_ll(:,:) !DRAGAN 21.09. write(9)plm_ll(:,:) write(9)sfclhx_ll(:,:) !DRAGAN 21.09. write(9)sfcshx_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)qwbs_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)twbs_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)grnflx_ll(:,:) !DRAGAN 21.09. write(9)vegfrc_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)IVGTYP_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)ISLTYP_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)glat_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)glon_ll(:,:) !DRAGAN 21.09. write(9)sm_ll(:,:) !DRAGAN 21.09. ! do l=1,lsm ! write(9)htm_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm ! write(9)vtm_ll(:,:,l) !DRAGAN 17.09. ! enddo write(9)sice_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)z0_ll(:,:) !DRAGAN 21.09. write(9)albedo_ll(:,:) !DRAGAN 21.09. write(9)albase_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)mxsnal_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)czen_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)czmean_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)fis_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)q10_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)th10_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)qz0_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)thz0_ll(:,:) !DRAGAN 21.09. !!!! write(9)u10_ll(:,:) !DRAGAN 21.09. !!!! write(9)v10_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)uz0_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)vz0_ll(:,:) !DRAGAN 21.09. !!!!!! write(9)ustar_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)ths_ll(:,:) !DRAGAN 21.09. !++++++++++++++++++++++++++++++++++++++++++++ !!!!!!!! write(9)RSWIN_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)RSWOUT_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)RSWTOA_ll(:,:) !DRAGAN 21.09. write(9)ASWIN_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)ASWOUT_ll(:,:) !DRAGAN 21.09. !!!avg2017 !!!!!!!! write(9)ASWTOA_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)rswnet_ll(:,:) !DRAGAN 21.09. !++++++++++++ !!!!!!!! write(9)RLWIN_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)RLWOUT_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)RLWTOA_ll(:,:) !DRAGAN 21.09. write(9)ALWIN_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)ALWOUT_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)ALWTOA_ll(:,:) !DRAGAN 21.09. !!!avg2017 !!!!!!!! write(9)rlwnet_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)epsr_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)radot_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)sigt4_ll(:,:) !DRAGAN 21.09. !++++++++++++ !!!avg2017 write(9)CUPPT_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cfracl_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cfracm_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cfrach_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)acfrcv_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)acfrst_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cnvtop_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cnvbot_ll(:,:) !DRAGAN 21.09. write(9)htop_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)hbot_ll(:,:) !DRAGAN 21.09. !!!avg2017 !++++++++++++ write(9)acprec_ll(:,:) !DRAGAN 21.09. write(9)cuprec_ll(:,:) !DRAGAN 21.09. write(9)aprec_ll(:,:) !DRAGAN 21.09. write(9)acsnow_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)acsnom_ll(:,:) !DRAGAN 21.09. write(9)prec_ll(:,:) !DRAGAN 21.09. !++++++++++++ write(9)stc_ll(:,:,1) !DRAGAN 21.09. write(9)stc_ll(:,:,4) !DRAGAN 21.09. write(9)smc_ll(:,:,1) !DRAGAN 21.09. write(9)smc_ll(:,:,4) !DRAGAN 21.09. write(9)SMSTAV_ll(:,:) !DRAGAN 21.09. write(9)BGROFF_ll(:,:) !DRAGAN 21.09. write(9)SSROFF_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)cmc_ll(:,:) !DRAGAN 21.09. write(9)potevp_ll(:,:) !DRAGAN 21.09. write(9)potflx_ll(:,:) !DRAGAN 21.09. write(9)sfcevp_ll(:,:) !DRAGAN 21.09. write(9)subshx_ll(:,:) !DRAGAN 21.09. write(9)snopcx_ll(:,:) !DRAGAN 21.09. write(9)si_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)sno_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)smstot_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)soiltb_ll(:,:) !DRAGAN 21.09. ! write(9)tg_ll(:,:) !DRAGAN 21.09. !!!avg2017 write(9)sr_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)akhs_ll(:,:) !DRAGAN 21.09. !!!!!!!! write(9)akms_ll(:,:) !DRAGAN 21.09. !++++++++++++ ! do l=1,lsm !!!!!!!! write(9)RSWTT_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm !!!!!!!! write(9)RLWTT_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm !!!!!!!! write(9)omgalf_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm !!!!!!!! write(9)q2_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm !!!!!!!! write(9)tcucn_ll(:,:,l) !DRAGAN 17.09. ! enddo ! do l=1,lsm !!!!!!!! write(9)train_ll(:,:,l) !DRAGAN 17.09. ! enddo !!!avg2017 do l=1,lsm !!!avg2017 write(9)div_ll(:,:,l) !DRAGAN 17.09. !!!avg2017 enddo !!!avg2017 do l=1,lsm !!!avg2017 write(9)rtop_ll(:,:,l) !DRAGAN 17.09. !!!avg2017 enddo !++++++++++++++ !!! do l=1,lsm !!! write(9)w_ll(:,:,l) !DRAGAN 2016 !!! enddo !!! do l=1,lsm !!! write(9)dwdt_ll(:,:,l) !DRAGAN 2016 !!! enddo !!! do l=1,lsm !!! write(9)z_ll(:,:,l) !DRAGAN 2016 !!! enddo do l=1,lsm write(9)hgt_ll(:,:,l) !DRAGAN 17.09. enddo do l=1,lsm write(9)tsl_ll(:,:,l) !DRAGAN 21.09. enddo !!!avg2017 do l=1,lsm !!!avg2017 write(9)fsl_ll(:,:,l) !DRAGAN 21.09. !!!avg2017 enddo do l=1,lsm write(9)qsl_ll(:,:,l) !DRAGAN 21.09. enddo !!!avg2017 do l=1,lsm !!!avg2017 write(9)q2sl_ll(:,:,l) !DRAGAN 21.09. !!!avg2017 enddo do l=1,lsm write(9)usp(:,:,l) !DRAGAN 21.09. enddo do l=1,lsm write(9)vsp(:,:,l) !DRAGAN 21.09. enddo != = = = = = = = = = = = = = = = = = = = = = = = = = ! write(10)hgt_ll(:,:,10) !DRAGAN 17.09 ! write(10)hgt_ll(:,:,4) !DRAGAN 17.09 write(9)tsl_ll_avg !DRAGAN 17.09. write(9)usp_avg !DRAGAN 20.09. print *,'igm=',igm, 'jgm=',jgm DEALLOCATE(fsl,usl,vsl,ull,vll,qsl,tsl,q2sl,tsfc) DEALLOCATE(fsl_ll,hgt_ll,usp,vsp,qsl_ll,tsl_ll,q2sl_ll) ! DEALLOCATE(usl,vsl,ull,vll,qsl,tsl,tsfc) ! DEALLOCATE(hgt_ll,usp,vsp,qsl_ll,tsl_ll) DEALLOCATE(SMSTAV_ll,BGROFF_ll,SSROFF_ll) print *, 'leaving qu2ll2 .... ' end subroutine qu2ll2