| 1 | c |
|---|
| 2 | c $Header$ |
|---|
| 3 | c |
|---|
| 4 | cIM on initialise les variables |
|---|
| 5 | c |
|---|
| 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 |
|---|
| 14 | CALL ini_undefSTD(itap,freq_outNMC) |
|---|
| 15 | c |
|---|
| 16 | cIM on interpole les champs sur les niveaux STD de pression |
|---|
| 17 | cIM a chaque pas de temps de la physique |
|---|
| 18 | c |
|---|
| 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 |
|---|
| 27 | c |
|---|
| 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) |
|---|
| 34 | c |
|---|
| 35 | |
|---|
| 36 | c |
|---|
| 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) |
|---|
| 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 |
|---|
| 49 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 50 | & zx_tmp_fi3d,uvSTD) |
|---|
| 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 |
|---|
| 57 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 58 | & zx_tmp_fi3d,vqSTD) |
|---|
| 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 |
|---|
| 65 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 66 | & zx_tmp_fi3d,vTSTD) |
|---|
| 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 |
|---|
| 73 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 74 | & zx_tmp_fi3d,wqSTD) |
|---|
| 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 |
|---|
| 81 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 82 | & zx_tmp_fi3d,vphiSTD) |
|---|
| 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 |
|---|
| 89 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 90 | & zx_tmp_fi3d,wTSTD) |
|---|
| 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 |
|---|
| 97 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 98 | & zx_tmp_fi3d,u2STD) |
|---|
| 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 |
|---|
| 105 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 106 | & zx_tmp_fi3d,v2STD) |
|---|
| 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 |
|---|
| 113 | CALL plevel_new(klon,klev,nlevSTD,.false.,pplay,rlevSTD, |
|---|
| 114 | & zx_tmp_fi3d,T2STD) |
|---|
| 115 | |
|---|
| 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 |
|---|
| 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) |
|---|
| 134 | c |
|---|
| 135 | cIM on somme les valeurs toutes les freq_calNMC secondes |
|---|
| 136 | c |
|---|
| 137 | CALL undefSTD(itap,freq_calNMC, read_climoz) |
|---|
| 138 | c |
|---|
| 139 | cIM on moyenne a la fin du mois ou du jour (toutes les freq_outNMC secondes) |
|---|
| 140 | c |
|---|
| 141 | CALL moy_undefSTD(itap,freq_outNMC,freq_moyNMC) |
|---|
| 142 | c |
|---|
| 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 |
|---|
| 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 |
|---|
| 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 |
|---|
