| 1 | ! Abderrahmane Idelkadi Septebmre 2009 |
|---|
| 2 | ! Sorties journalieres de COSP |
|---|
| 3 | ! Pour l'instant sorties Lidar et ISCCP |
|---|
| 4 | ! |
|---|
| 5 | ! sorties par jour |
|---|
| 6 | ! |
|---|
| 7 | !$OMP MASTER |
|---|
| 8 | zstoday = ecrit_day |
|---|
| 9 | zout = freq_COSP |
|---|
| 10 | ! |
|---|
| 11 | ! PRINT*, 'La frequence de sortie hf3d est de ', ecrit_hf |
|---|
| 12 | ! |
|---|
| 13 | |
|---|
| 14 | idayref = day_ref |
|---|
| 15 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
|---|
| 16 | |
|---|
| 17 | CALL histbeg_phy("histdayCOSP",itau_phy,zjulian,dtime,nhori,nid_day_cosp) |
|---|
| 18 | |
|---|
| 19 | ! Definition de l'axe vertical |
|---|
| 20 | CALL histvert(nid_day_cosp,"height","height","m",Nlevout,vgrid%z,nvert) |
|---|
| 21 | print*,'Ok height Nlevout, height =',Nlevout,vgrid%z |
|---|
| 22 | CALL histvert(nid_day_cosp,"height_mlev","height_mlev","m",Nlevlmdz,vgrid%mz,nvertm) |
|---|
| 23 | print*,'Ok height_mlev Nlevout, height_mlev =',Nlevout,vgrid%mz |
|---|
| 24 | ! CALL histvert(nid_day_cosp,"presnivs","Vertical levels","mb",Nlevout,presnivs,nvert) |
|---|
| 25 | |
|---|
| 26 | CALL histvert(nid_day_cosp,"sza","solar_zenith_angle","degrees",PARASOL_NREFL,PARASOL_SZA,nvertp) |
|---|
| 27 | |
|---|
| 28 | CALL histvert(nid_day_cosp,"pressure2","pressure","mb",7,ISCCP_PC,nvertisccp) |
|---|
| 29 | |
|---|
| 30 | CALL histvert(nid_day_cosp,"column","column","count",Ncolumns,column_ax,nvertcol) |
|---|
| 31 | |
|---|
| 32 | ! Sorties LIDAR |
|---|
| 33 | if (cfg%Llidar_sim) then |
|---|
| 34 | if (cfg%Lcllcalipso) then |
|---|
| 35 | CALL histdef(nid_day_cosp, "cllcalipso", & |
|---|
| 36 | "Lidar Low-level Cloud Fraction", "1", & |
|---|
| 37 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 38 | "ave(X)", zout,zstoday) |
|---|
| 39 | endif |
|---|
| 40 | if (cfg%Lclhcalipso) then |
|---|
| 41 | CALL histdef(nid_day_cosp, "clhcalipso", & |
|---|
| 42 | "Lidar High-level Cloud Fraction", "1", & |
|---|
| 43 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 44 | "ave(X)", zout,zstoday) |
|---|
| 45 | endif |
|---|
| 46 | if (cfg%Lclmcalipso) then |
|---|
| 47 | CALL histdef(nid_day_cosp, "clmcalipso", & |
|---|
| 48 | "Lidar Mid-level Cloud Fraction", "1", & |
|---|
| 49 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 50 | "ave(X)", zout,zstoday) |
|---|
| 51 | endif |
|---|
| 52 | if (cfg%Lcltcalipso) then |
|---|
| 53 | CALL histdef(nid_day_cosp, "cltcalipso", & |
|---|
| 54 | "Lidar Total Cloud Fraction", "1", & |
|---|
| 55 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 56 | "ave(X)", zout,zstoday) |
|---|
| 57 | endif |
|---|
| 58 | if (cfg%Lclcalipso) then |
|---|
| 59 | CALL histdef(nid_day_cosp, "clcalipso", & |
|---|
| 60 | "Lidar Cloud Fraction (532 nm)", "1", & |
|---|
| 61 | iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, & |
|---|
| 62 | "ave(X)", zout,zstoday) |
|---|
| 63 | endif |
|---|
| 64 | if (cfg%Lcfad_lidarsr532) then |
|---|
| 65 | do ii=1,SR_BINS |
|---|
| 66 | CALL histdef(nid_day_cosp, "cfad_lidarsr532_"//chcol(ii), & |
|---|
| 67 | "Lidar Scattering Ratio CFAD (532 nm)","1", & |
|---|
| 68 | iim,jj_nb,nhori, Nlevout,1,Nlevout,nvert, 32, & |
|---|
| 69 | "ave(X)", zout,zstoday) |
|---|
| 70 | enddo |
|---|
| 71 | endif |
|---|
| 72 | if (cfg%Lparasol_refl) then |
|---|
| 73 | CALL histdef(nid_day_cosp, "parasol_refl", & |
|---|
| 74 | "PARASOL-like mono-directional reflectance","1", & |
|---|
| 75 | iim,jj_nb,nhori, PARASOL_NREFL,1, PARASOL_NREFL, nvertp,32, & |
|---|
| 76 | "ave(X)", zout,zstoday) |
|---|
| 77 | endif |
|---|
| 78 | if (cfg%Latb532) then |
|---|
| 79 | do ii=1,Ncolumns |
|---|
| 80 | CALL histdef(nid_day_cosp, "atb532_"//chcol(ii), & |
|---|
| 81 | "Lidar Attenuated Total Backscatter (532 nm)","1", & |
|---|
| 82 | iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, & |
|---|
| 83 | "ave(X)", zout,zstoday) |
|---|
| 84 | enddo |
|---|
| 85 | endif |
|---|
| 86 | if (cfg%Lbeta_mol532) then |
|---|
| 87 | CALL histdef(nid_day_cosp, "beta_mol532", & |
|---|
| 88 | "Lidar Molecular Backscatter (532 nm)","m-1 sr-1", & |
|---|
| 89 | iim,jj_nb,nhori, Nlevlmdz,1,Nlevlmdz,nvertm, 32, & |
|---|
| 90 | "ave(X)", zout,zstoday) |
|---|
| 91 | endif |
|---|
| 92 | endif ! Lidar |
|---|
| 93 | |
|---|
| 94 | ! Sorties RADAR |
|---|
| 95 | !Attention A FAIRE |
|---|
| 96 | ! if (cfg%Lradar_sim) then |
|---|
| 97 | ! print*,'Ecriture sorties Radar' |
|---|
| 98 | ! if (cfg%Lcfad_dbze94) then |
|---|
| 99 | ! print*,'Ecriture de cfad_dbze94.nc ' |
|---|
| 100 | ! A revoir l axe vertical Nlvgrid |
|---|
| 101 | ! do ii=1,DBZE_BINS |
|---|
| 102 | ! dbze_ax(ii) = CFAD_ZE_MIN + CFAD_ZE_WIDTH*(ii - 0.5) |
|---|
| 103 | ! enddo |
|---|
| 104 | ! call write_netcdf4d('cfad_dbze94.nc',use_vgrid,nlon,nlat,Nlevout,DBZE_BINS, & |
|---|
| 105 | ! x,y,out_levs,dbze_ax,i,ndays,time,stradar%cfad_ze) |
|---|
| 106 | ! endif |
|---|
| 107 | ! if (cfg%Lclcalipso2) then |
|---|
| 108 | ! call write_netcdf3d('clcalipso2.nc',use_vgrid,'clcalipso2', & |
|---|
| 109 | ! nlon,nlat,Nlevout,x,y,out_levs,i,ndays,time,stradar%lidar_only_freq_cloud) |
|---|
| 110 | ! endif |
|---|
| 111 | ! if (cfg%Ldbze94) then |
|---|
| 112 | ! do ii=1,Ncolumns |
|---|
| 113 | ! xcol(ii)=float(i) |
|---|
| 114 | ! enddo |
|---|
| 115 | ! call write_netcdf4d('dbze94.nc',use_vgrid,nlon,nlat,Nlevout,Ncolumns, & |
|---|
| 116 | ! x,y,out_levs,xcol,i,ndays,time,sgradar%Ze_tot) |
|---|
| 117 | ! endif |
|---|
| 118 | ! if (cfg%Lcltlidarradar) then |
|---|
| 119 | ! call write_netcdf2d('cltlidarradar.nc','cltlidarradar', & |
|---|
| 120 | ! nlon,nlat,x,y,i,ndays,time,stradar%radar_lidar_tcc) |
|---|
| 121 | ! endif |
|---|
| 122 | ! endif ! Radar |
|---|
| 123 | |
|---|
| 124 | ! Sorties MISR |
|---|
| 125 | !Attention A FAIRE |
|---|
| 126 | ! if (cfg%Lmisr_sim) then |
|---|
| 127 | ! print*,'Ecriture sorties Misr' |
|---|
| 128 | ! call write_netcdf4d('clMISR.nc',use_vgrid,nlon,nlat,MISR_N_CTH,7, & |
|---|
| 129 | ! x,y,MISR_CTH,ISCCP_TAU,i,ndays,time,misr%fq_MISR) |
|---|
| 130 | ! endif |
|---|
| 131 | |
|---|
| 132 | ! Sorties ISCCP |
|---|
| 133 | if (cfg%Lisccp_sim) then |
|---|
| 134 | if (cfg%Lclisccp2) then |
|---|
| 135 | do ii=1,7 |
|---|
| 136 | CALL histdef(nid_day_cosp, "clisccp2_"//chcol(ii), & |
|---|
| 137 | "Cloud Fraction as Calculated by the ISCCP Simulator","1", & |
|---|
| 138 | iim,jj_nb,nhori,7,1,7,nvertisccp, 32, & |
|---|
| 139 | "ave(X)", zout,zstoday) |
|---|
| 140 | enddo |
|---|
| 141 | endif |
|---|
| 142 | if (cfg%Lboxtauisccp) then |
|---|
| 143 | CALL histdef(nid_day_cosp, "boxtauisccp", & |
|---|
| 144 | "Optical Depth in Each Column as Calculated by the ISCCP Simulator","1", & |
|---|
| 145 | iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, & |
|---|
| 146 | "ave(X)", zout,zstoday) |
|---|
| 147 | endif |
|---|
| 148 | if (cfg%Lboxptopisccp) then |
|---|
| 149 | CALL histdef(nid_day_cosp, "boxptopisccp", & |
|---|
| 150 | "Cloud Top Pressure in Each Column as Calculated by the ISCCP Simulator","Pa", & |
|---|
| 151 | iim,jj_nb,nhori,Ncolumns,1,Ncolumns,nvertcol, 32, & |
|---|
| 152 | "ave(X)", zout,zstoday) |
|---|
| 153 | endif |
|---|
| 154 | if (cfg%Ltclisccp) then |
|---|
| 155 | CALL histdef(nid_day_cosp, "tclisccp", & |
|---|
| 156 | "Total Cloud Fraction as Calculated by the ISCCP Simulator", "1", & |
|---|
| 157 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 158 | "ave(X)", zout,zstoday) |
|---|
| 159 | endif |
|---|
| 160 | if (cfg%Lctpisccp) then |
|---|
| 161 | CALL histdef(nid_day_cosp, "ctpisccp", & |
|---|
| 162 | "Mean Cloud Top Pressure as Calculated by the ISCCP Simulator", "Pa", & |
|---|
| 163 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 164 | "ave(X)", zout,zstoday) |
|---|
| 165 | endif |
|---|
| 166 | if (cfg%Ltauisccp) then |
|---|
| 167 | CALL histdef(nid_day_cosp, "tauisccp", & |
|---|
| 168 | "Optical Depth as Calculated by the ISCCP Simulator", "1", & |
|---|
| 169 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 170 | "ave(X)", zout,zstoday) |
|---|
| 171 | endif |
|---|
| 172 | if (cfg%Lalbisccp) then |
|---|
| 173 | CALL histdef(nid_day_cosp, "albisccp", & |
|---|
| 174 | "Mean Cloud Albedo as Calculated by the ISCCP Simulator", "1", & |
|---|
| 175 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 176 | "ave(X)", zout,zstoday) |
|---|
| 177 | endif |
|---|
| 178 | if (cfg%Lmeantbisccp) then |
|---|
| 179 | CALL histdef(nid_day_cosp, "meantbisccp", & |
|---|
| 180 | " Mean all-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", & |
|---|
| 181 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 182 | "ave(X)", zout,zstoday) |
|---|
| 183 | endif |
|---|
| 184 | if (cfg%Lmeantbclrisccp) then |
|---|
| 185 | CALL histdef(nid_day_cosp, "meantbclrisccp", & |
|---|
| 186 | "Mean clear-sky 10.5 micron brightness temperature as calculated by the ISCCP Simulator","K", & |
|---|
| 187 | iim, jj_nb,nhori,1,1,1,-99,32, & |
|---|
| 188 | "ave(X)", zout,zstoday) |
|---|
| 189 | endif |
|---|
| 190 | endif ! Isccp |
|---|
| 191 | |
|---|
| 192 | |
|---|
| 193 | CALL histend(nid_day_cosp) |
|---|
| 194 | !$OMP END MASTER |
|---|
| 195 | !$OMP BARRIER |
|---|
