| [524] | 1 | ! |
|---|
| 2 | ! $Header$ |
|---|
| 3 | ! |
|---|
| 4 | IF (ok_isccp) THEN |
|---|
| 5 | c |
|---|
| [735] | 6 | c$OMP MASTER |
|---|
| [684] | 7 | ndex2d = 0 |
|---|
| 8 | ndex3d = 0 |
|---|
| 9 | c |
|---|
| [524] | 10 | c pour les champs instantannes, il faut mettre la meme valeur pour |
|---|
| 11 | c zout et zsto. |
|---|
| 12 | c dtime est passe par ailleurs a histbeg |
|---|
| [644] | 13 | c zstophy = frequence de stockage des champs tous les pdt physiques |
|---|
| [524] | 14 | c zout = frequence d'ecriture des champs |
|---|
| [684] | 15 | cIM 300505 zstophy = dtime |
|---|
| 16 | c appel du simulateur toutes les 3heures |
|---|
| 17 | zcals(1) = dtime *6. !toutes les 3h (en s) |
|---|
| 18 | DO n=1, napisccp |
|---|
| 19 | zcalh(n) = zcals(n)/3600. !stoutes les Xh (en heures) |
|---|
| 20 | ENDDO !n |
|---|
| [524] | 21 | c |
|---|
| 22 | c ecriture 8 fois par jour |
|---|
| 23 | c zout = dtime * REAL(NINT(86400./dtime*ecrit_isccp)) |
|---|
| 24 | c ecriture toutes les 2h (12 fois par jour) |
|---|
| 25 | c zout = dtime * 4. |
|---|
| [644] | 26 | c ecriture toutes les 1/2 h (48 fois par jour) |
|---|
| [524] | 27 | c zout = dtime |
|---|
| 28 | c |
|---|
| [684] | 29 | c IF(freqout_isccp.EQ.1.) THEN |
|---|
| 30 | c ecriture jounaliere |
|---|
| 31 | zout_isccp(1) = ecrit_day !(en s) |
|---|
| [524] | 32 | c ecriture mensuelle |
|---|
| [684] | 33 | c zout = dtime * ecrit_mth !(en s) |
|---|
| 34 | DO n=1, napisccp |
|---|
| 35 | zoutj(n)=zout_isccp(n)/86400. !(en jours) |
|---|
| [524] | 36 | c |
|---|
| [684] | 37 | c le nombre de sous-colonnes ncol : ncol=(100.*zcalh)/zoutd |
|---|
| 38 | ncol(n)=NINT((100.*zcalh(n))/zoutj(n)) |
|---|
| 39 | PRINT*,'n ncol',n,ncol(n) |
|---|
| 40 | c |
|---|
| 41 | DO l=1, ncol(n) |
|---|
| 42 | vertlev(l,n)=float(l) |
|---|
| 43 | ENDDO !ncol |
|---|
| 44 | c |
|---|
| 45 | ENDDO !n |
|---|
| 46 | |
|---|
| [524] | 47 | c PRINT*, 'La frequence de sortie ISCCP est de ', ecrit_isccp |
|---|
| 48 | c |
|---|
| 49 | idayref = day_ref |
|---|
| 50 | CALL ymds2ju(annee_ref, 1, idayref, 0.0, zjulian) |
|---|
| [684] | 51 | c write(*,*)'ISCCP ', itau_phy, zjulian |
|---|
| [524] | 52 | c |
|---|
| 53 | c |
|---|
| 54 | c definition coordonnees lon,lat en globale |
|---|
| 55 | c |
|---|
| [704] | 56 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon) |
|---|
| 57 | cym DO i = 1, iim |
|---|
| 58 | cym zx_lon(i,1) = rlon(i+1) |
|---|
| 59 | cym zx_lon(i,jjmp1) = rlon(i+1) |
|---|
| 60 | cym ENDDO |
|---|
| [524] | 61 | |
|---|
| [704] | 62 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlat,zx_lat) |
|---|
| [524] | 63 | c |
|---|
| [684] | 64 | cIM BEG region |
|---|
| [735] | 65 | cym Desole dans un premier temps le mode region ne marchera pas |
|---|
| [704] | 66 | cym Il faudra voir dans un second temps pour l'implementer |
|---|
| [735] | 67 | cym Mais cela posera des problemes au niveau de la reconstruction |
|---|
| [704] | 68 | |
|---|
| [684] | 69 | imin_ins=1 |
|---|
| 70 | imax_ins=iim |
|---|
| 71 | jmin_ins=1 |
|---|
| 72 | jmax_ins=jjmp1 |
|---|
| [718] | 73 | cym do i=1,iim-1 |
|---|
| 74 | cym if(zx_lon(i,1).lt.lonmin_ins) imin_ins=i |
|---|
| 75 | cym if(zx_lon(i,1).le.lonmax_ins) imax_ins=i+1 |
|---|
| 76 | cym enddo |
|---|
| 77 | cym do j=1,jjmp1 |
|---|
| 78 | cym if(zx_lat(1,j).ge.latmin_ins) jmax_ins=j |
|---|
| 79 | cym if(zx_lat(1,j).gt.latmax_ins) jmin_ins=j |
|---|
| 80 | cym enddo |
|---|
| [684] | 81 | c |
|---|
| 82 | print*,'On stoke le fichier histISCCP instantanne sur ', |
|---|
| 83 | s imin_ins,imax_ins,jmin_ins,jmax_ins |
|---|
| [735] | 84 | cym print*,'On stoke le fichier histISCCP instantanne sur ', |
|---|
| 85 | cym s zx_lon(imin_ins,1),zx_lon(imax_ins,1), |
|---|
| 86 | cym s zx_lat(1,jmin_ins),zx_lat(1,jmax_ins) |
|---|
| [684] | 87 | cIM END region |
|---|
| 88 | c |
|---|
| 89 | IF(1.EQ.0) THEN |
|---|
| [704] | 90 | cym CALL histbeg("histISCCP.nc", iim,zx_lon(:,1),jjmp1,zx_lat(1,:), |
|---|
| 91 | cym . 1, iim, 1, jjmp1, |
|---|
| 92 | cym . itau_phy, zjulian, dtime, |
|---|
| 93 | cym . nhori, nid_isccp) |
|---|
| 94 | CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime, |
|---|
| [524] | 95 | . nhori, nid_isccp) |
|---|
| [704] | 96 | |
|---|
| [684] | 97 | ENDIF !(1.EQ.0) THEN |
|---|
| [524] | 98 | c |
|---|
| [704] | 99 | cym CALL histbeg("histISCCP.nc", iim,zx_lon(:,1), |
|---|
| 100 | cym . jjmp1,zx_lat(1,:), |
|---|
| 101 | cym . imin_ins,imax_ins-imin_ins+1, |
|---|
| 102 | cym . jmin_ins,jmax_ins-jmin_ins+1, |
|---|
| 103 | cym . itau_phy, zjulian, dtime, |
|---|
| 104 | cym . nhori, nid_isccp) |
|---|
| 105 | |
|---|
| 106 | CALL histbeg_phy("histISCCP.nc", itau_phy, zjulian, dtime, |
|---|
| [684] | 107 | . nhori, nid_isccp) |
|---|
| [704] | 108 | |
|---|
| [684] | 109 | c |
|---|
| 110 | IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN |
|---|
| [524] | 111 | CALL histvert(nid_isccp, "cldtopres","Cloud Top Pressure","mb", |
|---|
| 112 | . lmaxm1, cldtopres, nvert,'down') |
|---|
| [684] | 113 | ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN |
|---|
| 114 | CALL histvert(nid_isccp,"cldtopres3","Cloud Top Pressure","mb", |
|---|
| 115 | . lmax3, cldtopres3, nvert3,'down') |
|---|
| 116 | ENDIF |
|---|
| 117 | DO n=1, napisccp |
|---|
| 118 | CALL histvert(nid_isccp, "Nbcol"//verticaxe(n), |
|---|
| 119 | . "Nb of Column"//verticaxe(n),"1", |
|---|
| 120 | . ncol(n), vertlev(:,n), nvlev(n),'up') |
|---|
| 121 | ENDDO |
|---|
| [524] | 122 | c |
|---|
| [644] | 123 | IF(type_run.EQ."ENSP".OR.type_run.EQ."CLIM") THEN |
|---|
| 124 | c |
|---|
| [524] | 125 | c variables a ecrire |
|---|
| [684] | 126 | c |
|---|
| 127 | DO n=1, napisccp |
|---|
| [524] | 128 | c |
|---|
| 129 | DO k=1, kmaxm1 |
|---|
| [684] | 130 | CALL histdef(nid_isccp, "cldISCCP_"//taulev(k)//typinout(n), |
|---|
| [524] | 131 | . "LMDZ ISCCP cld", "%", |
|---|
| [704] | 132 | . iim, jjphy_nb,nhori,lmaxm1,1,lmaxm1,nvert,32, |
|---|
| [684] | 133 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| [524] | 134 | ENDDO |
|---|
| 135 | c |
|---|
| [684] | 136 | CALL histdef(nid_isccp, "nsunlit"//typinout(n), |
|---|
| [524] | 137 | . "Nb of calls with sunlit ", "%", |
|---|
| [704] | 138 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 139 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| [524] | 140 | c |
|---|
| [684] | 141 | ENDDO |
|---|
| 142 | c |
|---|
| [644] | 143 | ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN |
|---|
| 144 | c |
|---|
| [684] | 145 | DO n=1, napisccp |
|---|
| [644] | 146 | c |
|---|
| [684] | 147 | c print*,'n=',n,' avant histdef(..LMDZ ISCCP cld' |
|---|
| 148 | c |
|---|
| 149 | DO k=1, kmaxm1 |
|---|
| 150 | DO l=1, lmaxm1 |
|---|
| 151 | c |
|---|
| 152 | CALL histdef(nid_isccp, pclev(l)//taulev(k)//typinout(n), |
|---|
| [644] | 153 | . "LMDZ ISCCP cld "//cnameisccp(l,k), "%", |
|---|
| [704] | 154 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 155 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| [644] | 156 | c |
|---|
| [684] | 157 | ENDDO |
|---|
| [644] | 158 | ENDDO |
|---|
| 159 | c |
|---|
| [684] | 160 | c print*,'n=',n,' avant histdef(..Nb of calls sunlit' |
|---|
| 161 | CALL histdef(nid_isccp, "nsunlit"//typinout(n), |
|---|
| [644] | 162 | . "Nb of calls with sunlit ", "%", |
|---|
| [704] | 163 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 164 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| [644] | 165 | c |
|---|
| [684] | 166 | c 9types de nuages ISCCP-D2 |
|---|
| 167 | CALL histdef(nid_isccp, "cirr", |
|---|
| 168 | . "Cirrus lk ISCCP-D2", "%", |
|---|
| [704] | 169 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 170 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 171 | c |
|---|
| 172 | CALL histdef(nid_isccp, "cist", |
|---|
| 173 | . "CiSt lk ISCCP-D2", "%", |
|---|
| [704] | 174 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 175 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 176 | c |
|---|
| 177 | CALL histdef(nid_isccp, "deep", |
|---|
| 178 | . "Deep lk ISCCP-D2", "%", |
|---|
| [704] | 179 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 180 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 181 | c |
|---|
| 182 | CALL histdef(nid_isccp, "alcu", |
|---|
| 183 | . "AlCu lk ISCCP-D2", "%", |
|---|
| [704] | 184 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 185 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 186 | c |
|---|
| 187 | CALL histdef(nid_isccp, "alst", |
|---|
| 188 | . "AlSt lk ISCCP-D2", "%", |
|---|
| [704] | 189 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 190 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 191 | c |
|---|
| 192 | CALL histdef(nid_isccp, "nist", |
|---|
| 193 | . "NiSt lk ISCCP-D2", "%", |
|---|
| [704] | 194 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 195 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 196 | c |
|---|
| 197 | CALL histdef(nid_isccp, "cumu", |
|---|
| 198 | . "Cumu lk ISCCP-D2", "%", |
|---|
| [704] | 199 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 200 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 201 | c |
|---|
| 202 | CALL histdef(nid_isccp, "stcu", |
|---|
| 203 | . "StCu lk ISCCP-D2", "%", |
|---|
| [704] | 204 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 205 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 206 | c |
|---|
| 207 | CALL histdef(nid_isccp, "stra", |
|---|
| 208 | . "Stra lk ISCCP-D2", "%", |
|---|
| [704] | 209 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 210 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 211 | c |
|---|
| 212 | c 3_epaisseurs_optiques x3_pressions_au_sommet_des_nuages types de nuages |
|---|
| 213 | CALL histdef(nid_isccp, "thin", |
|---|
| 214 | . "Opt. thin ISCCP-D2 like clouds", "%", |
|---|
| [704] | 215 | . iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32, |
|---|
| [684] | 216 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 217 | c |
|---|
| 218 | CALL histdef(nid_isccp, "mid", |
|---|
| 219 | . "Opt. intermediate ISCCP-D2 like clouds", "%", |
|---|
| [704] | 220 | . iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32, |
|---|
| [684] | 221 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 222 | c |
|---|
| 223 | CALL histdef(nid_isccp, "thick", |
|---|
| 224 | . "Opt. thick ISCCP-D2 like clouds", "%", |
|---|
| [704] | 225 | . iim, jjphy_nb,nhori,lmax3,1,lmax3,nvert3,32, |
|---|
| [684] | 226 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 227 | c |
|---|
| 228 | c IF(1.EQ.0) THEN |
|---|
| 229 | c IF(n.EQ.3) THEN |
|---|
| 230 | c IF(n.EQ.1) THEN |
|---|
| 231 | c |
|---|
| 232 | cIM 070905 BEG |
|---|
| 233 | IF(1.EQ.0) THEN |
|---|
| 234 | print*,'n=',n,' avant histdef(..boxptop axe' |
|---|
| 235 | cIM verif boxptop |
|---|
| 236 | CALL histdef(nid_isccp,"boxptop"//verticaxe(n), |
|---|
| 237 | . "Boxptop axe"//verticaxe(n), "mb", |
|---|
| [704] | 238 | . iim, jjphy_nb,nhori, |
|---|
| [684] | 239 | . ncol(n),1,ncol(n),nvlev(n),32, |
|---|
| 240 | cIM . ncolmx,1,ncolmx,nvlev,32, |
|---|
| 241 | cIM . "inst(X)",dtime,dtime) |
|---|
| 242 | . "ave(X)",zcals(n),zout_isccp(n)) |
|---|
| 243 | ENDIF !(1.EQ.0) THEN |
|---|
| 244 | cIM 070905 END |
|---|
| 245 | c ENDIF !(n.EQ.3) THEN |
|---|
| 246 | c ENDIF !(1.EQ.0) THEN |
|---|
| 247 | c |
|---|
| 248 | c print*,'n=',n,' avant histdef(..seed axe' |
|---|
| 249 | CALL histdef(nid_isccp, "seed"//verticaxe(n), |
|---|
| 250 | . "seed axe"//verticaxe(n), "-", |
|---|
| [704] | 251 | . iim, jjphy_nb,nhori,1,1,1,-99,32, |
|---|
| [684] | 252 | cIM . "inst(X)", dtime,dtime) |
|---|
| 253 | . "ave(X)", zcals(n),zout_isccp(n)) |
|---|
| 254 | c |
|---|
| 255 | ENDDO !n |
|---|
| [644] | 256 | ENDIF |
|---|
| [524] | 257 | CALL histend(nid_isccp) |
|---|
| 258 | c |
|---|
| [735] | 259 | c$OMP END MASTER |
|---|
| [524] | 260 | ENDIF ! ok_isccp |
|---|
