[644] | 1 | c |
---|
| 2 | c $Header$ |
---|
| 3 | c |
---|
[1352] | 4 | cIM on initialise les variables |
---|
[644] | 5 | c |
---|
[1707] | 6 | missing_val=nf90_fill_real |
---|
| 7 | c |
---|
| 8 | cIM freq_moyNMC = frequences auxquelles on moyenne les champs accumules |
---|
| 9 | cIM sur les niveaux de pression standard du NMC |
---|
| 10 | DO n=1, nout |
---|
| 11 | freq_moyNMC(n)=freq_outNMC(n)/freq_calNMC(n) |
---|
| 12 | ENDDO |
---|
| 13 | c |
---|
[1352] | 14 | CALL ini_undefSTD(itap,freq_outNMC) |
---|
[644] | 15 | c |
---|
[1352] | 16 | cIM on interpole les champs sur les niveaux STD de pression |
---|
| 17 | cIM a chaque pas de temps de la physique |
---|
[644] | 18 | c |
---|
[1146] | 19 | c-------------------------------------------------------c |
---|
| 20 | c positionnement de l'argument logique a .false. c |
---|
| 21 | c pour ne pas recalculer deux fois la meme chose ! c |
---|
| 22 | c a cet effet un appel a plevel_new a ete deplace c |
---|
| 23 | c a la fin de la serie d'appels c |
---|
| 24 | c la boucle 'DO k=1, nlevSTD' a ete internalisee c |
---|
| 25 | c dans plevel_new, d'ou la creation de cette routine... c |
---|
| 26 | c-------------------------------------------------------c |
---|
[644] | 27 | c |
---|
[1146] | 28 | CALL plevel_new(klon,klev,nlevSTD,.true.,pplay,rlevSTD, |
---|
| 29 | & t_seri,tlevSTD) |
---|
| 30 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 31 | & u_seri,ulevSTD) |
---|
| 32 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 33 | & v_seri,vlevSTD) |
---|
[644] | 34 | c |
---|
[1146] | 35 | |
---|
[644] | 36 | c |
---|
[1146] | 37 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 38 | & zphi/RG,philevSTD) |
---|
| 39 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 40 | & qx(:,:,ivap),qlevSTD) |
---|
| 41 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 42 | & zx_rh*100.,rhlevSTD) |
---|
[644] | 43 | c |
---|
| 44 | DO l=1, klev |
---|
| 45 | DO i=1, klon |
---|
| 46 | zx_tmp_fi3d(i,l)=u_seri(i,l)*v_seri(i,l) |
---|
| 47 | ENDDO !i |
---|
| 48 | ENDDO !l |
---|
[1146] | 49 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 50 | & zx_tmp_fi3d,uvSTD) |
---|
[644] | 51 | c |
---|
| 52 | DO l=1, klev |
---|
| 53 | DO i=1, klon |
---|
| 54 | zx_tmp_fi3d(i,l)=v_seri(i,l)*q_seri(i,l) |
---|
| 55 | ENDDO !i |
---|
| 56 | ENDDO !l |
---|
[1146] | 57 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 58 | & zx_tmp_fi3d,vqSTD) |
---|
[644] | 59 | c |
---|
| 60 | DO l=1, klev |
---|
| 61 | DO i=1, klon |
---|
| 62 | zx_tmp_fi3d(i,l)=v_seri(i,l)*t_seri(i,l) |
---|
| 63 | ENDDO !i |
---|
| 64 | ENDDO !l |
---|
[1146] | 65 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 66 | & zx_tmp_fi3d,vTSTD) |
---|
[644] | 67 | c |
---|
| 68 | DO l=1, klev |
---|
| 69 | DO i=1, klon |
---|
| 70 | zx_tmp_fi3d(i,l)=omega(i,l)*qx(i,l,ivap) |
---|
| 71 | ENDDO !i |
---|
| 72 | ENDDO !l |
---|
[1146] | 73 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 74 | & zx_tmp_fi3d,wqSTD) |
---|
[644] | 75 | c |
---|
| 76 | DO l=1, klev |
---|
| 77 | DO i=1, klon |
---|
| 78 | zx_tmp_fi3d(i,l)=v_seri(i,l)*zphi(i,l)/RG |
---|
| 79 | ENDDO !i |
---|
| 80 | ENDDO !l |
---|
[1146] | 81 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 82 | & zx_tmp_fi3d,vphiSTD) |
---|
[644] | 83 | c |
---|
| 84 | DO l=1, klev |
---|
| 85 | DO i=1, klon |
---|
| 86 | zx_tmp_fi3d(i,l)=omega(i,l)*t_seri(i,l) |
---|
| 87 | ENDDO !i |
---|
| 88 | ENDDO !l |
---|
[1146] | 89 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 90 | & zx_tmp_fi3d,wTSTD) |
---|
[644] | 91 | c |
---|
| 92 | DO l=1, klev |
---|
| 93 | DO i=1, klon |
---|
| 94 | zx_tmp_fi3d(i,l)=u_seri(i,l)*u_seri(i,l) |
---|
| 95 | ENDDO !i |
---|
| 96 | ENDDO !l |
---|
[1146] | 97 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 98 | & zx_tmp_fi3d,u2STD) |
---|
[644] | 99 | c |
---|
| 100 | DO l=1, klev |
---|
| 101 | DO i=1, klon |
---|
| 102 | zx_tmp_fi3d(i,l)=v_seri(i,l)*v_seri(i,l) |
---|
| 103 | ENDDO !i |
---|
| 104 | ENDDO !l |
---|
[1146] | 105 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 106 | & zx_tmp_fi3d,v2STD) |
---|
[644] | 107 | c |
---|
| 108 | DO l=1, klev |
---|
| 109 | DO i=1, klon |
---|
| 110 | zx_tmp_fi3d(i,l)=t_seri(i,l)*t_seri(i,l) |
---|
| 111 | ENDDO !i |
---|
| 112 | ENDDO !l |
---|
[1146] | 113 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 114 | & zx_tmp_fi3d,T2STD) |
---|
| 115 | |
---|
[1398] | 116 | c |
---|
| 117 | zx_tmp_fi3d(:,:)=wo(:,:,1) * dobson_u * 1e3 / zmasse / rmo3 * rmd |
---|
| 118 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 119 | & zx_tmp_fi3d,O3STD) |
---|
| 120 | c |
---|
| 121 | if (read_climoz == 2) THEN |
---|
| 122 | zx_tmp_fi3d(:,:)=wo(:,:,2) * dobson_u * 1e3 / zmasse / rmo3 * rmd |
---|
| 123 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
---|
| 124 | & zx_tmp_fi3d,O3daySTD) |
---|
| 125 | endif |
---|
| 126 | c |
---|
[1146] | 127 | DO l=1, klev |
---|
| 128 | DO i=1, klon |
---|
| 129 | zx_tmp_fi3d(i,l)=paprs(i,l) |
---|
| 130 | ENDDO !i |
---|
| 131 | ENDDO !l |
---|
| 132 | CALL plevel_new(klon,klev,nlevSTD,.true.,zx_tmp_fi3d,rlevSTD, |
---|
| 133 | & omega,wlevSTD) |
---|
[644] | 134 | c |
---|
[1352] | 135 | cIM on somme les valeurs toutes les freq_calNMC secondes |
---|
[644] | 136 | c |
---|
[1418] | 137 | CALL undefSTD(itap,freq_calNMC, read_climoz) |
---|
[644] | 138 | c |
---|
[1352] | 139 | cIM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes) |
---|
[644] | 140 | c |
---|
[1352] | 141 | CALL moy_undefSTD(itap,freq_outNMC,freq_moyNMC) |
---|
[644] | 142 | c |
---|
[1352] | 143 | CALL plevel(klon,klev,.true.,pplay,50000., |
---|
| 144 | & zphi/RG,geo500) |
---|
| 145 | |
---|
| 146 | cIM on interpole a chaque pas de temps le SWup(clr) et SWdn(clr) a 200 hPa |
---|
[644] | 147 | c |
---|
| 148 | CALL plevel(klon,klevp1,.true.,paprs,20000., |
---|
| 149 | $ swdn0,SWdn200clr) |
---|
| 150 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 151 | $ swdn,SWdn200) |
---|
| 152 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 153 | $ swup0,SWup200clr) |
---|
| 154 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 155 | $ swup,SWup200) |
---|
| 156 | c |
---|
| 157 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 158 | $ lwdn0,LWdn200clr) |
---|
| 159 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 160 | $ lwdn,LWdn200) |
---|
| 161 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 162 | $ lwup0,LWup200clr) |
---|
| 163 | CALL plevel(klon,klevp1,.false.,paprs,20000., |
---|
| 164 | $ lwup,LWup200) |
---|
| 165 | c |
---|
[1707] | 166 | twriteSTD(:,:,1)=tsumSTD(:,:,1) |
---|
| 167 | qwriteSTD(:,:,1)=qsumSTD(:,:,1) |
---|
| 168 | rhwriteSTD(:,:,1)=rhsumSTD(:,:,1) |
---|
| 169 | phiwriteSTD(:,:,1)=phisumSTD(:,:,1) |
---|
| 170 | uwriteSTD(:,:,1)=usumSTD(:,:,1) |
---|
| 171 | vwriteSTD(:,:,1)=vsumSTD(:,:,1) |
---|
| 172 | wwriteSTD(:,:,1)=wsumSTD(:,:,1) |
---|
| 173 | |
---|
| 174 | twriteSTD(:,:,2)=tsumSTD(:,:,2) |
---|
| 175 | qwriteSTD(:,:,2)=qsumSTD(:,:,2) |
---|
| 176 | rhwriteSTD(:,:,2)=rhsumSTD(:,:,2) |
---|
| 177 | phiwriteSTD(:,:,2)=phisumSTD(:,:,2) |
---|
| 178 | uwriteSTD(:,:,2)=usumSTD(:,:,2) |
---|
| 179 | vwriteSTD(:,:,2)=vsumSTD(:,:,2) |
---|
| 180 | wwriteSTD(:,:,2)=wsumSTD(:,:,2) |
---|
| 181 | |
---|
| 182 | twriteSTD(:,:,3)=tlevSTD(:,:) |
---|
| 183 | qwriteSTD(:,:,3)=qlevSTD(:,:) |
---|
| 184 | rhwriteSTD(:,:,3)=rhlevSTD(:,:) |
---|
| 185 | phiwriteSTD(:,:,3)=philevSTD(:,:) |
---|
| 186 | uwriteSTD(:,:,3)=ulevSTD(:,:) |
---|
| 187 | vwriteSTD(:,:,3)=vlevSTD(:,:) |
---|
| 188 | wwriteSTD(:,:,3)=wlevSTD(:,:) |
---|
| 189 | |
---|
| 190 | twriteSTD(:,:,4)=tlevSTD(:,:) |
---|
| 191 | qwriteSTD(:,:,4)=qlevSTD(:,:) |
---|
| 192 | rhwriteSTD(:,:,4)=rhlevSTD(:,:) |
---|
| 193 | phiwriteSTD(:,:,4)=philevSTD(:,:) |
---|
| 194 | uwriteSTD(:,:,4)=ulevSTD(:,:) |
---|
| 195 | vwriteSTD(:,:,4)=vlevSTD(:,:) |
---|
| 196 | wwriteSTD(:,:,4)=wlevSTD(:,:) |
---|
| 197 | c |
---|
| 198 | cIM initialisation 5eme fichier de sortie |
---|
| 199 | twriteSTD(:,:,5)=tlevSTD(:,:) |
---|
| 200 | qwriteSTD(:,:,5)=qlevSTD(:,:) |
---|
| 201 | rhwriteSTD(:,:,5)=rhlevSTD(:,:) |
---|
| 202 | phiwriteSTD(:,:,5)=philevSTD(:,:) |
---|
| 203 | uwriteSTD(:,:,5)=ulevSTD(:,:) |
---|
| 204 | vwriteSTD(:,:,5)=vlevSTD(:,:) |
---|
| 205 | wwriteSTD(:,:,5)=wlevSTD(:,:) |
---|
| 206 | c |
---|
| 207 | cIM initialisation 6eme fichier de sortie |
---|
| 208 | twriteSTD(:,:,6)=tlevSTD(:,:) |
---|
| 209 | qwriteSTD(:,:,6)=qlevSTD(:,:) |
---|
| 210 | rhwriteSTD(:,:,6)=rhlevSTD(:,:) |
---|
| 211 | phiwriteSTD(:,:,6)=philevSTD(:,:) |
---|
| 212 | uwriteSTD(:,:,6)=ulevSTD(:,:) |
---|
| 213 | vwriteSTD(:,:,6)=vlevSTD(:,:) |
---|
| 214 | wwriteSTD(:,:,6)=wlevSTD(:,:) |
---|
| 215 | cIM for NMC files |
---|
| 216 | DO n=1, nlevSTD3 |
---|
| 217 | DO k=1, nlevSTD |
---|
| 218 | if(rlevSTD3(n).EQ.rlevSTD(k)) THEN |
---|
| 219 | twriteSTD3(:,n)=tlevSTD(:,k) |
---|
| 220 | qwriteSTD3(:,n)=qlevSTD(:,k) |
---|
| 221 | rhwriteSTD3(:,n)=rhlevSTD(:,k) |
---|
| 222 | phiwriteSTD3(:,n)=philevSTD(:,k) |
---|
| 223 | uwriteSTD3(:,n)=ulevSTD(:,k) |
---|
| 224 | vwriteSTD3(:,n)=vlevSTD(:,k) |
---|
| 225 | wwriteSTD3(:,n)=wlevSTD(:,k) |
---|
| 226 | endif !rlevSTD3(n).EQ.rlevSTD(k) |
---|
| 227 | ENDDO |
---|
| 228 | ENDDO |
---|
| 229 | c |
---|
| 230 | DO n=1, nlevSTD8 |
---|
| 231 | DO k=1, nlevSTD |
---|
| 232 | if(rlevSTD8(n).EQ.rlevSTD(k)) THEN |
---|
| 233 | tnondefSTD8(:,n)=tnondef(:,k,2) |
---|
| 234 | twriteSTD8(:,n)=tsumSTD(:,k,2) |
---|
| 235 | qwriteSTD8(:,n)=qsumSTD(:,k,2) |
---|
| 236 | rhwriteSTD8(:,n)=rhsumSTD(:,k,2) |
---|
| 237 | phiwriteSTD8(:,n)=phisumSTD(:,k,2) |
---|
| 238 | uwriteSTD8(:,n)=usumSTD(:,k,2) |
---|
| 239 | vwriteSTD8(:,n)=vsumSTD(:,k,2) |
---|
| 240 | wwriteSTD8(:,n)=wsumSTD(:,k,2) |
---|
| 241 | endif !rlevSTD8(n).EQ.rlevSTD(k) |
---|
| 242 | ENDDO |
---|
| 243 | ENDDO |
---|