[644] | 1 | c |
---|
| 2 | c $Header$ |
---|
| 3 | c |
---|
[524] | 4 | IF (ok_mensuel) THEN |
---|
| 5 | c |
---|
| 6 | ndex2d = 0 |
---|
| 7 | ndex3d = 0 |
---|
| 8 | |
---|
| 9 | itau_w = itau_phy + itap |
---|
| 10 | |
---|
| 11 | c |
---|
[644] | 12 | IF(type_run.EQ."CLIM".OR.type_run.EQ."ENSP") THEN |
---|
| 13 | c |
---|
[524] | 14 | IF(lev_histmth.GE.1) THEN |
---|
| 15 | c |
---|
| 16 | c Champs 2D: |
---|
| 17 | c |
---|
[766] | 18 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) |
---|
| 19 | CALL histwrite_phy(nid_mth,"phis",itau_w,pphis) |
---|
[524] | 20 | C |
---|
[644] | 21 | cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) |
---|
[766] | 22 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d) |
---|
| 23 | CALL histwrite_phy(nid_mth,"aire",itau_w,airephy) |
---|
[524] | 24 | c |
---|
[644] | 25 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100. |
---|
[766] | 26 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 27 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w, |
---|
| 28 | $ zx_tmp_fi2d) |
---|
[524] | 29 | c |
---|
[644] | 30 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter) |
---|
[766] | 31 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 32 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w, |
---|
| 33 | $ zx_tmp_fi2d) |
---|
[644] | 34 | c |
---|
[766] | 35 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d) |
---|
| 36 | CALL histwrite_phy(nid_mth,"slp",itau_w,slp) |
---|
[524] | 37 | |
---|
[766] | 38 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) |
---|
| 39 | CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol) |
---|
[524] | 40 | c |
---|
[766] | 41 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d) |
---|
| 42 | CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m) |
---|
[524] | 43 | c |
---|
| 44 | c ENSEMBLES BEG |
---|
[766] | 45 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 46 | CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m) |
---|
[524] | 47 | c |
---|
[766] | 48 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 49 | CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m) |
---|
[524] | 50 | c |
---|
| 51 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d) |
---|
| 52 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d) |
---|
| 53 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d) |
---|
[766] | 54 | c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d, |
---|
[524] | 55 | c . iim*jjmp1,ndex2d) |
---|
| 56 | c |
---|
| 57 | DO i=1, klon |
---|
| 58 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 59 | ENDDO |
---|
[766] | 60 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 61 | CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d) |
---|
[524] | 62 | c |
---|
| 63 | DO i = 1, klon |
---|
| 64 | zx_tmp_fi2d(i) = pctsrf(i,is_sic) |
---|
| 65 | ENDDO |
---|
[766] | 66 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 67 | CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d) |
---|
[524] | 68 | c |
---|
| 69 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
[766] | 70 | c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d, |
---|
[524] | 71 | c . iim*jjmp1*klev,ndex3d) |
---|
| 72 | c |
---|
| 73 | c ENSEMBLES END |
---|
| 74 | c |
---|
[766] | 75 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d) |
---|
| 76 | CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m) |
---|
[524] | 77 | c |
---|
[766] | 78 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d) |
---|
| 79 | CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m) |
---|
[524] | 80 | C |
---|
[766] | 81 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d) |
---|
| 82 | CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m) |
---|
[524] | 83 | C |
---|
| 84 | DO i = 1, klon |
---|
| 85 | zx_tmp_fi2d(i) = paprs(i,1) |
---|
| 86 | ENDDO |
---|
[766] | 87 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 88 | CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d) |
---|
[524] | 89 | c |
---|
[766] | 90 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d) |
---|
| 91 | CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf) |
---|
[524] | 92 | c |
---|
| 93 | if (.not. ok_veget) then |
---|
[766] | 94 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d) |
---|
| 95 | CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol) |
---|
[524] | 96 | endif |
---|
| 97 | c |
---|
| 98 | DO i = 1, klon |
---|
| 99 | zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) |
---|
| 100 | ENDDO |
---|
[766] | 101 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 102 | CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d) |
---|
[524] | 103 | c |
---|
[766] | 104 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d) |
---|
| 105 | CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain) |
---|
[524] | 106 | c |
---|
| 107 | DO i = 1, klon |
---|
| 108 | zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) |
---|
| 109 | ENDDO |
---|
[766] | 110 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 111 | CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d) |
---|
[524] | 112 | c |
---|
| 113 | DO i = 1, klon |
---|
| 114 | zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) |
---|
| 115 | ENDDO |
---|
[766] | 116 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 117 | CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d) |
---|
[524] | 118 | c |
---|
[766] | 119 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) |
---|
| 120 | CALL histwrite_phy(nid_mth,"snow",itau_w,snow_fall) |
---|
[524] | 121 | c |
---|
[644] | 122 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d) |
---|
| 123 | c CALL histwrite(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, |
---|
| 124 | c . ndex2d) |
---|
[524] | 125 | c |
---|
[766] | 126 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) |
---|
| 127 | CALL histwrite_phy(nid_mth,"evap",itau_w,evap) |
---|
[524] | 128 | c |
---|
[766] | 129 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) |
---|
| 130 | CALL histwrite_phy(nid_mth,"tops",itau_w,topsw) |
---|
[524] | 131 | c |
---|
[766] | 132 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) |
---|
| 133 | CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0) |
---|
[524] | 134 | c |
---|
[766] | 135 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) |
---|
| 136 | CALL histwrite_phy(nid_mth,"topl",itau_w,toplw) |
---|
[524] | 137 | c |
---|
[766] | 138 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) |
---|
| 139 | CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0) |
---|
[524] | 140 | c |
---|
| 141 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) |
---|
[766] | 142 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 143 | CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d) |
---|
[524] | 144 | c |
---|
[766] | 145 | c |
---|
[524] | 146 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) |
---|
[766] | 147 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 148 | CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d) |
---|
[524] | 149 | c |
---|
| 150 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) |
---|
[766] | 151 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 152 | CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d) |
---|
[524] | 153 | c |
---|
| 154 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) |
---|
[766] | 155 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 156 | CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d) |
---|
[524] | 157 | c |
---|
[766] | 158 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d) |
---|
| 159 | CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200) |
---|
[524] | 160 | c |
---|
[766] | 161 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d) |
---|
| 162 | CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr) |
---|
[524] | 163 | c |
---|
[766] | 164 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d) |
---|
| 165 | CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200) |
---|
[524] | 166 | c |
---|
[766] | 167 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d) |
---|
| 168 | CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr) |
---|
[524] | 169 | c |
---|
[766] | 170 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d) |
---|
| 171 | CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200) |
---|
[524] | 172 | c |
---|
[766] | 173 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d) |
---|
| 174 | CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr) |
---|
[524] | 175 | c |
---|
| 176 | zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon) |
---|
[766] | 177 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 178 | CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d) |
---|
[524] | 179 | c |
---|
| 180 | zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon) |
---|
[766] | 181 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 182 | CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d) |
---|
[524] | 183 | c |
---|
| 184 | c200 END |
---|
[766] | 185 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) |
---|
| 186 | CALL histwrite_phy(nid_mth,"sols",itau_w,solsw) |
---|
[524] | 187 | c |
---|
[766] | 188 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) |
---|
| 189 | CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0) |
---|
[524] | 190 | c |
---|
[766] | 191 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) |
---|
| 192 | CALL histwrite_phy(nid_mth,"soll",itau_w,sollw) |
---|
[524] | 193 | c |
---|
[766] | 194 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d) |
---|
| 195 | CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol) |
---|
[684] | 196 | c |
---|
[766] | 197 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) |
---|
| 198 | CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0) |
---|
[524] | 199 | c |
---|
| 200 | c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1) |
---|
| 201 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) |
---|
[766] | 202 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 203 | CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d) |
---|
[524] | 204 | c |
---|
| 205 | c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1) |
---|
| 206 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) |
---|
[766] | 207 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 208 | CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
[524] | 209 | c |
---|
| 210 | c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1) |
---|
| 211 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) |
---|
[766] | 212 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 213 | CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d) |
---|
[524] | 214 | c |
---|
| 215 | c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1) |
---|
| 216 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) |
---|
[766] | 217 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 218 | CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d) |
---|
[524] | 219 | c |
---|
[644] | 220 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon) |
---|
| 221 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon) |
---|
| 222 | zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) |
---|
[766] | 223 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 224 | CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d) |
---|
[524] | 225 | c |
---|
[644] | 226 | DO i=1, klon |
---|
| 227 | sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 228 | ENDDO !i=1, klon |
---|
| 229 | |
---|
| 230 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon) |
---|
| 231 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon) |
---|
| 232 | zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) |
---|
[766] | 233 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 234 | CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
[524] | 235 | c |
---|
[766] | 236 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) |
---|
| 237 | CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown) |
---|
[524] | 238 | c |
---|
[644] | 239 | c DO i=1, klon |
---|
| 240 | cIM 080904 sollwdownclr(i)=lwdn0(i,1) |
---|
| 241 | c sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 242 | c ENDDO !i=1, klon |
---|
[524] | 243 | c |
---|
[766] | 244 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d) |
---|
| 245 | CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr) |
---|
[524] | 246 | c |
---|
[896] | 247 | c |
---|
| 248 | c effets des aerosols |
---|
| 249 | c |
---|
| 250 | IF (ok_ade.OR.ok_aie) THEN |
---|
| 251 | zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon) |
---|
[925] | 252 | CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d) |
---|
[896] | 253 | c |
---|
| 254 | zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon) |
---|
[925] | 255 | CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d) |
---|
[896] | 256 | c |
---|
| 257 | zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon) |
---|
[925] | 258 | CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d) |
---|
[896] | 259 | c |
---|
| 260 | zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon) |
---|
[925] | 261 | CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d) |
---|
[896] | 262 | endif |
---|
| 263 | |
---|
| 264 | c |
---|
[766] | 265 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) |
---|
| 266 | CALL histwrite_phy(nid_mth,"bils",itau_w,bils) |
---|
[524] | 267 | c |
---|
| 268 | zx_tmp_fi2d(1:klon)=-1*sens(1:klon) |
---|
| 269 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) |
---|
[766] | 270 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 271 | CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d) |
---|
[524] | 272 | c |
---|
[766] | 273 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) |
---|
| 274 | CALL histwrite_phy(nid_mth,"fder",itau_w,fder) |
---|
[524] | 275 | c |
---|
| 276 | c |
---|
| 277 | c DO i = 1, klon |
---|
| 278 | c zx_tmp_fi2d(i) = fluxu(i,1) |
---|
| 279 | c ENDDO |
---|
| 280 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 281 | c CALL histwrite(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 282 | c |
---|
| 283 | c DO i = 1, klon |
---|
| 284 | c zx_tmp_fi2d(i) = fluxv(i,1) |
---|
| 285 | c ENDDO |
---|
| 286 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 287 | c CALL histwrite(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
| 288 | c |
---|
| 289 | c |
---|
[766] | 290 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d) |
---|
| 291 | CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte) |
---|
[524] | 292 | c |
---|
[766] | 293 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d) |
---|
| 294 | CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving) |
---|
| 295 | |
---|
| 296 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqfonte,zx_tmp_2d) |
---|
| 297 | CALL histwrite_phy(nid_mth,"fqfonte",itau_w,zxfqfonte) |
---|
| 298 | |
---|
[524] | 299 | DO nsrf = 1, nbsrf |
---|
| 300 | zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) |
---|
[766] | 301 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 302 | CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w, |
---|
| 303 | $ zx_tmp_fi2d) |
---|
[524] | 304 | C |
---|
| 305 | zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) |
---|
[766] | 306 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 307 | CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w, |
---|
| 308 | $ zx_tmp_fi2d) |
---|
[524] | 309 | ENDDO |
---|
| 310 | c |
---|
| 311 | DO nsrf = 1, nbsrf |
---|
| 312 | CYYY |
---|
[892] | 313 | IF(nsrf.GE.2) THEN |
---|
[644] | 314 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100. |
---|
[766] | 315 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 316 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(nsrf),itau_w, |
---|
| 317 | & zx_tmp_fi2d) |
---|
[644] | 318 | c |
---|
| 319 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf) |
---|
[766] | 320 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 321 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(nsrf),itau_w, |
---|
| 322 | $ zx_tmp_fi2d) |
---|
[897] | 323 | ENDIF !nsrf.GE.2 |
---|
[524] | 324 | C |
---|
| 325 | zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) |
---|
[766] | 326 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 327 | CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w, |
---|
| 328 | $ zx_tmp_fi2d) |
---|
[524] | 329 | C |
---|
| 330 | zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) |
---|
[766] | 331 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 332 | CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w, |
---|
| 333 | $ zx_tmp_fi2d) |
---|
[524] | 334 | C |
---|
| 335 | zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) |
---|
[766] | 336 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 337 | CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w, |
---|
| 338 | $ zx_tmp_fi2d) |
---|
[524] | 339 | C |
---|
| 340 | zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) |
---|
[766] | 341 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 342 | CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w, |
---|
| 343 | $ zx_tmp_fi2d) |
---|
[524] | 344 | c |
---|
| 345 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) |
---|
[766] | 346 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 347 | CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w, |
---|
| 348 | $ zx_tmp_fi2d) |
---|
[524] | 349 | C |
---|
| 350 | zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) |
---|
[766] | 351 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 352 | CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w, |
---|
| 353 | $ zx_tmp_fi2d) |
---|
[524] | 354 | C |
---|
[684] | 355 | zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) |
---|
[766] | 356 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 357 | CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w, |
---|
| 358 | $ zx_tmp_fi2d) |
---|
[684] | 359 | C |
---|
[644] | 360 | IF(1.EQ.0) THEN |
---|
| 361 | C |
---|
| 362 | zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf) |
---|
[766] | 363 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 364 | CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w, |
---|
| 365 | $ zx_tmp_fi2d) |
---|
[644] | 366 | C |
---|
| 367 | zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon)) |
---|
| 368 | $ *pctsrf(1:klon,nsrf) |
---|
[766] | 369 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 370 | CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w, |
---|
| 371 | $ zx_tmp_fi2d) |
---|
[644] | 372 | C |
---|
| 373 | ENDIF !(1.EQ.0) THEN |
---|
| 374 | C |
---|
[524] | 375 | END DO |
---|
| 376 | c |
---|
[766] | 377 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) |
---|
| 378 | CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm) |
---|
[524] | 379 | c |
---|
[766] | 380 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) |
---|
| 381 | CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh) |
---|
[524] | 382 | c |
---|
[766] | 383 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) |
---|
| 384 | CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl) |
---|
[524] | 385 | c |
---|
[766] | 386 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) |
---|
| 387 | CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm) |
---|
[524] | 388 | c |
---|
[766] | 389 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) |
---|
| 390 | CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh) |
---|
[524] | 391 | c |
---|
[766] | 392 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) |
---|
| 393 | CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt) |
---|
[524] | 394 | c |
---|
[766] | 395 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) |
---|
| 396 | CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq) |
---|
[524] | 397 | c |
---|
| 398 | zx_tmp_fi2d(1:klon) = flwp(1:klon) |
---|
[766] | 399 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 400 | CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d) |
---|
[524] | 401 | c |
---|
| 402 | zx_tmp_fi2d(1:klon) = fiwp(1:klon) |
---|
[766] | 403 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 404 | CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d) |
---|
[524] | 405 | c |
---|
[766] | 406 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) |
---|
| 407 | CALL histwrite_phy(nid_mth,"ue",itau_w,ue) |
---|
[524] | 408 | c |
---|
[766] | 409 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) |
---|
| 410 | CALL histwrite_phy(nid_mth,"ve",itau_w,ve) |
---|
[524] | 411 | c |
---|
[766] | 412 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) |
---|
| 413 | CALL histwrite_phy(nid_mth,"uq",itau_w,uq) |
---|
[524] | 414 | c |
---|
[766] | 415 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) |
---|
| 416 | CALL histwrite_phy(nid_mth,"vq",itau_w,vq) |
---|
[524] | 417 | cKE43 |
---|
| 418 | IF (iflag_con.GE.3) THEN ! sb |
---|
| 419 | c |
---|
[766] | 420 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d) |
---|
| 421 | CALL histwrite_phy(nid_mth,"cape",itau_w,cape) |
---|
[524] | 422 | c |
---|
[766] | 423 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d) |
---|
| 424 | CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase) |
---|
[524] | 425 | c |
---|
[766] | 426 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d) |
---|
| 427 | CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct) |
---|
[524] | 428 | c |
---|
[766] | 429 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d) |
---|
| 430 | CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf) |
---|
[524] | 431 | c |
---|
| 432 | c34EK |
---|
| 433 | c |
---|
[766] | 434 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d) |
---|
| 435 | CALL histwrite_phy(nid_mth,"prw",itau_w,prw) |
---|
[524] | 436 | c |
---|
[644] | 437 | ENDIF !iflag_con.GE.3 |
---|
[524] | 438 | c |
---|
| 439 | c |
---|
[644] | 440 | cIM cf. AM 081204 BEG |
---|
| 441 | cHBTM2 |
---|
[524] | 442 | c |
---|
[644] | 443 | c |
---|
[766] | 444 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d) |
---|
| 445 | CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh) |
---|
[644] | 446 | c |
---|
[766] | 447 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d) |
---|
| 448 | CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt) |
---|
[644] | 449 | c |
---|
[766] | 450 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d) |
---|
| 451 | CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl) |
---|
[644] | 452 | c |
---|
[766] | 453 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d) |
---|
| 454 | CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL) |
---|
[644] | 455 | c |
---|
[766] | 456 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d) |
---|
| 457 | CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL) |
---|
[644] | 458 | c |
---|
[766] | 459 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d) |
---|
| 460 | CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL) |
---|
[644] | 461 | c |
---|
[766] | 462 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d) |
---|
| 463 | CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm) |
---|
[644] | 464 | c |
---|
[766] | 465 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d) |
---|
| 466 | CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1) |
---|
[644] | 467 | c |
---|
[766] | 468 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d) |
---|
| 469 | CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2) |
---|
[644] | 470 | c |
---|
[766] | 471 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d) |
---|
| 472 | CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3) |
---|
[644] | 473 | c |
---|
| 474 | cIM cf. AM 081204 END |
---|
| 475 | c |
---|
[524] | 476 | c Champs interpolles sur des niveaux de pression |
---|
[644] | 477 | DO k=1, nlevSTD |
---|
[524] | 478 | c |
---|
[644] | 479 | c bb=clevSTD(k) |
---|
| 480 | c |
---|
| 481 | c IF(k.GE.2) THEN |
---|
| 482 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 483 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 484 | c aa=clevSTD(k) |
---|
| 485 | c bb=aa(1:lnblnk1(aa)) |
---|
| 486 | c ENDIF |
---|
| 487 | c |
---|
| 488 | IF(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
| 489 | $ bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
| 490 | c |
---|
[766] | 491 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d) |
---|
| 492 | CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2)) |
---|
[524] | 493 | c |
---|
[766] | 494 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d) |
---|
| 495 | CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2)) |
---|
[524] | 496 | c |
---|
[766] | 497 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d) |
---|
| 498 | CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2)) |
---|
[524] | 499 | c |
---|
[766] | 500 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d) |
---|
| 501 | CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2)) |
---|
[524] | 502 | c |
---|
[644] | 503 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR. |
---|
| 504 | ENDDO !nlevSTD |
---|
| 505 | c |
---|
| 506 | DO i=1, klon |
---|
| 507 | IF (pctsrf(i,is_oce).GT.epsfra) THEN |
---|
| 508 | zx_tmp_fi2d(i) = fluxo(i) |
---|
| 509 | ELSE |
---|
[684] | 510 | zx_tmp_fi2d(i) = 0. |
---|
[644] | 511 | ENDIF |
---|
| 512 | ENDDO |
---|
| 513 | c |
---|
[766] | 514 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 515 | CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d) |
---|
[644] | 516 | c |
---|
| 517 | DO i=1, klon |
---|
| 518 | IF (pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 519 | zx_tmp_fi2d(i) = fluxg(i) |
---|
| 520 | ELSE |
---|
[684] | 521 | zx_tmp_fi2d(i) = 0. |
---|
[644] | 522 | ENDIF |
---|
| 523 | ENDDO |
---|
| 524 | c |
---|
[766] | 525 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 526 | CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d) |
---|
[644] | 527 | c |
---|
[684] | 528 | DO i=1, klon |
---|
| 529 | IF (pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 530 | . pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 531 | zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ |
---|
| 532 | . ftsol(i, is_sic) * pctsrf(i,is_sic))/ |
---|
| 533 | . (pctsrf(i,is_oce)+pctsrf(i,is_sic)) |
---|
| 534 | ELSE |
---|
| 535 | zx_tmp_fi2d(i) = 273.15 |
---|
| 536 | ENDIF |
---|
| 537 | ENDDO |
---|
| 538 | c |
---|
[766] | 539 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 540 | CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d) |
---|
[684] | 541 | c |
---|
[644] | 542 | IF (OCEAN.EQ.'force ') THEN |
---|
| 543 | c |
---|
| 544 | DO i=1, klon |
---|
| 545 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 546 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 547 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 548 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 549 | ELSE |
---|
| 550 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 551 | ENDIF |
---|
| 552 | ENDDO |
---|
[766] | 553 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 554 | CALL histwrite_phy(nid_mth,"lmt_bils",itau_w,zx_tmp_fi2d) |
---|
[644] | 555 | c |
---|
| 556 | ELSE IF (OCEAN.EQ.'slab ') THEN |
---|
| 557 | c |
---|
| 558 | DO i=1, klon |
---|
| 559 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 560 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 561 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 562 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 563 | ELSE |
---|
| 564 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 565 | ENDIF |
---|
| 566 | ENDDO |
---|
[766] | 567 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 568 | CALL histwrite_phy(nid_mth,"slab_bils",itau_w,zx_tmp_fi2d) |
---|
[644] | 569 | c |
---|
| 570 | DO i=1, klon |
---|
| 571 | IF(pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 572 | $ pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 573 | zx_tmp_fi2d(i)=tslab(i) |
---|
| 574 | ELSE |
---|
| 575 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 576 | ENDIF |
---|
| 577 | ENDDO !i=1, klon |
---|
[766] | 578 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 579 | CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d) |
---|
[644] | 580 | c |
---|
[766] | 581 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d) |
---|
| 582 | CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice) |
---|
[644] | 583 | c |
---|
[766] | 584 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d) |
---|
| 585 | CALL histwrite_phy(nid_mth,"siceh",itau_w, seaice/1000.) |
---|
[644] | 586 | c |
---|
| 587 | ENDIF !OCEAN.EQ.force/slab |
---|
| 588 | c |
---|
[524] | 589 | ENDIF !lev_histmth.GE.1 |
---|
| 590 | c |
---|
| 591 | IF(lev_histmth.GE.2) THEN |
---|
| 592 | c |
---|
| 593 | c Champs 3D: |
---|
| 594 | C |
---|
| 595 | c |
---|
[766] | 596 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d) |
---|
| 597 | CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc) |
---|
[524] | 598 | c |
---|
[766] | 599 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d) |
---|
| 600 | CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc) |
---|
[524] | 601 | c |
---|
[766] | 602 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 603 | CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri) |
---|
[524] | 604 | c |
---|
| 605 | c#ifdef histmthNMC |
---|
| 606 | c#include "write_histmthNMC.h" |
---|
| 607 | c#endif |
---|
| 608 | c |
---|
[766] | 609 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) |
---|
| 610 | CALL histwrite_phy(nid_mth,"ovap",itau_w,qx(:,:,ivap)) |
---|
[524] | 611 | c |
---|
| 612 | c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/ |
---|
| 613 | c . (1-qx(1:klon,1:klev,ivap)) |
---|
| 614 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 615 | c CALL histwrite(nid_mth,"wvap",itau_w,zx_tmp_3d, |
---|
| 616 | c . iim*jjmp1*klev,ndex3d) |
---|
| 617 | c |
---|
[766] | 618 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) |
---|
| 619 | CALL histwrite_phy(nid_mth,"geop",itau_w,zphi) |
---|
[524] | 620 | c |
---|
[766] | 621 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) |
---|
| 622 | CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri) |
---|
[524] | 623 | c |
---|
[766] | 624 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) |
---|
| 625 | CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri) |
---|
[524] | 626 | c |
---|
[766] | 627 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) |
---|
| 628 | CALL histwrite_phy(nid_mth,"vitw",itau_w,omega) |
---|
[524] | 629 | c |
---|
[766] | 630 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) |
---|
| 631 | CALL histwrite_phy(nid_mth,"pres",itau_w,pplay) |
---|
[524] | 632 | c |
---|
[766] | 633 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra, zx_tmp_3d) |
---|
| 634 | CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra) |
---|
[524] | 635 | c |
---|
[766] | 636 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d) |
---|
| 637 | CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon) |
---|
[524] | 638 | c |
---|
[766] | 639 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) |
---|
| 640 | CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh) |
---|
| 641 | |
---|
[644] | 642 | DO k=1, klev |
---|
| 643 | DO i=1, klon |
---|
| 644 | cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968 |
---|
| 645 | c $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 646 | c $ *(paprs(i,1)/101325.0) |
---|
| 647 | c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968 |
---|
| 648 | zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968 |
---|
| 649 | $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 650 | $ *(paprs(i,1)/101325.0) |
---|
| 651 | ENDDO !i |
---|
| 652 | ENDDO !k |
---|
| 653 | c |
---|
| 654 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) |
---|
| 655 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d) |
---|
[766] | 656 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 657 | CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d) |
---|
[524] | 658 | c |
---|
[766] | 659 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) |
---|
| 660 | CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t) |
---|
[524] | 661 | c |
---|
[766] | 662 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap), |
---|
| 663 | cym . zx_tmp_3d) |
---|
| 664 | CALL histwrite_phy(nid_mth,"dqphy",itau_w, d_qx(:,:,ivap)) |
---|
| 665 | cIM 22.03.04 END |
---|
[524] | 666 | c |
---|
| 667 | ENDIF !lev_histmth.GE.2 |
---|
| 668 | c |
---|
| 669 | IF(lev_histmth.GE.3) THEN |
---|
| 670 | c |
---|
| 671 | DO nsrf=1, nbsrf |
---|
| 672 | c |
---|
| 673 | zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) |
---|
[766] | 674 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 675 | CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w, |
---|
| 676 | $ zx_tmp_fi2d) |
---|
[524] | 677 | C |
---|
| 678 | zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) |
---|
[766] | 679 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 680 | CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w, |
---|
| 681 | $ zx_tmp_fi2d) |
---|
[524] | 682 | c |
---|
| 683 | zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) |
---|
[766] | 684 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) |
---|
| 685 | CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w |
---|
| 686 | $ ,zx_tmp_fi2d) |
---|
[524] | 687 | c |
---|
| 688 | ENDDO !nsrf=1, nbsrf |
---|
| 689 | c |
---|
[766] | 690 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) |
---|
| 691 | CALL histwrite_phy(nid_mth,"albs",itau_w,albsol) |
---|
[524] | 692 | c |
---|
[766] | 693 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) |
---|
| 694 | CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw) |
---|
[524] | 695 | c |
---|
| 696 | ENDIF !lev_histmth.GE.3 |
---|
| 697 | c |
---|
| 698 | c |
---|
| 699 | IF(lev_histmth.GE.4) THEN |
---|
| 700 | c |
---|
[766] | 701 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d) |
---|
| 702 | CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0) |
---|
| 703 | |
---|
[653] | 704 | IF (iflag_con.GE.3) THEN |
---|
| 705 | c |
---|
[766] | 706 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d) |
---|
| 707 | CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma) |
---|
[524] | 708 | c |
---|
[766] | 709 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d) |
---|
| 710 | CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd) |
---|
[684] | 711 | c |
---|
[766] | 712 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d) |
---|
| 713 | CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd) |
---|
[524] | 714 | c |
---|
[766] | 715 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d) |
---|
| 716 | CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0) |
---|
[524] | 717 | c |
---|
[653] | 718 | ENDIF !iflag_con.GE.3 |
---|
[524] | 719 | c |
---|
[766] | 720 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) |
---|
| 721 | CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn) |
---|
[524] | 722 | c |
---|
[766] | 723 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) |
---|
| 724 | CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn) |
---|
[524] | 725 | c |
---|
| 726 | cIM: 101003 : K/30min ==> K/s |
---|
| 727 | zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys |
---|
[766] | 728 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 729 | CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d) |
---|
[524] | 730 | c |
---|
[644] | 731 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
[766] | 732 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 733 | CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d) |
---|
[524] | 734 | c |
---|
| 735 | cIM: 101003 : K/30min ==> K/s |
---|
| 736 | zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys |
---|
[766] | 737 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 738 | CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d) |
---|
| 739 | cIM: 071003 |
---|
[524] | 740 | cIM: 101003 : K/30min ==> K/s |
---|
| 741 | zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ |
---|
| 742 | . d_t_eva(1:klon,1:klev))/pdtphys |
---|
[766] | 743 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 744 | CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d) |
---|
[524] | 745 | c |
---|
[644] | 746 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys |
---|
[766] | 747 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 748 | CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d) |
---|
[524] | 749 | c |
---|
| 750 | cIM: 101003 : K/30min ==> K/s |
---|
| 751 | zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 752 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 753 | CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d) |
---|
[524] | 754 | c |
---|
[644] | 755 | zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 756 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 757 | CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d) |
---|
[524] | 758 | c |
---|
| 759 | cIM: 101003 : K/30min ==> K/s |
---|
| 760 | zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys |
---|
[766] | 761 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 762 | CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d) |
---|
[524] | 763 | c |
---|
[644] | 764 | zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys |
---|
[766] | 765 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 766 | CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d) |
---|
[524] | 767 | c |
---|
| 768 | zpt_conv = 0. |
---|
| 769 | where (ptconv) zpt_conv = 1. |
---|
[766] | 770 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) |
---|
| 771 | CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv) |
---|
[524] | 772 | c |
---|
[766] | 773 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d) |
---|
| 774 | CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs) |
---|
[524] | 775 | c |
---|
| 776 | cIM: 101003 : K/30min ==> K/s |
---|
| 777 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys |
---|
[766] | 778 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 779 | CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d) |
---|
[524] | 780 | c |
---|
[644] | 781 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys |
---|
[766] | 782 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 783 | CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d) |
---|
[524] | 784 | c |
---|
| 785 | cIM: 101003 : K/day ==> K/s |
---|
| 786 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY |
---|
[766] | 787 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 788 | CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d) |
---|
[524] | 789 | c |
---|
| 790 | cIM: 101003 : K/day ==> K/s |
---|
| 791 | zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY |
---|
[766] | 792 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 793 | CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d) |
---|
[524] | 794 | c |
---|
| 795 | cIM: 101003 : K/day ==> K/s |
---|
| 796 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY |
---|
[766] | 797 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 798 | CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d) |
---|
[524] | 799 | c |
---|
| 800 | cIM: 101003 : K/day ==> K/s |
---|
| 801 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY |
---|
[766] | 802 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 803 | CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d) |
---|
[524] | 804 | c |
---|
| 805 | cIM: 101003 : deja en K/s |
---|
| 806 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) |
---|
[766] | 807 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 808 | CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d) |
---|
[524] | 809 | c |
---|
[644] | 810 | zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 811 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 812 | CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d) |
---|
[524] | 813 | c |
---|
[644] | 814 | zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 815 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 816 | CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d) |
---|
[524] | 817 | c |
---|
| 818 | IF (ok_orodr) THEN |
---|
[644] | 819 | zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys |
---|
[766] | 820 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 821 | CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d) |
---|
[524] | 822 | c |
---|
[644] | 823 | zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys |
---|
[766] | 824 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 825 | CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d) |
---|
[524] | 826 | c |
---|
| 827 | ENDIF |
---|
| 828 | C |
---|
| 829 | IF (ok_orolf) THEN |
---|
[644] | 830 | zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys |
---|
[766] | 831 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 832 | CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d) |
---|
[524] | 833 | c |
---|
[644] | 834 | zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys |
---|
[766] | 835 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 836 | CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d) |
---|
[524] | 837 | ENDIF |
---|
| 838 | C |
---|
[862] | 839 | |
---|
| 840 | #if defined(INCA_AER) && defined(CPP_COUPLE) |
---|
| 841 | CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau) |
---|
| 842 | CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1)) |
---|
| 843 | CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2)) |
---|
| 844 | CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi) |
---|
| 845 | CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi) |
---|
| 846 | CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1)) |
---|
| 847 | CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2)) |
---|
| 848 | CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3)) |
---|
| 849 | CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:)) |
---|
| 850 | CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:)) |
---|
| 851 | CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1)) |
---|
| 852 | CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2)) |
---|
| 853 | CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:)) |
---|
| 854 | CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3)) |
---|
| 855 | #endif |
---|
| 856 | #ifndef INCA |
---|
[644] | 857 | IF (nqmax.GE.3) THEN |
---|
[653] | 858 | DO iq=3,nqmax |
---|
[766] | 859 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d) |
---|
| 860 | CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq)) |
---|
[644] | 861 | ENDDO |
---|
| 862 | ENDIF |
---|
[862] | 863 | #endif |
---|
[644] | 864 | ENDIF !lev_histmth.GE.4 |
---|
[524] | 865 | c |
---|
[644] | 866 | ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN |
---|
| 867 | c |
---|
| 868 | IF(lev_histmth.GE.1) THEN |
---|
| 869 | c |
---|
| 870 | c Champs 2D: |
---|
| 871 | c |
---|
[766] | 872 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis/RG,zx_tmp_2d) |
---|
| 873 | CALL histwrite_phy(nid_mth,"phis",itau_w,pphis/RG) |
---|
[644] | 874 | C |
---|
| 875 | cIM 310804 CALL gr_fi_ecrit(1,klon,iim,jjmp1,paire,zx_tmp_2d) |
---|
[766] | 876 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d) |
---|
| 877 | CALL histwrite_phy(nid_mth,"aire",itau_w,airephy) |
---|
[644] | 878 | c |
---|
| 879 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter)* 100. |
---|
[766] | 880 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 881 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_ter),itau_w, |
---|
| 882 | $ zx_tmp_fi2d) |
---|
[644] | 883 | c |
---|
| 884 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_ter) |
---|
[766] | 885 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 886 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_ter),itau_w, |
---|
| 887 | $ zx_tmp_fi2d) |
---|
[644] | 888 | c |
---|
| 889 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic)* 100. |
---|
[766] | 890 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 891 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_lic),itau_w, |
---|
| 892 | $ zx_tmp_fi2d) |
---|
[644] | 893 | c |
---|
| 894 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_lic) |
---|
[766] | 895 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d, zx_tmp_2d) |
---|
| 896 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_lic),itau_w, |
---|
| 897 | $ zx_tmp_fi2d) |
---|
[644] | 898 | c |
---|
| 899 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic)* 100. |
---|
[766] | 900 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 901 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_sic),itau_w, |
---|
| 902 | $ zx_tmp_fi2d) |
---|
[644] | 903 | c |
---|
| 904 | zx_tmp_fi2d(1 : klon) = pctsrf(1 : klon, is_sic) |
---|
[766] | 905 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 906 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_sic),itau_w, |
---|
| 907 | $ zx_tmp_fi2d) |
---|
[644] | 908 | c |
---|
[766] | 909 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfluxlat,zx_tmp_2d) |
---|
| 910 | CALL histwrite_phy(nid_mth,"flat",itau_w,zxfluxlat) |
---|
[644] | 911 | c |
---|
| 912 | zx_tmp_fi2d(1 : klon) = topsw(1:klon)-toplw(1:klon) |
---|
[766] | 913 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 914 | CALL histwrite_phy(nid_mth,"bilTOA",itau_w,zx_tmp_fi2d) |
---|
[644] | 915 | c |
---|
| 916 | c IF(itap.EQ.ecrit_mth) THEN |
---|
| 917 | c DO i=1, klon |
---|
| 918 | c IF(nday_clr(i).GT.0.) THEN |
---|
| 919 | c tops01(i)=tops01(i)/nday_clr(i) |
---|
| 920 | c ELSE |
---|
| 921 | c tops01(i)=1.e+20 |
---|
| 922 | c ENDIF |
---|
| 923 | c ENDDO |
---|
| 924 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, tops01,zx_tmp_2d) |
---|
[766] | 925 | c CALL histwrite_phy(nid_mth,"tops01",itau_w,zx_tmp_2d, |
---|
[644] | 926 | c . iim*jjmp1,ndex2d) |
---|
| 927 | c ENDIF !itap.EQ.ecrit_mth |
---|
| 928 | c |
---|
| 929 | c220404 BEG |
---|
| 930 | c |
---|
| 931 | c temperature tendency due to total diabatic heating |
---|
| 932 | DO l=1, klev |
---|
| 933 | DO i=1, klon |
---|
| 934 | d_t(i,l) = (t_seri(i,l)-t(i,l))/dtime |
---|
| 935 | ENDDO !i |
---|
| 936 | ENDDO !l |
---|
| 937 | c |
---|
[766] | 938 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t,zx_tmp_3d) |
---|
| 939 | CALL histwrite_phy(nid_mth,"dtdia",itau_w,d_t) |
---|
[644] | 940 | c |
---|
| 941 | c temperature tendency due to SW radiation |
---|
| 942 | DO l=1, klev |
---|
| 943 | DO i=1, klon |
---|
| 944 | zx_tmp_fi3d(i,l)=heat(i,l)/RDAY |
---|
| 945 | ENDDO !i |
---|
| 946 | ENDDO !l |
---|
| 947 | c |
---|
[766] | 948 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 949 | CALL histwrite_phy(nid_mth,"dtsw",itau_w,zx_tmp_fi3d) |
---|
[644] | 950 | c |
---|
| 951 | c temperature tendency due to LW radiation |
---|
| 952 | DO l=1, klev |
---|
| 953 | DO i=1, klon |
---|
| 954 | zx_tmp_fi3d(i,l)=-1*cool(i,l)/RDAY |
---|
| 955 | ENDDO !i |
---|
| 956 | ENDDO !l |
---|
| 957 | c |
---|
[766] | 958 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 959 | CALL histwrite_phy(nid_mth,"dtlw",itau_w,zx_tmp_fi3d) |
---|
[644] | 960 | c |
---|
| 961 | c temperature tendency due to moist convective processes |
---|
| 962 | DO l=1, klev |
---|
| 963 | DO i=1, klon |
---|
| 964 | zx_tmp_fi3d(i,l)=d_t_con(i,l)/pdtphys |
---|
| 965 | ENDDO !i |
---|
| 966 | ENDDO !l |
---|
| 967 | c |
---|
[766] | 968 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 969 | CALL histwrite_phy(nid_mth,"dtcon",itau_w,zx_tmp_fi3d) |
---|
[644] | 970 | c |
---|
| 971 | c temperature tendency due to dry convective processes |
---|
| 972 | DO l=1, klev |
---|
| 973 | DO i=1, klon |
---|
| 974 | zx_tmp_fi3d(i,l)=d_t_ajs(i,l)/pdtphys |
---|
| 975 | ENDDO !i |
---|
| 976 | ENDDO !l |
---|
| 977 | c |
---|
[766] | 978 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 979 | CALL histwrite_phy(nid_mth,"dtajs",itau_w,zx_tmp_fi3d) |
---|
[644] | 980 | c |
---|
| 981 | c temperature tendency due to large scale precipitation |
---|
| 982 | DO l=1, klev |
---|
| 983 | DO i=1, klon |
---|
| 984 | zx_tmp_fi3d(i,l)=d_t_lsc(i,l)/pdtphys |
---|
| 985 | ENDDO !i |
---|
| 986 | ENDDO !l |
---|
| 987 | c |
---|
[766] | 988 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 989 | CALL histwrite_phy(nid_mth,"dtlsc",itau_w,zx_tmp_fi3d) |
---|
[644] | 990 | c |
---|
| 991 | c total moisture tendency due to diabatic processes |
---|
| 992 | DO l=1, klev |
---|
| 993 | DO i=1, klon |
---|
| 994 | d_qx(i,l,ivap) = ( q_seri(i,l) - qx(i,l,ivap) ) / dtime |
---|
| 995 | zx_tmp_fi3d(i,l)=d_qx(i,l,ivap) |
---|
| 996 | ENDDO !i |
---|
| 997 | ENDDO !l |
---|
| 998 | c |
---|
[766] | 999 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1, zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1000 | CALL histwrite_phy(nid_mth,"dqdia",itau_w,zx_tmp_fi3d) |
---|
[644] | 1001 | c |
---|
| 1002 | cBEG 210404 BEG |
---|
| 1003 | c |
---|
| 1004 | IF(ok_orodr.AND.ok_orolf) THEN |
---|
| 1005 | DO k = 1, klev |
---|
| 1006 | DO i = 1, klon |
---|
| 1007 | d_u_oli(i,k) = (d_u_oro(i,k) + d_u_lif(i,k))/pdtphys |
---|
| 1008 | d_v_oli(i,k) = (d_v_oro(i,k) + d_v_lif(i,k))/pdtphys |
---|
| 1009 | ENDDO |
---|
| 1010 | ENDDO |
---|
| 1011 | ENDIF !(ok_orodr.AND.ok_orolf) THEN |
---|
| 1012 | c |
---|
[766] | 1013 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u_oli,zx_tmp_3d) |
---|
| 1014 | CALL histwrite_phy(nid_mth,"duoli",itau_w,d_u_oli) |
---|
[644] | 1015 | c |
---|
[766] | 1016 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v_oli,zx_tmp_3d) |
---|
| 1017 | CALL histwrite_phy(nid_mth,"dvoli",itau_w,d_v_oli) |
---|
[644] | 1018 | c |
---|
| 1019 | zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys |
---|
[766] | 1020 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1021 | CALL histwrite_phy(nid_mth,"ducon",itau_w,zx_tmp_fi3d) |
---|
[644] | 1022 | c |
---|
| 1023 | zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys |
---|
[766] | 1024 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1025 | CALL histwrite_phy(nid_mth,"dvcon",itau_w,zx_tmp_fi3d) |
---|
[644] | 1026 | c |
---|
[766] | 1027 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_u,zx_tmp_3d) |
---|
| 1028 | CALL histwrite_phy(nid_mth,"dutot",itau_w,d_u) |
---|
[644] | 1029 | c |
---|
[766] | 1030 | cym CALL gr_fi_ecrit(klev, klon,iim,jjmp1,d_v,zx_tmp_3d) |
---|
| 1031 | CALL histwrite_phy(nid_mth,"dvtot",itau_w,d_v) |
---|
[644] | 1032 | |
---|
| 1033 | cEND 210404 END |
---|
| 1034 | c |
---|
| 1035 | c cldtrue |
---|
| 1036 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldtrue,zx_tmp_3d) |
---|
[766] | 1037 | c CALL histwrite_phy(nid_mth,"cldtrue",itau_w,zx_tmp_3d, |
---|
[644] | 1038 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1039 | c |
---|
| 1040 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldTOA,zx_tmp_3d) |
---|
[766] | 1041 | c CALL histwrite_phy(nid_mth,"cldTOA",itau_w,zx_tmp_3d, |
---|
[644] | 1042 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1043 | c |
---|
| 1044 | c CALL gr_fi_ecrit(klev, klon,iim,jjmp1,cldSFC,zx_tmp_3d) |
---|
[766] | 1045 | c CALL histwrite_phy(nid_mth,"cldSFC",itau_w,zx_tmp_3d, |
---|
[644] | 1046 | c $ iim*jjmp1*klev,ndex3d) |
---|
| 1047 | c |
---|
| 1048 | DO i=1, klon |
---|
| 1049 | cIM 151004 zx_tmp_fi2d(i)=GWDtaux1(i,1)-GWDtaux0(i,1) |
---|
| 1050 | zx_tmp_fi2d(i)=zustrdr(i)+zustrli(i) |
---|
| 1051 | ENDDO |
---|
[766] | 1052 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1053 | CALL histwrite_phy(nid_mth,"GWDtaux",itau_w,zx_tmp_fi2d) |
---|
[644] | 1054 | c |
---|
| 1055 | DO i=1, klon |
---|
| 1056 | cIM 151004 zx_tmp_fi2d(i)=GWDtauy1(i,1)-GWDtauy0(i,1) |
---|
| 1057 | zx_tmp_fi2d(i)=zvstrdr(i)+zvstrli(i) |
---|
| 1058 | ENDDO |
---|
[766] | 1059 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1060 | CALL histwrite_phy(nid_mth,"GWDtauy",itau_w,zx_tmp_fi2d) |
---|
[644] | 1061 | c |
---|
| 1062 | c220404 END |
---|
| 1063 | cAMIP2 END |
---|
| 1064 | c |
---|
[766] | 1065 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, slp,zx_tmp_2d) |
---|
| 1066 | CALL histwrite_phy(nid_mth,"slp",itau_w,slp) |
---|
[644] | 1067 | |
---|
[766] | 1068 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) |
---|
| 1069 | CALL histwrite_phy(nid_mth,"tsol",itau_w,zxtsol) |
---|
[644] | 1070 | c |
---|
[766] | 1071 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zt2m,zx_tmp_2d) |
---|
| 1072 | CALL histwrite_phy(nid_mth,"t2m",itau_w,zt2m) |
---|
[644] | 1073 | c |
---|
| 1074 | c ENSEMBLES BEG |
---|
[766] | 1075 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 1076 | CALL histwrite_phy(nid_mth,"t2m_min",itau_w,zt2m) |
---|
[644] | 1077 | c |
---|
[766] | 1078 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zt2m,zx_tmp_2d) |
---|
| 1079 | CALL histwrite_phy(nid_mth,"t2m_max",itau_w,zt2m) |
---|
[644] | 1080 | c |
---|
| 1081 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsoil(:,1,is_ter),zx_tmp_2d) |
---|
| 1082 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ftsol(:,is_ter),zx_tmp_2d) |
---|
| 1083 | c CALL gr_fi_ecrit(1,klon,iim,jjmp1,ztsol,zx_tmp_2d) |
---|
[766] | 1084 | c CALL histwrite_phy(nid_mth,"tsoil",itau_w,zx_tmp_2d, |
---|
[644] | 1085 | c . iim*jjmp1,ndex2d) |
---|
| 1086 | c |
---|
| 1087 | DO i=1, klon |
---|
| 1088 | zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i)) |
---|
| 1089 | ENDDO |
---|
[766] | 1090 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1091 | CALL histwrite_phy(nid_mth,"wind10m",itau_w,zx_tmp_fi2d) |
---|
[644] | 1092 | c |
---|
| 1093 | DO i = 1, klon |
---|
| 1094 | zx_tmp_fi2d(i) = pctsrf(i,is_sic) |
---|
| 1095 | ENDDO |
---|
[766] | 1096 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1097 | CALL histwrite_phy(nid_mth,"sicf",itau_w,zx_tmp_fi2d) |
---|
[644] | 1098 | c |
---|
| 1099 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
[766] | 1100 | c CALL histwrite_phy(nid_mth,"temp",itau_w,zx_tmp_3d, |
---|
[644] | 1101 | c . iim*jjmp1*klev,ndex3d) |
---|
| 1102 | c |
---|
| 1103 | c ENSEMBLES END |
---|
| 1104 | c |
---|
[766] | 1105 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zq2m,zx_tmp_2d) |
---|
| 1106 | CALL histwrite_phy(nid_mth,"q2m",itau_w,zq2m) |
---|
[644] | 1107 | c |
---|
[766] | 1108 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zu10m,zx_tmp_2d) |
---|
| 1109 | CALL histwrite_phy(nid_mth,"u10m",itau_w,zu10m) |
---|
[644] | 1110 | C |
---|
[766] | 1111 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zv10m,zx_tmp_2d) |
---|
| 1112 | CALL histwrite_phy(nid_mth,"v10m",itau_w,zv10m) |
---|
[644] | 1113 | C |
---|
| 1114 | DO i = 1, klon |
---|
| 1115 | zx_tmp_fi2d(i) = paprs(i,1) |
---|
| 1116 | ENDDO |
---|
[766] | 1117 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1118 | CALL histwrite_phy(nid_mth,"psol",itau_w,zx_tmp_fi2d) |
---|
[644] | 1119 | c |
---|
[766] | 1120 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxqsurf,zx_tmp_2d) |
---|
| 1121 | CALL histwrite_phy(nid_mth,"qsurf",itau_w,zxqsurf) |
---|
[644] | 1122 | c |
---|
| 1123 | if (.not. ok_veget) then |
---|
[766] | 1124 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, qsol,zx_tmp_2d) |
---|
| 1125 | CALL histwrite_phy(nid_mth,"qsol",itau_w,qsol) |
---|
[644] | 1126 | endif |
---|
| 1127 | c |
---|
| 1128 | DO i = 1, klon |
---|
| 1129 | zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i) |
---|
| 1130 | ENDDO |
---|
[766] | 1131 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1132 | CALL histwrite_phy(nid_mth,"precip",itau_w,zx_tmp_fi2d) |
---|
[644] | 1133 | c |
---|
[766] | 1134 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, nday_rain,zx_tmp_2d) |
---|
| 1135 | CALL histwrite_phy(nid_mth,"ndayrain",itau_w,nday_rain) |
---|
[644] | 1136 | c |
---|
| 1137 | DO i = 1, klon |
---|
| 1138 | zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i) |
---|
| 1139 | ENDDO |
---|
[766] | 1140 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1141 | CALL histwrite_phy(nid_mth,"plul",itau_w,zx_tmp_fi2d) |
---|
[644] | 1142 | c |
---|
| 1143 | DO i = 1, klon |
---|
| 1144 | zx_tmp_fi2d(i) = rain_con(i) + snow_con(i) |
---|
| 1145 | ENDDO |
---|
[766] | 1146 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1147 | CALL histwrite_phy(nid_mth,"pluc",itau_w,zx_tmp_fi2d) |
---|
[644] | 1148 | c |
---|
[766] | 1149 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) |
---|
| 1150 | CALL histwrite_phy(nid_mth,"snowf",itau_w,snow_fall) |
---|
[644] | 1151 | c |
---|
| 1152 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zsnow_mass,zx_tmp_2d) |
---|
[766] | 1153 | c CALL histwrite_phy(nid_mth,"snow_mass",itau_w,zx_tmp_2d,iim*jjmp1, |
---|
[644] | 1154 | c . ndex2d) |
---|
| 1155 | c |
---|
[766] | 1156 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) |
---|
| 1157 | CALL histwrite_phy(nid_mth,"evap",itau_w,evap) |
---|
[644] | 1158 | c |
---|
[766] | 1159 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) |
---|
| 1160 | CALL histwrite_phy(nid_mth,"tops",itau_w,topsw) |
---|
[644] | 1161 | c |
---|
[766] | 1162 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw0,zx_tmp_2d) |
---|
| 1163 | CALL histwrite_phy(nid_mth,"tops0",itau_w,topsw0) |
---|
[644] | 1164 | c |
---|
[766] | 1165 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) |
---|
| 1166 | CALL histwrite_phy(nid_mth,"topl",itau_w,toplw) |
---|
[644] | 1167 | c |
---|
[766] | 1168 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw0,zx_tmp_2d) |
---|
| 1169 | CALL histwrite_phy(nid_mth,"topl0",itau_w,toplw0) |
---|
[644] | 1170 | c |
---|
| 1171 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, klevp1 ) |
---|
[766] | 1172 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1173 | CALL histwrite_phy(nid_mth, "SWupTOA",itau_w,zx_tmp_fi2d) |
---|
[644] | 1174 | c |
---|
| 1175 | c |
---|
| 1176 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, klevp1 ) |
---|
[766] | 1177 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1178 | CALL histwrite_phy(nid_mth, "SWupTOAclr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1179 | c |
---|
| 1180 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, klevp1 ) |
---|
[766] | 1181 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1182 | CALL histwrite_phy(nid_mth, "SWdnTOA",itau_w,zx_tmp_fi2d) |
---|
[644] | 1183 | c |
---|
| 1184 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, klevp1 ) |
---|
[766] | 1185 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1186 | CALL histwrite_phy(nid_mth, "SWdnTOAclr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1187 | c |
---|
[766] | 1188 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200,zx_tmp_2d) |
---|
| 1189 | CALL histwrite_phy(nid_mth,"SWup200",itau_w,SWup200) |
---|
[644] | 1190 | c |
---|
[766] | 1191 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWup200clr,zx_tmp_2d) |
---|
| 1192 | CALL histwrite_phy(nid_mth,"SWup200clr",itau_w,SWup200clr) |
---|
[644] | 1193 | c |
---|
[766] | 1194 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200,zx_tmp_2d) |
---|
| 1195 | CALL histwrite_phy(nid_mth,"SWdn200",itau_w,SWdn200) |
---|
[644] | 1196 | c |
---|
[766] | 1197 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, SWdn200clr,zx_tmp_2d) |
---|
| 1198 | CALL histwrite_phy(nid_mth,"SWdn200clr",itau_w,SWdn200clr) |
---|
[644] | 1199 | c |
---|
[766] | 1200 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200,zx_tmp_2d) |
---|
| 1201 | CALL histwrite_phy(nid_mth,"LWup200",itau_w,LWup200) |
---|
[644] | 1202 | c |
---|
[766] | 1203 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, LWup200clr,zx_tmp_2d) |
---|
| 1204 | CALL histwrite_phy(nid_mth,"LWup200clr",itau_w,LWup200clr) |
---|
[644] | 1205 | c |
---|
| 1206 | zx_tmp_fi2d(1:klon)=-1*LWdn200(1:klon) |
---|
[766] | 1207 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1208 | CALL histwrite_phy(nid_mth,"LWdn200",itau_w,zx_tmp_fi2d) |
---|
[644] | 1209 | c |
---|
| 1210 | zx_tmp_fi2d(1:klon)=-1*LWdn200clr(1:klon) |
---|
[766] | 1211 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1212 | CALL histwrite_phy(nid_mth,"LWdn200clr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1213 | c |
---|
| 1214 | c200 END |
---|
[766] | 1215 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) |
---|
| 1216 | CALL histwrite_phy(nid_mth,"sols",itau_w,solsw) |
---|
[644] | 1217 | c |
---|
[766] | 1218 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw0,zx_tmp_2d) |
---|
| 1219 | CALL histwrite_phy(nid_mth,"sols0",itau_w,solsw0) |
---|
[644] | 1220 | c |
---|
[766] | 1221 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) |
---|
| 1222 | CALL histwrite_phy(nid_mth,"soll",itau_w,sollw) |
---|
[644] | 1223 | c |
---|
[766] | 1224 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, radsol,zx_tmp_2d) |
---|
| 1225 | CALL histwrite_phy(nid_mth,"radsol",itau_w,radsol) |
---|
[684] | 1226 | c |
---|
[766] | 1227 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw0,zx_tmp_2d) |
---|
| 1228 | CALL histwrite_phy(nid_mth,"soll0",itau_w,sollw0) |
---|
[644] | 1229 | c |
---|
| 1230 | c zx_tmp_fi2d(1 : klon) = ZFSUP( 1 : klon, 1) |
---|
| 1231 | zx_tmp_fi2d(1 : klon) = swup ( 1 : klon, 1 ) |
---|
[766] | 1232 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1233 | CALL histwrite_phy(nid_mth, "SWupSFC",itau_w,zx_tmp_fi2d) |
---|
[644] | 1234 | c |
---|
| 1235 | c zx_tmp_fi2d(1 : klon) = ZFSUP0( 1 : klon, 1) |
---|
| 1236 | zx_tmp_fi2d(1 : klon) = swup0 ( 1 : klon, 1 ) |
---|
[766] | 1237 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1238 | CALL histwrite_phy(nid_mth, "SWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1239 | c |
---|
| 1240 | c zx_tmp_fi2d(1 : klon) = ZFSDN( 1 : klon, 1) |
---|
| 1241 | zx_tmp_fi2d(1 : klon) = swdn ( 1 : klon, 1 ) |
---|
[766] | 1242 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1243 | CALL histwrite_phy(nid_mth, "SWdnSFC",itau_w,zx_tmp_fi2d) |
---|
[644] | 1244 | c |
---|
| 1245 | c zx_tmp_fi2d(1 : klon) = ZFSDN0( 1 : klon, 1) |
---|
| 1246 | zx_tmp_fi2d(1 : klon) = swdn0 ( 1 : klon, 1 ) |
---|
[766] | 1247 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1248 | CALL histwrite_phy(nid_mth, "SWdnSFCclr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1249 | c |
---|
| 1250 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw(1:klon)+sollwdown(1:klon) |
---|
| 1251 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw(1:klon)-sollwdown(1:klon) |
---|
| 1252 | zx_tmp_fi2d(1:klon)=sollwdown(1:klon)-sollw(1:klon) |
---|
[766] | 1253 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 1254 | CALL histwrite_phy(nid_mth,"LWupSFC",itau_w,zx_tmp_fi2d) |
---|
[644] | 1255 | c |
---|
| 1256 | DO i=1, klon |
---|
| 1257 | sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 1258 | ENDDO !i=1, klon |
---|
| 1259 | c |
---|
| 1260 | cIM 080904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)+sollwdownclr(1:klon) |
---|
| 1261 | cIM 130904 zx_tmp_fi2d(1:klon)=sollw0(1:klon)-sollwdownclr(1:klon) |
---|
| 1262 | zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon) |
---|
[766] | 1263 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d, zx_tmp_2d) |
---|
| 1264 | CALL histwrite_phy(nid_mth,"LWupSFCclr",itau_w,zx_tmp_fi2d) |
---|
[644] | 1265 | c |
---|
[766] | 1266 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) |
---|
| 1267 | CALL histwrite_phy(nid_mth,"LWdnSFC",itau_w,sollwdown) |
---|
[644] | 1268 | c |
---|
| 1269 | c DO i=1, klon |
---|
| 1270 | cIM 08004 sollwdownclr(i)=lwdn0(i,1) |
---|
| 1271 | c sollwdownclr(i)=-1.*lwdn0(i,1) |
---|
| 1272 | c ENDDO !i=1, klon |
---|
| 1273 | c |
---|
[766] | 1274 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdownclr,zx_tmp_2d) |
---|
| 1275 | CALL histwrite_phy(nid_mth,"LWdnSFCclr",itau_w,sollwdownclr) |
---|
[644] | 1276 | c |
---|
[766] | 1277 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) |
---|
| 1278 | CALL histwrite_phy(nid_mth,"bils",itau_w,bils) |
---|
[644] | 1279 | c |
---|
| 1280 | zx_tmp_fi2d(1:klon)=-1*sens(1:klon) |
---|
| 1281 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) |
---|
[766] | 1282 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1283 | CALL histwrite_phy(nid_mth,"sens",itau_w,zx_tmp_fi2d) |
---|
[644] | 1284 | c |
---|
[766] | 1285 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) |
---|
| 1286 | CALL histwrite_phy(nid_mth,"fder",itau_w,fder) |
---|
[644] | 1287 | c |
---|
[897] | 1288 | c Effets des aerosols |
---|
[644] | 1289 | c |
---|
[897] | 1290 | IF (ok_ade.OR.ok_aie) THEN |
---|
| 1291 | zx_tmp_fi2d(1:klon) = topswai(1:klon) - topswad(1:klon) |
---|
| 1292 | CALL histwrite_phy(nid_mth,"topsad",itau_w,zx_tmp_fi2d) |
---|
| 1293 | c |
---|
| 1294 | zx_tmp_fi2d(1:klon) = solswai(1:klon) - solswad(1:klon) |
---|
| 1295 | CALL histwrite_phy(nid_mth,"solsad",itau_w,zx_tmp_fi2d) |
---|
| 1296 | c |
---|
| 1297 | zx_tmp_fi2d(1:klon) = topsw(1:klon) - topswai(1:klon) |
---|
| 1298 | CALL histwrite_phy(nid_mth,"topsai",itau_w,zx_tmp_fi2d) |
---|
| 1299 | c |
---|
| 1300 | zx_tmp_fi2d(1:klon) = solsw(1:klon) - solswai(1:klon) |
---|
| 1301 | CALL histwrite_phy(nid_mth,"solsai",itau_w,zx_tmp_fi2d) |
---|
| 1302 | ENDIF |
---|
| 1303 | c |
---|
| 1304 | c |
---|
[644] | 1305 | c DO i = 1, klon |
---|
| 1306 | c zx_tmp_fi2d(i) = fluxu(i,1) |
---|
| 1307 | c ENDDO |
---|
| 1308 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
[766] | 1309 | c CALL histwrite_phy(nid_mth,"frtu",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
[644] | 1310 | c |
---|
| 1311 | c DO i = 1, klon |
---|
| 1312 | c zx_tmp_fi2d(i) = fluxv(i,1) |
---|
| 1313 | c ENDDO |
---|
| 1314 | c CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
[766] | 1315 | c CALL histwrite_phy(nid_mth,"frtv",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) |
---|
[644] | 1316 | c |
---|
| 1317 | c |
---|
[766] | 1318 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxffonte,zx_tmp_2d) |
---|
| 1319 | CALL histwrite_phy(nid_mth,"ffonte",itau_w,zxffonte) |
---|
[644] | 1320 | c |
---|
[766] | 1321 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxfqcalving,zx_tmp_2d) |
---|
| 1322 | CALL histwrite_phy(nid_mth,"fqcalving",itau_w,zxfqcalving) |
---|
[644] | 1323 | c |
---|
| 1324 | DO nsrf = 1, nbsrf |
---|
| 1325 | zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) |
---|
[766] | 1326 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1327 | CALL histwrite_phy(nid_mth,"taux_"//clnsurf(nsrf),itau_w, |
---|
| 1328 | $ zx_tmp_fi2d) |
---|
[644] | 1329 | C |
---|
| 1330 | zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) |
---|
[766] | 1331 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1332 | CALL histwrite_phy(nid_mth,"tauy_"//clnsurf(nsrf),itau_w, |
---|
| 1333 | $ zx_tmp_fi2d) |
---|
[644] | 1334 | ENDDO |
---|
| 1335 | c |
---|
| 1336 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce)*100. |
---|
[766] | 1337 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1338 | CALL histwrite_phy(nid_mth,"pourc_"//clnsurf(is_oce),itau_w, |
---|
| 1339 | $ zx_tmp_fi2d) |
---|
[644] | 1340 | c |
---|
| 1341 | zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, is_oce) |
---|
[766] | 1342 | cym CALL gr_fi_ecrit(1,klon,iim,jjmp1,zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1343 | CALL histwrite_phy(nid_mth,"fract_"//clnsurf(is_oce),itau_w, |
---|
| 1344 | $ zx_tmp_fi2d) |
---|
[644] | 1345 | c |
---|
| 1346 | DO nsrf = 1, nbsrf |
---|
| 1347 | CYYY |
---|
| 1348 | C |
---|
| 1349 | zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) |
---|
[766] | 1350 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1351 | CALL histwrite_phy(nid_mth,"tsol_"//clnsurf(nsrf),itau_w, |
---|
| 1352 | $ zx_tmp_fi2d) |
---|
[644] | 1353 | C |
---|
| 1354 | zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) |
---|
[766] | 1355 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1356 | CALL histwrite_phy(nid_mth,"sens_"//clnsurf(nsrf),itau_w, |
---|
| 1357 | $ zx_tmp_fi2d) |
---|
[644] | 1358 | C |
---|
| 1359 | zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) |
---|
[766] | 1360 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1361 | CALL histwrite_phy(nid_mth,"lat_"//clnsurf(nsrf),itau_w, |
---|
| 1362 | $ zx_tmp_fi2d) |
---|
[644] | 1363 | C |
---|
| 1364 | zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf) |
---|
[766] | 1365 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1366 | CALL histwrite_phy(nid_mth,"flw_"//clnsurf(nsrf),itau_w, |
---|
| 1367 | $ zx_tmp_fi2d) |
---|
[644] | 1368 | c |
---|
| 1369 | zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf) |
---|
[766] | 1370 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1371 | CALL histwrite_phy(nid_mth,"fsw_"//clnsurf(nsrf),itau_w, |
---|
| 1372 | $ zx_tmp_fi2d) |
---|
[644] | 1373 | C |
---|
| 1374 | zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf) |
---|
[766] | 1375 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1376 | CALL histwrite_phy(nid_mth,"wbils_"//clnsurf(nsrf),itau_w, |
---|
| 1377 | $ zx_tmp_fi2d) |
---|
[644] | 1378 | C |
---|
[684] | 1379 | zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf) |
---|
[766] | 1380 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1381 | CALL histwrite_phy(nid_mth,"wbilo_"//clnsurf(nsrf),itau_w, |
---|
| 1382 | $ zx_tmp_fi2d) |
---|
[684] | 1383 | C |
---|
[644] | 1384 | IF(1.EQ.0) THEN |
---|
| 1385 | C |
---|
| 1386 | zx_tmp_fi2d(1:klon) = fevap(1:klon,nsrf)*pctsrf(1:klon,nsrf) |
---|
[766] | 1387 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1388 | CALL histwrite_phy(nid_mth,"evapour_"//clnsurf(nsrf),itau_w, |
---|
| 1389 | $ zx_tmp_fi2d) |
---|
[644] | 1390 | C |
---|
| 1391 | zx_tmp_fi2d(1:klon) = (rain_fall(1:klon) + snow_fall(1:klon)) |
---|
| 1392 | $ *pctsrf(1:klon,nsrf) |
---|
[766] | 1393 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1394 | CALL histwrite_phy(nid_mth,"prepour_"//clnsurf(nsrf),itau_w, |
---|
| 1395 | $ zx_tmp_fi2d) |
---|
[644] | 1396 | C |
---|
| 1397 | ENDIF !(1.EQ.0) THEN |
---|
| 1398 | C |
---|
| 1399 | END DO |
---|
| 1400 | c |
---|
[766] | 1401 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) |
---|
| 1402 | CALL histwrite_phy(nid_mth,"cdrm",itau_w,cdragm) |
---|
[644] | 1403 | c |
---|
[766] | 1404 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) |
---|
| 1405 | CALL histwrite_phy(nid_mth,"cdrh",itau_w,cdragh) |
---|
[644] | 1406 | c |
---|
[766] | 1407 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) |
---|
| 1408 | CALL histwrite_phy(nid_mth,"cldl",itau_w,cldl) |
---|
[644] | 1409 | c |
---|
[766] | 1410 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) |
---|
| 1411 | CALL histwrite_phy(nid_mth,"cldm",itau_w,cldm) |
---|
[644] | 1412 | c |
---|
[766] | 1413 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) |
---|
| 1414 | CALL histwrite_phy(nid_mth,"cldh",itau_w,cldh) |
---|
[644] | 1415 | c |
---|
[766] | 1416 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt*100.,zx_tmp_2d) |
---|
| 1417 | CALL histwrite_phy(nid_mth,"cldt",itau_w,cldt*100.) |
---|
[644] | 1418 | c |
---|
[766] | 1419 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) |
---|
| 1420 | CALL histwrite_phy(nid_mth,"cldq",itau_w,cldq) |
---|
[644] | 1421 | c |
---|
| 1422 | zx_tmp_fi2d(1:klon) = flwp(1:klon) |
---|
[766] | 1423 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1424 | CALL histwrite_phy(nid_mth,"lwp",itau_w,zx_tmp_fi2d) |
---|
[644] | 1425 | c |
---|
| 1426 | zx_tmp_fi2d(1:klon) = fiwp(1:klon) |
---|
[766] | 1427 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1428 | CALL histwrite_phy(nid_mth,"iwp",itau_w,zx_tmp_fi2d) |
---|
[644] | 1429 | c |
---|
[766] | 1430 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) |
---|
| 1431 | CALL histwrite_phy(nid_mth,"ue",itau_w,ue) |
---|
[644] | 1432 | c |
---|
[766] | 1433 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) |
---|
| 1434 | CALL histwrite_phy(nid_mth,"ve",itau_w,ve) |
---|
[644] | 1435 | c |
---|
[766] | 1436 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, uq,zx_tmp_2d) |
---|
| 1437 | CALL histwrite_phy(nid_mth,"uq",itau_w,uq) |
---|
[644] | 1438 | c |
---|
[766] | 1439 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, vq,zx_tmp_2d) |
---|
| 1440 | CALL histwrite_phy(nid_mth,"vq",itau_w,vq) |
---|
[644] | 1441 | cKE43 |
---|
| 1442 | IF (iflag_con.GE.3) THEN ! sb |
---|
| 1443 | c |
---|
[766] | 1444 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, cape,zx_tmp_2d) |
---|
| 1445 | CALL histwrite_phy(nid_mth,"cape",itau_w,cape) |
---|
[644] | 1446 | c |
---|
[766] | 1447 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,pbase,zx_tmp_2d) |
---|
| 1448 | CALL histwrite_phy(nid_mth,"pbase",itau_w,pbase) |
---|
[644] | 1449 | c |
---|
[766] | 1450 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_pct,zx_tmp_2d) |
---|
| 1451 | CALL histwrite_phy(nid_mth,"ptop",itau_w,ema_pct) |
---|
[644] | 1452 | c |
---|
[766] | 1453 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,ema_cbmf,zx_tmp_2d) |
---|
| 1454 | CALL histwrite_phy(nid_mth,"fbase",itau_w,ema_cbmf) |
---|
[644] | 1455 | c |
---|
| 1456 | c34EK |
---|
| 1457 | c |
---|
[766] | 1458 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, prw,zx_tmp_2d) |
---|
| 1459 | CALL histwrite_phy(nid_mth,"prw",itau_w,prw) |
---|
[644] | 1460 | c |
---|
| 1461 | ENDIF !iflag_con.GE.3 |
---|
| 1462 | c |
---|
| 1463 | cIM cf. AM 081204 BEG |
---|
| 1464 | cHBTM2 |
---|
| 1465 | c |
---|
| 1466 | c |
---|
[766] | 1467 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblh,zx_tmp_2d) |
---|
| 1468 | CALL histwrite_phy(nid_mth,"s_pblh",itau_w,s_pblh) |
---|
[644] | 1469 | c |
---|
[766] | 1470 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_pblt,zx_tmp_2d) |
---|
| 1471 | CALL histwrite_phy(nid_mth,"s_pblt",itau_w,s_pblt) |
---|
[644] | 1472 | c |
---|
[766] | 1473 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_lcl,zx_tmp_2d) |
---|
| 1474 | CALL histwrite_phy(nid_mth,"s_lcl",itau_w,s_lcl) |
---|
[644] | 1475 | c |
---|
[766] | 1476 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_capCL,zx_tmp_2d) |
---|
| 1477 | CALL histwrite_phy(nid_mth,"s_capCL",itau_w,s_capCL) |
---|
[644] | 1478 | c |
---|
[766] | 1479 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_oliqCL,zx_tmp_2d) |
---|
| 1480 | CALL histwrite_phy(nid_mth,"s_oliqCL",itau_w,s_oliqCL) |
---|
[644] | 1481 | c |
---|
[766] | 1482 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_cteiCL,zx_tmp_2d) |
---|
| 1483 | CALL histwrite_phy(nid_mth,"s_cteiCL",itau_w,s_cteiCL) |
---|
[644] | 1484 | c |
---|
[766] | 1485 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_therm,zx_tmp_2d) |
---|
| 1486 | CALL histwrite_phy(nid_mth,"s_therm",itau_w,s_therm) |
---|
[644] | 1487 | c |
---|
[766] | 1488 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb1,zx_tmp_2d) |
---|
| 1489 | CALL histwrite_phy(nid_mth,"s_trmb1",itau_w,s_trmb1) |
---|
[644] | 1490 | c |
---|
[766] | 1491 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb2,zx_tmp_2d) |
---|
| 1492 | CALL histwrite_phy(nid_mth,"s_trmb2",itau_w,s_trmb2) |
---|
[644] | 1493 | c |
---|
[766] | 1494 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, s_trmb3,zx_tmp_2d) |
---|
| 1495 | CALL histwrite_phy(nid_mth,"s_trmb3",itau_w,s_trmb3) |
---|
[644] | 1496 | c |
---|
| 1497 | cIM cf. AM 081204 END |
---|
| 1498 | c |
---|
| 1499 | c Champs interpolles sur des niveaux de pression |
---|
| 1500 | c DO k=1, nlevENS |
---|
| 1501 | DO k=1, nlevSTD |
---|
| 1502 | c |
---|
| 1503 | c bb=clevSTD(k) |
---|
| 1504 | c |
---|
| 1505 | c IF(k.GE.2) THEN |
---|
| 1506 | IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) |
---|
| 1507 | IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) |
---|
| 1508 | c aa=clevSTD(k) |
---|
| 1509 | c bb=aa(1:lnblnk1(aa)) |
---|
| 1510 | c ENDIF |
---|
| 1511 | c |
---|
| 1512 | IF(bb2.EQ."850".OR.bb2.EQ."700" |
---|
| 1513 | $ .OR.bb2.EQ."500".OR.bb2.EQ."200") THEN |
---|
| 1514 | c |
---|
[766] | 1515 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,usumSTD(:,k,2),zx_tmp_2d) |
---|
| 1516 | CALL histwrite_phy(nid_mth,"u"//bb2,itau_w,usumSTD(:,k,2)) |
---|
[644] | 1517 | c |
---|
[766] | 1518 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,vsumSTD(:,k,2),zx_tmp_2d) |
---|
| 1519 | CALL histwrite_phy(nid_mth,"v"//bb2,itau_w,vsumSTD(:,k,2)) |
---|
[644] | 1520 | c |
---|
[766] | 1521 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,wsumSTD(:,k,2),zx_tmp_2d) |
---|
| 1522 | CALL histwrite_phy(nid_mth,"w"//bb2,itau_w,wsumSTD(:,k,2)) |
---|
[644] | 1523 | c |
---|
[766] | 1524 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,phisumSTD(:,k,2),zx_tmp_2d) |
---|
| 1525 | CALL histwrite_phy(nid_mth,"phi"//bb2,itau_w,phisumSTD(:,k,2)) |
---|
[644] | 1526 | c |
---|
| 1527 | ENDIF !(bb2.EQ."850".OR.bb2.EQ."700" |
---|
| 1528 | ENDDO !nlevSTD |
---|
| 1529 | c |
---|
| 1530 | DO i=1, klon |
---|
| 1531 | IF (pctsrf(i,is_oce).GT.epsfra) THEN |
---|
| 1532 | zx_tmp_fi2d(i) = fluxo(i) |
---|
| 1533 | ELSE |
---|
[684] | 1534 | zx_tmp_fi2d(i) = 0. |
---|
[644] | 1535 | ENDIF |
---|
| 1536 | ENDDO |
---|
| 1537 | c |
---|
[766] | 1538 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1539 | CALL histwrite_phy(nid_mth,"fluxo",itau_w,zx_tmp_fi2d) |
---|
[644] | 1540 | c |
---|
| 1541 | DO i=1, klon |
---|
| 1542 | IF (pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1543 | zx_tmp_fi2d(i) = fluxg(i) |
---|
| 1544 | ELSE |
---|
[684] | 1545 | zx_tmp_fi2d(i) = 0. |
---|
[644] | 1546 | ENDIF |
---|
| 1547 | ENDDO |
---|
| 1548 | c |
---|
[766] | 1549 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1550 | CALL histwrite_phy(nid_mth,"fluxg",itau_w,zx_tmp_fi2d) |
---|
[644] | 1551 | c |
---|
[684] | 1552 | DO i=1, klon |
---|
| 1553 | IF (pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 1554 | . pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1555 | zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ |
---|
| 1556 | . ftsol(i, is_sic) * pctsrf(i,is_sic))/ |
---|
| 1557 | . (pctsrf(i,is_oce)+pctsrf(i,is_sic)) |
---|
| 1558 | ELSE |
---|
| 1559 | zx_tmp_fi2d(i) = 273.15 |
---|
| 1560 | ENDIF |
---|
| 1561 | ENDDO |
---|
| 1562 | c |
---|
[766] | 1563 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1564 | CALL histwrite_phy(nid_mth,"t_oce_sic",itau_w,zx_tmp_fi2d) |
---|
[684] | 1565 | c |
---|
[644] | 1566 | IF (OCEAN.EQ.'force ') THEN |
---|
| 1567 | c |
---|
| 1568 | DO i=1, klon |
---|
| 1569 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 1570 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 1571 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 1572 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 1573 | ELSE |
---|
| 1574 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1575 | ENDIF |
---|
| 1576 | ENDDO |
---|
[766] | 1577 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1578 | CALL histwrite_phy(nid_mth,"lmt_bils",itau_w, |
---|
| 1579 | . zx_tmp_fi2d) |
---|
[644] | 1580 | c |
---|
| 1581 | ELSE IF (OCEAN.EQ.'slab ') THEN |
---|
| 1582 | c |
---|
| 1583 | c |
---|
| 1584 | DO i=1, klon |
---|
| 1585 | IF((pctsrf(i,is_oce).GT.epsfra).OR. |
---|
| 1586 | . (pctsrf(i,is_sic).GT.epsfra)) THEN |
---|
| 1587 | zx_tmp_fi2d(i) = (radsol(i) + fluxo(i))*pctsrf(i,is_oce)+ |
---|
| 1588 | . fluxg(i)*pctsrf(i,is_sic) |
---|
| 1589 | ELSE |
---|
| 1590 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1591 | ENDIF |
---|
| 1592 | ENDDO |
---|
[766] | 1593 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1594 | CALL histwrite_phy(nid_mth,"slab_bils",itau_w, |
---|
| 1595 | . zx_tmp_fi2d) |
---|
[644] | 1596 | c |
---|
| 1597 | DO i=1, klon |
---|
| 1598 | IF(pctsrf(i,is_oce).GT.epsfra.OR. |
---|
| 1599 | $ pctsrf(i,is_sic).GT.epsfra) THEN |
---|
| 1600 | zx_tmp_fi2d(i)=tslab(i) |
---|
| 1601 | ELSE |
---|
| 1602 | zx_tmp_fi2d(i) = 1.E+20 |
---|
| 1603 | ENDIF |
---|
| 1604 | ENDDO !i=1, klon |
---|
[766] | 1605 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) |
---|
| 1606 | CALL histwrite_phy(nid_mth,"tslab",itau_w,zx_tmp_fi2d) |
---|
[644] | 1607 | c |
---|
[766] | 1608 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice,zx_tmp_2d) |
---|
| 1609 | CALL histwrite_phy(nid_mth,"seaice",itau_w,seaice) |
---|
[644] | 1610 | c |
---|
[766] | 1611 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, seaice/1000.,zx_tmp_2d) |
---|
| 1612 | CALL histwrite_phy(nid_mth,"siceh",itau_w,seaice/1000.) |
---|
[644] | 1613 | c |
---|
| 1614 | ENDIF !OCEAN.EQ.force/slab |
---|
| 1615 | c |
---|
| 1616 | ENDIF !lev_histmth.GE.1 |
---|
| 1617 | c |
---|
| 1618 | IF(lev_histmth.GE.2) THEN |
---|
| 1619 | c |
---|
| 1620 | c Champs 3D: |
---|
| 1621 | C |
---|
| 1622 | c |
---|
[766] | 1623 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, flwc,zx_tmp_3d) |
---|
| 1624 | CALL histwrite_phy(nid_mth,"lwcon",itau_w,flwc) |
---|
[644] | 1625 | c |
---|
[766] | 1626 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d) |
---|
| 1627 | CALL histwrite_phy(nid_mth,"iwcon",itau_w,fiwc) |
---|
[644] | 1628 | c |
---|
[766] | 1629 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) |
---|
| 1630 | CALL histwrite_phy(nid_mth,"temp",itau_w,t_seri) |
---|
[524] | 1631 | c |
---|
[644] | 1632 | c#ifdef histmthNMC |
---|
| 1633 | c#include "write_histmthNMC.h" |
---|
| 1634 | c#endif |
---|
| 1635 | c |
---|
[766] | 1636 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) |
---|
| 1637 | CALL histwrite_phy(nid_mth,"ovap",itau_w, qx(:,:,ivap)) |
---|
[524] | 1638 | c |
---|
[644] | 1639 | c zx_tmp_fi3d(1:klon,1:klev)=qx(1:klon,1:klev,ivap)/ |
---|
| 1640 | c . (1-qx(1:klon,1:klev,ivap)) |
---|
| 1641 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
[766] | 1642 | c CALL histwrite_phy(nid_mth,"wvap",itau_w,zx_tmp_3d, |
---|
[644] | 1643 | c . iim*jjmp1*klev,ndex3d) |
---|
| 1644 | c |
---|
[766] | 1645 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) |
---|
| 1646 | CALL histwrite_phy(nid_mth,"geop",itau_w,zphi) |
---|
[524] | 1647 | c |
---|
[766] | 1648 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) |
---|
| 1649 | CALL histwrite_phy(nid_mth,"vitu",itau_w,u_seri) |
---|
[524] | 1650 | c |
---|
[766] | 1651 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) |
---|
| 1652 | CALL histwrite_phy(nid_mth,"vitv",itau_w,v_seri) |
---|
[644] | 1653 | c |
---|
[766] | 1654 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) |
---|
| 1655 | CALL histwrite_phy(nid_mth,"vitw",itau_w,omega) |
---|
[644] | 1656 | c |
---|
[766] | 1657 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) |
---|
| 1658 | CALL histwrite_phy(nid_mth,"pres",itau_w,pplay) |
---|
[644] | 1659 | c |
---|
[766] | 1660 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldfra*100., zx_tmp_3d) |
---|
| 1661 | CALL histwrite_phy(nid_mth,"rneb",itau_w,cldfra*100.) |
---|
[644] | 1662 | c |
---|
[766] | 1663 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, rnebcon, zx_tmp_3d) |
---|
| 1664 | CALL histwrite_phy(nid_mth,"rnebcon",itau_w,rnebcon) |
---|
[644] | 1665 | c |
---|
[766] | 1666 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_rh, zx_tmp_3d) |
---|
| 1667 | CALL histwrite_phy(nid_mth,"rhum",itau_w,zx_rh) |
---|
[644] | 1668 | c |
---|
| 1669 | DO k=1, klev |
---|
| 1670 | DO i=1, klon |
---|
| 1671 | cIM 151004 zx_tmp_fi3d(i,k)=MAX(wo(i,k),1.0e-12)*RG/46.6968 |
---|
| 1672 | c $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 1673 | c $ *(paprs(i,1)/101325.0) |
---|
| 1674 | c zx_tmp2_fi3d(i,k)=MAX(wo(i,k)*1.0e+6,1.0e-6)*RG/46.6968 |
---|
| 1675 | zx_tmp_fi3d(i,k)=wo(i,k)*RG/46.6968 |
---|
| 1676 | $ /(paprs(i,k)-paprs(i,k+1)) |
---|
| 1677 | $ *(paprs(i,1)/101325.0) |
---|
| 1678 | ENDDO !i |
---|
| 1679 | ENDDO !k |
---|
| 1680 | c |
---|
| 1681 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, wo, zx_tmp_3d) |
---|
| 1682 | c CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp2_fi3d,zx_tmp_3d) |
---|
[766] | 1683 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1684 | CALL histwrite_phy(nid_mth,"ozone",itau_w,zx_tmp_fi3d) |
---|
[644] | 1685 | c |
---|
[766] | 1686 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) |
---|
| 1687 | CALL histwrite_phy(nid_mth,"dtphy",itau_w,d_t) |
---|
[644] | 1688 | c |
---|
[766] | 1689 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_qx(:,:,ivap),zx_tmp_3d) |
---|
| 1690 | CALL histwrite_phy(nid_mth,"dqphy",itau_w,d_qx(:,:,ivap)) |
---|
[644] | 1691 | c |
---|
| 1692 | ENDIF !lev_histmth.GE.2 |
---|
| 1693 | c |
---|
| 1694 | IF(lev_histmth.GE.3) THEN |
---|
| 1695 | c |
---|
| 1696 | DO nsrf=1, nbsrf |
---|
| 1697 | c |
---|
| 1698 | zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) |
---|
[766] | 1699 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1700 | CALL histwrite_phy(nid_mth,"albe_"//clnsurf(nsrf),itau_w, |
---|
| 1701 | $ zx_tmp_fi2d) |
---|
[644] | 1702 | C |
---|
| 1703 | zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) |
---|
[766] | 1704 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d , zx_tmp_2d) |
---|
| 1705 | CALL histwrite_phy(nid_mth,"rugs_"//clnsurf(nsrf),itau_w, |
---|
| 1706 | $ zx_tmp_fi2d) |
---|
[644] | 1707 | c |
---|
| 1708 | zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf) |
---|
[766] | 1709 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, agesno,zx_tmp_2d) |
---|
| 1710 | CALL histwrite_phy(nid_mth,"ages_"//clnsurf(nsrf),itau_w |
---|
| 1711 | $ ,agesno) |
---|
[644] | 1712 | c |
---|
| 1713 | ENDDO !nsrf=1, nbsrf |
---|
| 1714 | c |
---|
[766] | 1715 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) |
---|
| 1716 | CALL histwrite_phy(nid_mth,"albs",itau_w,albsol) |
---|
[644] | 1717 | c |
---|
[766] | 1718 | cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) |
---|
| 1719 | CALL histwrite_phy(nid_mth,"albslw",itau_w,albsollw) |
---|
[644] | 1720 | c |
---|
| 1721 | ENDIF !lev_histmth.GE.3 |
---|
| 1722 | c |
---|
| 1723 | c |
---|
| 1724 | IF(lev_histmth.GE.4) THEN |
---|
| 1725 | c |
---|
[766] | 1726 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, clwcon0, zx_tmp_3d) |
---|
| 1727 | CALL histwrite_phy(nid_mth,"clwcon",itau_w,clwcon0) |
---|
[644] | 1728 | c |
---|
[653] | 1729 | IF (iflag_con.GE.3) THEN |
---|
| 1730 | c |
---|
[766] | 1731 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, Ma, zx_tmp_3d) |
---|
| 1732 | CALL histwrite_phy(nid_mth,"Ma",itau_w,Ma) |
---|
[644] | 1733 | c |
---|
[766] | 1734 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, upwd, zx_tmp_3d) |
---|
| 1735 | CALL histwrite_phy(nid_mth,"upwd",itau_w,upwd) |
---|
[644] | 1736 | c |
---|
[766] | 1737 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd, zx_tmp_3d) |
---|
| 1738 | CALL histwrite_phy(nid_mth,"dnwd",itau_w,dnwd) |
---|
[644] | 1739 | c |
---|
[766] | 1740 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, dnwd0, zx_tmp_3d) |
---|
| 1741 | CALL histwrite_phy(nid_mth,"dnwd0",itau_w,dnwd0) |
---|
[644] | 1742 | c |
---|
[653] | 1743 | ENDIF !iflag_con.GE.3 |
---|
[644] | 1744 | c |
---|
[766] | 1745 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) |
---|
| 1746 | CALL histwrite_phy(nid_mth,"dtdyn",itau_w,d_t_dyn) |
---|
[644] | 1747 | c |
---|
[766] | 1748 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_q_dyn, zx_tmp_3d) |
---|
| 1749 | CALL histwrite_phy(nid_mth,"dqdyn",itau_w,d_q_dyn) |
---|
[644] | 1750 | c |
---|
| 1751 | cIM: 101003 : K/30min ==> K/s |
---|
| 1752 | c |
---|
| 1753 | zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys |
---|
[766] | 1754 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1755 | CALL histwrite_phy(nid_mth,"dqcon",itau_w,zx_tmp_fi3d) |
---|
[644] | 1756 | c |
---|
| 1757 | cIM: 101003 : K/30min ==> K/s |
---|
| 1758 | zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ |
---|
| 1759 | . d_t_eva(1:klon,1:klev))/pdtphys |
---|
[766] | 1760 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zx_tmp_fi3d, zx_tmp_3d) |
---|
| 1761 | CALL histwrite_phy(nid_mth,"dtlschr",itau_w,zx_tmp_fi3d) |
---|
[644] | 1762 | c |
---|
| 1763 | zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys |
---|
[766] | 1764 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1765 | CALL histwrite_phy(nid_mth,"dqlsc",itau_w,zx_tmp_fi3d) |
---|
[644] | 1766 | c |
---|
| 1767 | cIM: 101003 : K/30min ==> K/s |
---|
| 1768 | zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 1769 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1770 | CALL histwrite_phy(nid_mth,"dtvdf",itau_w,zx_tmp_fi3d) |
---|
[644] | 1771 | c |
---|
| 1772 | zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 1773 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1774 | CALL histwrite_phy(nid_mth,"dqvdf",itau_w,zx_tmp_fi3d) |
---|
[644] | 1775 | c |
---|
| 1776 | cIM: 101003 : K/30min ==> K/s |
---|
| 1777 | zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys |
---|
[766] | 1778 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1779 | CALL histwrite_phy(nid_mth,"dteva",itau_w,zx_tmp_fi3d) |
---|
[644] | 1780 | c |
---|
| 1781 | zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys |
---|
[766] | 1782 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1783 | CALL histwrite_phy(nid_mth,"dqeva",itau_w,zx_tmp_fi3d) |
---|
[644] | 1784 | c |
---|
| 1785 | zpt_conv = 0. |
---|
| 1786 | where (ptconv) zpt_conv = 1. |
---|
[766] | 1787 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zpt_conv, zx_tmp_3d) |
---|
| 1788 | CALL histwrite_phy(nid_mth,"ptconv",itau_w,zpt_conv) |
---|
[644] | 1789 | c |
---|
[766] | 1790 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, ratqs, zx_tmp_3d) |
---|
| 1791 | CALL histwrite_phy(nid_mth,"ratqs",itau_w,ratqs) |
---|
[644] | 1792 | c |
---|
| 1793 | zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys |
---|
[766] | 1794 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1795 | CALL histwrite_phy(nid_mth,"dqajs",itau_w,zx_tmp_fi3d) |
---|
[644] | 1796 | c |
---|
| 1797 | cIM: 101003 : K/day ==> K/s |
---|
| 1798 | zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY |
---|
[766] | 1799 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1800 | CALL histwrite_phy(nid_mth,"dtswr",itau_w,zx_tmp_fi3d) |
---|
[644] | 1801 | c |
---|
| 1802 | cIM: 101003 : K/day ==> K/s |
---|
| 1803 | zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY |
---|
[766] | 1804 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1805 | CALL histwrite_phy(nid_mth,"dtsw0",itau_w,zx_tmp_fi3d) |
---|
[644] | 1806 | c |
---|
| 1807 | cIM: 101003 : K/day ==> K/s |
---|
| 1808 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY |
---|
[766] | 1809 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1810 | CALL histwrite_phy(nid_mth,"dtlwr",itau_w,zx_tmp_fi3d) |
---|
[644] | 1811 | c |
---|
| 1812 | cIM: 101003 : K/day ==> K/s |
---|
| 1813 | zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY |
---|
[766] | 1814 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1815 | CALL histwrite_phy(nid_mth,"dtlw0",itau_w,zx_tmp_fi3d) |
---|
[644] | 1816 | c |
---|
| 1817 | cIM: 101003 : deja en K/s |
---|
| 1818 | zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) |
---|
[766] | 1819 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1820 | CALL histwrite_phy(nid_mth,"dtec",itau_w,zx_tmp_fi3d) |
---|
[644] | 1821 | c |
---|
| 1822 | zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 1823 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1824 | CALL histwrite_phy(nid_mth,"duvdf",itau_w,zx_tmp_fi3d) |
---|
[644] | 1825 | c |
---|
| 1826 | zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys |
---|
[766] | 1827 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1828 | CALL histwrite_phy(nid_mth,"dvvdf",itau_w,zx_tmp_fi3d) |
---|
[644] | 1829 | c |
---|
| 1830 | IF (ok_orodr) THEN |
---|
| 1831 | zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys |
---|
[766] | 1832 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1833 | CALL histwrite_phy(nid_mth,"duoro",itau_w,zx_tmp_fi3d) |
---|
[644] | 1834 | c |
---|
| 1835 | zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys |
---|
[766] | 1836 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1837 | CALL histwrite_phy(nid_mth,"dvoro",itau_w,zx_tmp_fi3d) |
---|
[644] | 1838 | c |
---|
| 1839 | ENDIF |
---|
| 1840 | C |
---|
| 1841 | IF (ok_orolf) THEN |
---|
| 1842 | zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys |
---|
[766] | 1843 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1844 | CALL histwrite_phy(nid_mth,"dulif",itau_w,zx_tmp_fi3d) |
---|
[644] | 1845 | c |
---|
| 1846 | zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys |
---|
[766] | 1847 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) |
---|
| 1848 | CALL histwrite_phy(nid_mth,"dvlif",itau_w,zx_tmp_fi3d) |
---|
[644] | 1849 | ENDIF |
---|
| 1850 | C |
---|
[862] | 1851 | #if defined(INCA_AER) && defined(CPP_COUPLE) |
---|
| 1852 | CALL histwrite_phy(nid_mth,"tau",itau_w,cldtau) |
---|
| 1853 | CALL histwrite_phy(nid_mth,"CCM1",itau_w,ccm(:,:,1)) |
---|
| 1854 | CALL histwrite_phy(nid_mth,"CCM2",itau_w,ccm(:,:,2)) |
---|
| 1855 | CALL histwrite_phy(nid_mth,"taupi",itau_w,cldtaupi) |
---|
| 1856 | CALL histwrite_phy(nid_mth,"cldemi",itau_w,cldemi) |
---|
| 1857 | CALL histwrite_phy(nid_mth,"SWTOAAS_ZERO",itau_w,topsw_inca(:,1)) |
---|
| 1858 | CALL histwrite_phy(nid_mth,"SWTOAAS_AER",itau_w,topsw_inca(:,2)) |
---|
| 1859 | CALL histwrite_phy(nid_mth,"SWTOAAS_NAT",itau_w,topsw_inca(:,3)) |
---|
| 1860 | CALL histwrite_phy(nid_mth,"SWTOAAS_AD",itau_w,topswad_inca(:)) |
---|
| 1861 | CALL histwrite_phy(nid_mth,"SWTOAAS_AI",itau_w,topswai_inca(:)) |
---|
| 1862 | CALL histwrite_phy(nid_mth,"SWTOACS_ZERO",itau_w,topsw0_inca(:,1)) |
---|
| 1863 | CALL histwrite_phy(nid_mth,"SWTOACS_AER",itau_w,topsw0_inca(:,2)) |
---|
| 1864 | CALL histwrite_phy(nid_mth,"SWTOACS_AD",itau_w,topswad0_inca(:)) |
---|
| 1865 | CALL histwrite_phy(nid_mth,"SWTOACS_NAT",itau_w,topsw0_inca(:,3)) |
---|
| 1866 | #endif |
---|
| 1867 | #ifndef INCA |
---|
[524] | 1868 | IF (nqmax.GE.3) THEN |
---|
[653] | 1869 | DO iq=3,nqmax |
---|
[766] | 1870 | cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,iq), zx_tmp_3d) |
---|
| 1871 | CALL histwrite_phy(nid_mth,tnom(iq),itau_w,qx(:,:,iq)) |
---|
[524] | 1872 | ENDDO |
---|
| 1873 | ENDIF |
---|
| 1874 | c |
---|
[862] | 1875 | #endif |
---|
[524] | 1876 | ENDIF !lev_histmth.GE.4 |
---|
| 1877 | c |
---|
[644] | 1878 | ENDIF !(type_run.EQ."AMIP") |
---|
| 1879 | c |
---|
[524] | 1880 | if (ok_sync) then |
---|
[766] | 1881 | c$OMP MASTER |
---|
[524] | 1882 | call histsync(nid_mth) |
---|
[766] | 1883 | c$OMP END MASTER |
---|
[524] | 1884 | endif |
---|
[644] | 1885 | c |
---|
| 1886 | ENDIF !ok_mensuel |
---|