Changeset 97 for trunk/libf/phytitan/write_histday.h
- Timestamp:
- Mar 22, 2011, 5:25:44 PM (14 years ago)
- File:
-
- 1 edited
-
trunk/libf/phytitan/write_histday.h (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/libf/phytitan/write_histday.h
r3 r97 7 7 ndex3d = 0 8 8 zx_tmp_2d = 0. 9 z4_tmp_2d = 0.10 9 zx_tmp_3d = 0. 11 z4_tmp_3d = 0.12 10 zx_tmp_fi2d=0. 13 11 zx_tmp_fi3d=0. … … 24 22 c 25 23 CALL gr_fi_ecrit(1,klon,iim,jjmp1,pphis,zx_tmp_2d) 26 z4_tmp_2d=zx_tmp_2d 27 CALL histwrite(nid_day,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 24 CALL histwrite(nid_day,"phis",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 28 25 C 29 26 CALL gr_fi_ecrit(1,klon,iim,jjmp1,airephy,zx_tmp_2d) 30 z4_tmp_2d=zx_tmp_2d 31 CALL histwrite(nid_day,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 27 CALL histwrite(nid_day,"aire",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 32 28 c 33 29 ccccccc axe Ls 34 30 do j=1,jjmp1 35 31 do i=1,iim 36 z 4_tmp_2d(i,j)=zls*180./RPI ! zls est en radians !!32 zx_tmp_2d(i,j)=zls*180./RPI ! zls est en radians !! 37 33 enddo 38 34 enddo 39 CALL histwrite(nid_day,"ls",itau_w,z 4_tmp_2d,iim*jjmp1,ndex2d)35 CALL histwrite(nid_day,"ls",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 40 36 c 41 37 ccccccccccccc 2D fields, variables 42 38 c 43 39 CALL gr_fi_ecrit(1, klon,iim,jjmp1, ftsol,zx_tmp_2d) 44 z4_tmp_2d=zx_tmp_2d 45 CALL histwrite(nid_day,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 40 CALL histwrite(nid_day,"tsol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 46 41 c 47 42 DO i = 1, klon … … 49 44 ENDDO 50 45 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 51 z4_tmp_2d=zx_tmp_2d 52 CALL histwrite(nid_day,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 46 CALL histwrite(nid_day,"psol",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 53 47 c 54 48 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, ue,zx_tmp_2d) 55 c z4_tmp_2d=zx_tmp_2d 56 c CALL histwrite(nid_day,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 49 c CALL histwrite(nid_day,"ue",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 57 50 c 58 51 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, ve,zx_tmp_2d) 59 c z4_tmp_2d=zx_tmp_2d 60 c CALL histwrite(nid_day,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 52 c CALL histwrite(nid_day,"ve",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 61 53 c 62 54 ENDIF !lev_histday.GE.1 … … 68 60 c 69 61 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, t_seri, zx_tmp_3d) 70 z4_tmp_3d=zx_tmp_3d 71 CALL histwrite(nid_day,"temp",itau_w,z4_tmp_3d, 62 CALL histwrite(nid_day,"temp",itau_w,zx_tmp_3d, 72 63 . iim*jjmp1*klev,ndex3d) 73 64 c 74 65 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 75 z4_tmp_3d=zx_tmp_3d 76 CALL histwrite(nid_day,"pres",itau_w,z4_tmp_3d, 66 CALL histwrite(nid_day,"pres",itau_w,zx_tmp_3d, 77 67 . iim*jjmp1*klev,ndex3d) 78 68 c 79 69 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 80 z4_tmp_3d=zx_tmp_3d 81 CALL histwrite(nid_day,"geop",itau_w,z4_tmp_3d, 70 CALL histwrite(nid_day,"geop",itau_w,zx_tmp_3d, 82 71 . iim*jjmp1*klev,ndex3d) 83 72 c 84 73 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 85 z4_tmp_3d=zx_tmp_3d 86 CALL histwrite(nid_day,"vitu",itau_w,z4_tmp_3d, 74 CALL histwrite(nid_day,"vitu",itau_w,zx_tmp_3d, 87 75 . iim*jjmp1*klev,ndex3d) 88 76 c 89 77 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 90 z4_tmp_3d=zx_tmp_3d 91 CALL histwrite(nid_day,"vitv",itau_w,z4_tmp_3d, 78 CALL histwrite(nid_day,"vitv",itau_w,zx_tmp_3d, 92 79 . iim*jjmp1*klev,ndex3d) 93 80 c 94 81 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 95 z4_tmp_3d=zx_tmp_3d 96 CALL histwrite(nid_day,"vitw",itau_w,z4_tmp_3d, 82 CALL histwrite(nid_day,"vitw",itau_w,zx_tmp_3d, 97 83 . iim*jjmp1*klev,ndex3d) 98 84 c 99 85 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 100 z4_tmp_2d=zx_tmp_2d 101 CALL histwrite(nid_day,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 86 CALL histwrite(nid_day,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 102 87 c 103 88 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_dyn, zx_tmp_3d) 104 z4_tmp_3d=zx_tmp_3d 105 CALL histwrite(nid_day,"dudyn",itau_w,z4_tmp_3d, 89 CALL histwrite(nid_day,"dudyn",itau_w,zx_tmp_3d, 106 90 . iim*jjmp1*klev,ndex3d) 107 91 c 108 92 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) 109 z4_tmp_3d=zx_tmp_3d 110 CALL histwrite(nid_day,"duvdf",itau_w,z4_tmp_3d, 93 CALL histwrite(nid_day,"duvdf",itau_w,zx_tmp_3d, 111 94 . iim*jjmp1*klev,ndex3d) 112 95 c … … 115 98 DO iq=1,nmicro 116 99 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d) 117 z4_tmp_3d=zx_tmp_3d 118 CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d, 100 CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d, 119 101 . iim*jjmp1*klev,ndex3d) 120 102 ENDDO … … 123 105 DO iq=nmicro+1,nqmax 124 106 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d) 125 z4_tmp_3d=zx_tmp_3d 126 CALL histwrite(nid_day,tname(iq),itau_w,z4_tmp_3d, 107 CALL histwrite(nid_day,tname(iq),itau_w,zx_tmp_3d, 127 108 . iim*jjmp1*klev,ndex3d) 128 109 ENDDO … … 140 121 c 141 122 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 142 z4_tmp_2d=zx_tmp_2d 143 CALL histwrite(nid_day,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 123 CALL histwrite(nid_day,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 144 124 c 145 125 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 146 z4_tmp_2d=zx_tmp_2d 147 CALL histwrite(nid_day,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 126 CALL histwrite(nid_day,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 148 127 c 149 128 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 150 z4_tmp_2d=zx_tmp_2d 151 CALL histwrite(nid_day,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 129 CALL histwrite(nid_day,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 152 130 c 153 131 c 3D … … 155 133 zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev) 156 134 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 157 z4_tmp_3d=zx_tmp_3d 158 CALL histwrite(nid_day,"SWnet",itau_w,z4_tmp_3d, 135 CALL histwrite(nid_day,"SWnet",itau_w,zx_tmp_3d, 159 136 . iim*jjmp1*klev,ndex3d) 160 137 c 161 138 zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev) 162 139 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 163 z4_tmp_3d=zx_tmp_3d 164 CALL histwrite(nid_day,"LWnet",itau_w,z4_tmp_3d, 140 CALL histwrite(nid_day,"LWnet",itau_w,zx_tmp_3d, 165 141 . iim*jjmp1*klev,ndex3d) 166 142 c … … 171 147 c 172 148 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 173 z4_tmp_3d=zx_tmp_3d 174 CALL histwrite(nid_day,"dtdyn",itau_w,z4_tmp_3d, 149 CALL histwrite(nid_day,"dtdyn",itau_w,zx_tmp_3d, 175 150 . iim*jjmp1*klev,ndex3d) 176 151 c 177 152 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) 178 z4_tmp_3d=zx_tmp_3d 179 CALL histwrite(nid_day,"dtphy",itau_w,z4_tmp_3d, 153 CALL histwrite(nid_day,"dtphy",itau_w,zx_tmp_3d, 180 154 . iim*jjmp1*klev,ndex3d) 181 155 c K/s 182 156 zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev) 183 157 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 184 z4_tmp_3d=zx_tmp_3d 185 CALL histwrite(nid_day,"dtvdf",itau_w,z4_tmp_3d, 158 CALL histwrite(nid_day,"dtvdf",itau_w,zx_tmp_3d, 186 159 . iim*jjmp1*klev,ndex3d) 187 160 c … … 189 162 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev) 190 163 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 191 z4_tmp_3d=zx_tmp_3d 192 CALL histwrite(nid_day,"dtajs",itau_w,z4_tmp_3d, 164 CALL histwrite(nid_day,"dtajs",itau_w,zx_tmp_3d, 193 165 . iim*jjmp1*klev,ndex3d) 194 166 c … … 196 168 zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev) 197 169 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 198 z4_tmp_3d=zx_tmp_3d 199 CALL histwrite(nid_day,"dtswr",itau_w,z4_tmp_3d, 170 CALL histwrite(nid_day,"dtswr",itau_w,zx_tmp_3d, 200 171 . iim*jjmp1*klev,ndex3d) 201 172 c … … 203 174 zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev) 204 175 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 205 z4_tmp_3d=zx_tmp_3d 206 CALL histwrite(nid_day,"dtlwr",itau_w,z4_tmp_3d, 176 CALL histwrite(nid_day,"dtlwr",itau_w,zx_tmp_3d, 207 177 . iim*jjmp1*klev,ndex3d) 208 178 c K/s 209 179 zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) 210 180 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 211 z4_tmp_3d=zx_tmp_3d 212 CALL histwrite(nid_day,"dtec",itau_w,z4_tmp_3d, 181 CALL histwrite(nid_day,"dtec",itau_w,zx_tmp_3d, 213 182 . iim*jjmp1*klev,ndex3d) 214 183 c … … 220 189 c 221 190 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d) 222 c z4_tmp_2d=zx_tmp_2d223 191 c CALL histwrite(nid_day,"taux_",itau_w, 224 c $ z 4_tmp_2d,iim*jjmp1,ndex2d)192 c $ zx_tmp_2d,iim*jjmp1,ndex2d) 225 193 c 226 194 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d) 227 c z4_tmp_2d=zx_tmp_2d228 195 c CALL histwrite(nid_day,"tauy_",itau_w, 229 c $ z 4_tmp_2d,iim*jjmp1,ndex2d)196 c $ zx_tmp_2d,iim*jjmp1,ndex2d) 230 197 c 231 198 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 232 c z4_tmp_2d=zx_tmp_2d 233 c CALL histwrite(nid_day,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 199 c CALL histwrite(nid_day,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 234 200 c 235 201 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 236 c z4_tmp_2d=zx_tmp_2d 237 c CALL histwrite(nid_day,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 202 c CALL histwrite(nid_day,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 238 203 c 239 204 ENDIF !lev_histday.GE.5
Note: See TracChangeset
for help on using the changeset viewer.
