| 1 | ! Ecriture des fichiers de sorties COSP |
|---|
| 2 | ! Sorties journalierres |
|---|
| 3 | ! Abderrahmane Idelkadi Septembre 2009 |
|---|
| 4 | |
|---|
| 5 | IF (MOD(itap,NINT(freq_COSP/dtime)).EQ.0) THEN |
|---|
| 6 | |
|---|
| 7 | itau_wcosp = itau_phy + itap |
|---|
| 8 | |
|---|
| 9 | ! Sorties LIDAR |
|---|
| 10 | if (cfg%Llidar_sim) then |
|---|
| 11 | if (cfg%Lcllcalipso) then |
|---|
| 12 | CALL histwrite_phy(nid_hf_cosp,"cllcalipso",itau_wcosp,stlidar%cldlayer(:,1)) |
|---|
| 13 | endif |
|---|
| 14 | if (cfg%Lclhcalipso) then |
|---|
| 15 | CALL histwrite_phy(nid_hf_cosp,"clhcalipso",itau_wcosp,stlidar%cldlayer(:,3)) |
|---|
| 16 | endif |
|---|
| 17 | if (cfg%Lclmcalipso) then |
|---|
| 18 | CALL histwrite_phy(nid_hf_cosp,"clmcalipso",itau_wcosp,stlidar%cldlayer(:,2)) |
|---|
| 19 | endif |
|---|
| 20 | if (cfg%Lcltcalipso) then |
|---|
| 21 | CALL histwrite_phy(nid_hf_cosp,"cltcalipso",itau_wcosp,stlidar%cldlayer(:,4)) |
|---|
| 22 | endif |
|---|
| 23 | if (cfg%Lclcalipso) then |
|---|
| 24 | CALL histwrite_phy(nid_hf_cosp,"clcalipso",itau_wcosp,stlidar%lidarcld) |
|---|
| 25 | endif |
|---|
| 26 | if (cfg%Lcfad_lidarsr532) then |
|---|
| 27 | do ii=1,SR_BINS |
|---|
| 28 | CALL histwrite_phy(nid_hf_cosp,"cfad_lidarsr532_"//chcol(ii),itau_wcosp,stlidar%cfad_sr(:,ii,:)) |
|---|
| 29 | enddo |
|---|
| 30 | endif |
|---|
| 31 | if (cfg%Lparasol_refl) then |
|---|
| 32 | CALL histwrite_phy(nid_hf_cosp,"parasol_refl",itau_wcosp,stlidar%parasolrefl) |
|---|
| 33 | endif |
|---|
| 34 | if (cfg%Latb532) then |
|---|
| 35 | do ii=1,Ncolumns |
|---|
| 36 | CALL histwrite_phy(nid_hf_cosp,"atb532_"//chcol(ii),itau_wcosp,sglidar%beta_tot(:,ii,:)) |
|---|
| 37 | enddo |
|---|
| 38 | endif |
|---|
| 39 | if (cfg%Lbeta_mol532) then |
|---|
| 40 | CALL histwrite_phy(nid_hf_cosp,"beta_mol532",itau_wcosp,sglidar%beta_mol) |
|---|
| 41 | endif |
|---|
| 42 | endif ! Lidar |
|---|
| 43 | |
|---|
| 44 | ! Sorties RADAR |
|---|
| 45 | !Attention A FAIRE |
|---|
| 46 | ! if (cfg%Lradar_sim) then |
|---|
| 47 | ! print*,'Ecriture sorties Radar' |
|---|
| 48 | ! if (cfg%Lcfad_dbze94) then |
|---|
| 49 | ! print*,'Ecriture de cfad_dbze94.nc ' |
|---|
| 50 | ! A revoir l axe vertical Nlvgrid |
|---|
| 51 | ! do ii=1,DBZE_BINS |
|---|
| 52 | ! dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5) |
|---|
| 53 | ! enddo |
|---|
| 54 | ! call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, & |
|---|
| 55 | ! x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze) |
|---|
| 56 | ! endif |
|---|
| 57 | ! if (cfg%Lclcalipso2) then |
|---|
| 58 | ! call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', & |
|---|
| 59 | ! nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud) |
|---|
| 60 | ! endif |
|---|
| 61 | ! if (cfg%Ldbze94) then |
|---|
| 62 | ! do ii=1,Ncolumns |
|---|
| 63 | ! xcol(ii)=float(i) |
|---|
| 64 | ! enddo |
|---|
| 65 | ! call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, & |
|---|
| 66 | ! x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot) |
|---|
| 67 | ! endif |
|---|
| 68 | ! if (cfg%Lcltlidarradar) then |
|---|
| 69 | ! call write_netcdf2d('cltlidarradar.nc','cltlidarradar', & |
|---|
| 70 | ! nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc) |
|---|
| 71 | ! endif |
|---|
| 72 | ! endif ! Radar |
|---|
| 73 | |
|---|
| 74 | ! Sorties MISR |
|---|
| 75 | !Attention A FAIRE |
|---|
| 76 | ! if (cfg%Lmisr_sim) then |
|---|
| 77 | ! print*,'Ecriture sorties Misr' |
|---|
| 78 | ! call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, & |
|---|
| 79 | ! x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR) |
|---|
| 80 | ! endif |
|---|
| 81 | |
|---|
| 82 | ! Sorties ISCCP |
|---|
| 83 | if (cfg%Lisccp_sim) then |
|---|
| 84 | if (cfg%Lclisccp2) then |
|---|
| 85 | do ii=1,7 |
|---|
| 86 | CALL histwrite_phy(nid_hf_cosp,"clisccp2_"//chcol(ii),itau_wcosp,isccp%fq_isccp(:,ii,:)) |
|---|
| 87 | enddo |
|---|
| 88 | endif |
|---|
| 89 | if (cfg%Lboxtauisccp) then |
|---|
| 90 | CALL histwrite_phy(nid_hf_cosp,"boxtauisccp",itau_wcosp,isccp%boxtau) |
|---|
| 91 | endif |
|---|
| 92 | if (cfg%Lboxptopisccp) then |
|---|
| 93 | CALL histwrite_phy(nid_hf_cosp,"boxptopisccp",itau_wcosp,isccp%boxptop) |
|---|
| 94 | endif |
|---|
| 95 | if (cfg%Ltclisccp) then |
|---|
| 96 | CALL histwrite_phy(nid_hf_cosp,"tclisccp",itau_wcosp,isccp%totalcldarea) |
|---|
| 97 | endif |
|---|
| 98 | if (cfg%Lctpisccp) then |
|---|
| 99 | CALL histwrite_phy(nid_hf_cosp,"ctpisccp",itau_wcosp,isccp%meanptop) |
|---|
| 100 | |
|---|
| 101 | endif |
|---|
| 102 | if (cfg%Ltauisccp) then |
|---|
| 103 | CALL histwrite_phy(nid_hf_cosp,"tauisccp",itau_wcosp,isccp%meantaucld) |
|---|
| 104 | endif |
|---|
| 105 | if (cfg%Lalbisccp) then |
|---|
| 106 | CALL histwrite_phy(nid_hf_cosp,"albisccp",itau_wcosp,isccp%meanalbedocld) |
|---|
| 107 | endif |
|---|
| 108 | if (cfg%Lmeantbisccp) then |
|---|
| 109 | CALL histwrite_phy(nid_hf_cosp,"meantbisccp",itau_wcosp,isccp%meantb) |
|---|
| 110 | endif |
|---|
| 111 | if (cfg%Lmeantbclrisccp) then |
|---|
| 112 | CALL histwrite_phy(nid_hf_cosp,"meantbclrisccp",itau_wcosp,isccp%meantbclr) |
|---|
| 113 | endif |
|---|
| 114 | endif ! Isccp |
|---|
| 115 | |
|---|
| 116 | ! if (ok_sync) then |
|---|
| 117 | !$OMP MASTER |
|---|
| 118 | call histsync(nid_hf_cosp) |
|---|
| 119 | !$OMP END MASTER |
|---|
| 120 | ! endif |
|---|
| 121 | |
|---|
| 122 | ENDIF ! if freq_COSP |
|---|
