| 1 | SUBROUTINE aerave_new ( ndata, |
|---|
| 2 | & longdata,epdata,omegdata,gdata, |
|---|
| 3 | & longref,epref,temp,nir,longir |
|---|
| 4 | & ,epir,omegir,gir,qref) |
|---|
| 5 | c ,omegaref ) |
|---|
| 6 | |
|---|
| 7 | |
|---|
| 8 | IMPLICIT NONE |
|---|
| 9 | c....................................................................... |
|---|
| 10 | c |
|---|
| 11 | c R.Fournier 02/1996 |
|---|
| 12 | c (modif F.Forget 02/1996) |
|---|
| 13 | c le spectre est decoupe en "nir" bandes et cette routine calcule |
|---|
| 14 | c les donnees radiatives moyenne sur chaque bande : l'optimisation |
|---|
| 15 | c est faite pour une temperature au sol "temp" et une epaisseur |
|---|
| 16 | c optique de l'atmosphere "epref" a la longueur d'onde "longref" |
|---|
| 17 | c |
|---|
| 18 | c dans la version actuelle, les ponderations sont independantes de |
|---|
| 19 | c l'epaisseur optique : c'est a dire que "omegir", "gir" |
|---|
| 20 | c et "epir/epre" sont independants de "epref". |
|---|
| 21 | c en effet les ponderations sont choisies pour une solution exacte |
|---|
| 22 | c en couche mince et milieu isotherme. |
|---|
| 23 | c |
|---|
| 24 | c entree |
|---|
| 25 | c |
|---|
| 26 | c ndata : taille des champs data |
|---|
| 27 | c longdata,epdata,omegdata,gdata : proprietes radiative de l'aerosol |
|---|
| 28 | c (longdata longueur d'onde en METRES) |
|---|
| 29 | c * longref : longueur d'onde a laquelle l'epaisseur optique |
|---|
| 30 | c est connue |
|---|
| 31 | c * epref : epaisseur optique a longref |
|---|
| 32 | c * temp : temperature choisie pour la ponderation (Planck) |
|---|
| 33 | c * nir : nombre d'intervals dans la discretisation spectrale |
|---|
| 34 | c du GCM |
|---|
| 35 | c * longir : longueurs d'onde definissant ces intervals |
|---|
| 36 | c |
|---|
| 37 | c sortie |
|---|
| 38 | c |
|---|
| 39 | c * epir : epaisseur optique moyenne pour chaque interval |
|---|
| 40 | c * omegir : "scattering albedo" moyen pour chaque interval |
|---|
| 41 | c * gir : "assymetry factor" moyen pour chaque interval |
|---|
| 42 | c * qref : extinction coefficient at reference wavelength |
|---|
| 43 | c * omegaref : single scat. albedo at reference wavelength |
|---|
| 44 | c |
|---|
| 45 | c....................................................................... |
|---|
| 46 | c |
|---|
| 47 | REAL longref |
|---|
| 48 | REAL epref |
|---|
| 49 | REAL temp |
|---|
| 50 | INTEGER nir |
|---|
| 51 | REAL*8 longir(nir+1) |
|---|
| 52 | REAL epir(nir) |
|---|
| 53 | REAL omegir(nir) |
|---|
| 54 | REAL gir(nir) |
|---|
| 55 | c |
|---|
| 56 | c....................................................................... |
|---|
| 57 | c |
|---|
| 58 | INTEGER iir,nirmx |
|---|
| 59 | PARAMETER (nirmx=100) |
|---|
| 60 | INTEGER idata,ndata |
|---|
| 61 | c |
|---|
| 62 | c....................................................................... |
|---|
| 63 | c |
|---|
| 64 | REAL emit |
|---|
| 65 | REAL totalemit(nirmx) |
|---|
| 66 | REAL longdata(ndata),epdata(ndata) |
|---|
| 67 | & ,omegdata(ndata),gdata(ndata) |
|---|
| 68 | REAL qextcorrdata(ndata) |
|---|
| 69 | INTEGER ibande,nbande |
|---|
| 70 | PARAMETER (nbande=1000) |
|---|
| 71 | REAL long,deltalong |
|---|
| 72 | INTEGER ilong |
|---|
| 73 | INTEGER i1,i2 |
|---|
| 74 | REAL c1,c2 |
|---|
| 75 | REAL factep,qextcorr,omeg,g |
|---|
| 76 | REAL qref,omegaref |
|---|
| 77 | c |
|---|
| 78 | c....................................................................... |
|---|
| 79 | c |
|---|
| 80 | DOUBLE PRECISION tmp1 |
|---|
| 81 | REAL tmp2,tmp3 |
|---|
| 82 | c |
|---|
| 83 | c |
|---|
| 84 | long=longref |
|---|
| 85 | |
|---|
| 86 | |
|---|
| 87 | !if(nir.eq.27)then |
|---|
| 88 | !print*,'long',long |
|---|
| 89 | !print*,'longdata',longdata |
|---|
| 90 | !print*,'epdata',epdata |
|---|
| 91 | !print*,'omegdata',omegdata |
|---|
| 92 | !print*,'gdata',gdata |
|---|
| 93 | !print*,'data looks aok!' |
|---|
| 94 | |
|---|
| 95 | !print*,'ndata=',ndata |
|---|
| 96 | !print*,'longdata',shape(longdata) |
|---|
| 97 | !print*,'epdata',shape(epdata) |
|---|
| 98 | !print*,'omegdata',shape(omegdata) |
|---|
| 99 | !print*,'gdata',shape(gdata) |
|---|
| 100 | ! print*,'longref',longref |
|---|
| 101 | !print*,'epref',epref |
|---|
| 102 | !print*,'temp',temp |
|---|
| 103 | !print*,'nir',nir |
|---|
| 104 | !print*,'longir',longir |
|---|
| 105 | !print*,'epir',epir |
|---|
| 106 | !print*,'omegir',gir |
|---|
| 107 | !print*,'qref',qref |
|---|
| 108 | !print*,'longir=',longir |
|---|
| 109 | !stop |
|---|
| 110 | !endif |
|---|
| 111 | |
|---|
| 112 | |
|---|
| 113 | c******************************************************** |
|---|
| 114 | c interpolation |
|---|
| 115 | ilong=1 |
|---|
| 116 | DO idata=2,ndata |
|---|
| 117 | IF (long.gt.longdata(idata)) ilong=idata |
|---|
| 118 | ENDDO |
|---|
| 119 | i1=ilong |
|---|
| 120 | i2=ilong+1 |
|---|
| 121 | IF (i2.gt.ndata) i2=ndata |
|---|
| 122 | IF (long.lt.longdata(1)) i2=1 |
|---|
| 123 | IF (i1.eq.i2) THEN |
|---|
| 124 | c1=1.E+0 |
|---|
| 125 | c2=0.E+0 |
|---|
| 126 | ELSE |
|---|
| 127 | c1=(longdata(i2)-long) / (longdata(i2)-longdata(i1)) |
|---|
| 128 | c2=(longdata(i1)-long) / (longdata(i1)-longdata(i2)) |
|---|
| 129 | ENDIF |
|---|
| 130 | c******************************************************** |
|---|
| 131 | c |
|---|
| 132 | qref=c1*epdata(i1)+c2*epdata(i2) |
|---|
| 133 | omegaref=c1*omegdata(i1)+c2*omegdata(i2) |
|---|
| 134 | factep=qref/epref |
|---|
| 135 | DO idata=1,ndata |
|---|
| 136 | qextcorrdata(idata)=epdata(idata)/factep |
|---|
| 137 | ENDDO |
|---|
| 138 | c |
|---|
| 139 | c....................................................................... |
|---|
| 140 | c |
|---|
| 141 | DO iir=1,nir |
|---|
| 142 | c |
|---|
| 143 | c....................................................................... |
|---|
| 144 | c |
|---|
| 145 | deltalong=(longir(iir+1)-longir(iir)) / nbande |
|---|
| 146 | totalemit(iir)=0.E+0 |
|---|
| 147 | epir(iir)=0.E+0 |
|---|
| 148 | omegir(iir)=0.E+0 |
|---|
| 149 | gir(iir)=0.E+0 |
|---|
| 150 | c |
|---|
| 151 | c....................................................................... |
|---|
| 152 | c |
|---|
| 153 | DO ibande=1,nbande |
|---|
| 154 | c |
|---|
| 155 | c....................................................................... |
|---|
| 156 | c |
|---|
| 157 | long=longir(iir) + (ibande-0.5E+0) * deltalong |
|---|
| 158 | CALL blackl(DBLE(long),DBLE(temp),tmp1) |
|---|
| 159 | emit=REAL(tmp1) |
|---|
| 160 | c |
|---|
| 161 | c....................................................................... |
|---|
| 162 | c |
|---|
| 163 | c******************************************************** |
|---|
| 164 | c interpolation |
|---|
| 165 | ilong=1 |
|---|
| 166 | DO idata=2,ndata |
|---|
| 167 | IF (long.gt.longdata(idata)) ilong=idata |
|---|
| 168 | ENDDO |
|---|
| 169 | i1=ilong |
|---|
| 170 | i2=ilong+1 |
|---|
| 171 | IF (i2.gt.ndata) i2=ndata |
|---|
| 172 | IF (long.lt.longdata(1)) i2=1 |
|---|
| 173 | IF (i1.eq.i2) THEN |
|---|
| 174 | c1=1.E+0 |
|---|
| 175 | c2=0.E+0 |
|---|
| 176 | ELSE |
|---|
| 177 | c1=(longdata(i2)-long) / (longdata(i2)-longdata(i1)) |
|---|
| 178 | c2=(longdata(i1)-long) / (longdata(i1)-longdata(i2)) |
|---|
| 179 | ENDIF |
|---|
| 180 | c******************************************************** |
|---|
| 181 | c |
|---|
| 182 | qextcorr=c1*qextcorrdata(i1)+c2*qextcorrdata(i2) |
|---|
| 183 | omeg=c1*omegdata(i1)+c2*omegdata(i2) |
|---|
| 184 | g=c1*gdata(i1)+c2*gdata(i2) |
|---|
| 185 | c |
|---|
| 186 | c....................................................................... |
|---|
| 187 | c |
|---|
| 188 | totalemit(iir)=totalemit(iir)+deltalong*emit |
|---|
| 189 | epir(iir)=epir(iir)+deltalong*emit*qextcorr |
|---|
| 190 | omegir(iir)=omegir(iir)+deltalong*emit*omeg*qextcorr |
|---|
| 191 | gir(iir)=gir(iir)+deltalong*emit*omeg*qextcorr*g |
|---|
| 192 | c |
|---|
| 193 | c....................................................................... |
|---|
| 194 | c |
|---|
| 195 | ENDDO |
|---|
| 196 | c |
|---|
| 197 | c....................................................................... |
|---|
| 198 | c |
|---|
| 199 | gir(iir)=gir(iir)/omegir(iir) |
|---|
| 200 | omegir(iir)=omegir(iir)/epir(iir) |
|---|
| 201 | epir(iir)=epir(iir)/totalemit(iir) |
|---|
| 202 | c |
|---|
| 203 | c....................................................................... |
|---|
| 204 | c |
|---|
| 205 | ENDDO |
|---|
| 206 | c |
|---|
| 207 | c...................................................................... |
|---|
| 208 | c |
|---|
| 209 | c Diagnostic de controle si on moyenne sur tout le spectre vis ou IR : |
|---|
| 210 | c ~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~ |
|---|
| 211 | c tmp2=0.E+0 |
|---|
| 212 | c DO iir=1,nir |
|---|
| 213 | c tmp2=tmp2+totalemit(iir) |
|---|
| 214 | c ENDDO |
|---|
| 215 | c tmp3=5.67E-8 * temp**4 |
|---|
| 216 | c IF (abs((tmp2-tmp3)/tmp3).gt.0.05E+0) THEN |
|---|
| 217 | c PRINT *,'!!!! <---> il manque du Planck (voir moyenne.F)' |
|---|
| 218 | c PRINT *,'somme des bandes :',tmp2,'--- Planck:',tmp3 |
|---|
| 219 | c ENDIF |
|---|
| 220 | c |
|---|
| 221 | c...................................................................... |
|---|
| 222 | c |
|---|
| 223 | RETURN |
|---|
| 224 | END |
|---|
