Changeset 97 for trunk/libf/phytitan/write_histins.h
- Timestamp:
- Mar 22, 2011, 5:25:44 PM (14 years ago)
- File:
-
- 1 edited
-
trunk/libf/phytitan/write_histins.h (modified) (13 diffs)
Legend:
- Unmodified
- Added
- Removed
-
trunk/libf/phytitan/write_histins.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_ins,"phis",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 24 CALL histwrite(nid_ins,"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_ins,"aire",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 27 CALL histwrite(nid_ins,"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_ins,"ls",itau_w,z 4_tmp_2d,iim*jjmp1,ndex2d)35 CALL histwrite(nid_ins,"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_ins,"tsol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 40 CALL histwrite(nid_ins,"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_ins,"psol",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 46 CALL histwrite(nid_ins,"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_ins,"ue",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 49 c CALL histwrite(nid_ins,"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_ins,"ve",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 52 c CALL histwrite(nid_ins,"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_ins,"temp",itau_w,z4_tmp_3d, 62 CALL histwrite(nid_ins,"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_ins,"pres",itau_w,z4_tmp_3d, 66 CALL histwrite(nid_ins,"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_ins,"geop",itau_w,z4_tmp_3d, 70 CALL histwrite(nid_ins,"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_ins,"vitu",itau_w,z4_tmp_3d, 74 CALL histwrite(nid_ins,"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_ins,"vitv",itau_w,z4_tmp_3d, 78 CALL histwrite(nid_ins,"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_ins,"vitw",itau_w,z4_tmp_3d, 82 CALL histwrite(nid_ins,"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_ins,"tops",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 86 CALL histwrite(nid_ins,"tops",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 102 87 c 103 88 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_u_vdf, zx_tmp_3d) 104 z4_tmp_3d=zx_tmp_3d 105 CALL histwrite(nid_ins,"duvdf",itau_w,z4_tmp_3d, 89 CALL histwrite(nid_ins,"duvdf",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_dyn, zx_tmp_3d) 109 z4_tmp_3d=zx_tmp_3d 110 CALL histwrite(nid_ins,"dudyn",itau_w,z4_tmp_3d, 93 CALL histwrite(nid_ins,"dudyn",itau_w,zx_tmp_3d, 111 94 . iim*jjmp1*klev,ndex3d) 112 95 c … … 122 105 DO iq=1,nmicro 123 106 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qaer(1,1,iq), zx_tmp_3d) 124 z4_tmp_3d=zx_tmp_3d 125 CALL histwrite(nid_ins,tname(iq),itau_w,z4_tmp_3d, 107 CALL histwrite(nid_ins,tname(iq),itau_w,zx_tmp_3d, 126 108 . iim*jjmp1*klev,ndex3d) 127 109 ENDDO … … 130 112 DO iq=nmicro+1,nqmax 131 113 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,tr_seri(1,1,iq),zx_tmp_3d) 132 z4_tmp_3d=zx_tmp_3d 133 CALL histwrite(nid_ins,tname(iq),itau_w,z4_tmp_3d, 114 CALL histwrite(nid_ins,tname(iq),itau_w,zx_tmp_3d, 134 115 . iim*jjmp1*klev,ndex3d) 135 116 ENDDO … … 142 123 c 143 124 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 144 z4_tmp_2d=zx_tmp_2d 145 CALL histwrite(nid_ins,"topl",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 125 CALL histwrite(nid_ins,"topl",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 146 126 c 147 127 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 148 z4_tmp_2d=zx_tmp_2d 149 CALL histwrite(nid_ins,"sols",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 128 CALL histwrite(nid_ins,"sols",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 150 129 c 151 130 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 152 z4_tmp_2d=zx_tmp_2d 153 CALL histwrite(nid_ins,"soll",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 131 CALL histwrite(nid_ins,"soll",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 154 132 c 155 133 c 3D … … 157 135 zx_tmp_fi3d(1:klon,1:klev)=swnet(1:klon,1:klev) 158 136 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 159 z4_tmp_3d=zx_tmp_3d 160 CALL histwrite(nid_ins,"SWnet",itau_w,z4_tmp_3d, 137 CALL histwrite(nid_ins,"SWnet",itau_w,zx_tmp_3d, 161 138 . iim*jjmp1*klev,ndex3d) 162 139 c 163 140 zx_tmp_fi3d(1:klon,1:klev)=lwnet(1:klon,1:klev) 164 141 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 165 z4_tmp_3d=zx_tmp_3d 166 CALL histwrite(nid_ins,"LWnet",itau_w,z4_tmp_3d, 142 CALL histwrite(nid_ins,"LWnet",itau_w,zx_tmp_3d, 167 143 . iim*jjmp1*klev,ndex3d) 168 144 c … … 173 149 c 174 150 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t_dyn, zx_tmp_3d) 175 z4_tmp_3d=zx_tmp_3d 176 CALL histwrite(nid_ins,"dtdyn",itau_w,z4_tmp_3d, 151 CALL histwrite(nid_ins,"dtdyn",itau_w,zx_tmp_3d, 177 152 . iim*jjmp1*klev,ndex3d) 178 153 c 179 154 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_t, zx_tmp_3d) 180 z4_tmp_3d=zx_tmp_3d 181 CALL histwrite(nid_ins,"dtphy",itau_w,z4_tmp_3d, 155 CALL histwrite(nid_ins,"dtphy",itau_w,zx_tmp_3d, 182 156 . iim*jjmp1*klev,ndex3d) 183 157 c K/s 184 158 zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev) 185 159 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 186 z4_tmp_3d=zx_tmp_3d 187 CALL histwrite(nid_ins,"dtvdf",itau_w,z4_tmp_3d, 160 CALL histwrite(nid_ins,"dtvdf",itau_w,zx_tmp_3d, 188 161 . iim*jjmp1*klev,ndex3d) 189 162 c … … 191 164 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev) 192 165 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 193 z4_tmp_3d=zx_tmp_3d 194 CALL histwrite(nid_ins,"dtajs",itau_w,z4_tmp_3d, 166 CALL histwrite(nid_ins,"dtajs",itau_w,zx_tmp_3d, 195 167 . iim*jjmp1*klev,ndex3d) 196 168 c … … 198 170 zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev) 199 171 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 200 z4_tmp_3d=zx_tmp_3d 201 CALL histwrite(nid_ins,"dtswr",itau_w,z4_tmp_3d, 172 CALL histwrite(nid_ins,"dtswr",itau_w,zx_tmp_3d, 202 173 . iim*jjmp1*klev,ndex3d) 203 174 c … … 205 176 zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev) 206 177 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 207 z4_tmp_3d=zx_tmp_3d 208 CALL histwrite(nid_ins,"dtlwr",itau_w,z4_tmp_3d, 178 CALL histwrite(nid_ins,"dtlwr",itau_w,zx_tmp_3d, 209 179 . iim*jjmp1*klev,ndex3d) 210 180 c K/s 211 181 zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev) 212 182 CALL gr_fi_ecrit(klev,klon,iim,jjmp1,zx_tmp_fi3d,zx_tmp_3d) 213 z4_tmp_3d=zx_tmp_3d 214 CALL histwrite(nid_ins,"dtec",itau_w,z4_tmp_3d, 183 CALL histwrite(nid_ins,"dtec",itau_w,zx_tmp_3d, 215 184 . iim*jjmp1*klev,ndex3d) 216 185 c 217 186 c CALL gr_fi_ecrit(klev,klon,iim,jjmp1, d_v_vdf, zx_tmp_3d) 218 c z4_tmp_3d=zx_tmp_3d 219 c CALL histwrite(nid_ins,"dvvdf",itau_w,z4_tmp_3d, 187 c CALL histwrite(nid_ins,"dvvdf",itau_w,zx_tmp_3d, 220 188 c . iim*jjmp1*klev,ndex3d) 221 189 c … … 227 195 c 228 196 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxu , zx_tmp_2d) 229 c z4_tmp_2d=zx_tmp_2d230 197 c CALL histwrite(nid_ins,"taux_",itau_w, 231 c $ z 4_tmp_2d,iim*jjmp1,ndex2d)198 c $ zx_tmp_2d,iim*jjmp1,ndex2d) 232 199 c 233 200 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, fluxv , zx_tmp_2d) 234 c z4_tmp_2d=zx_tmp_2d235 201 c CALL histwrite(nid_ins,"tauy_",itau_w, 236 c $ z 4_tmp_2d,iim*jjmp1,ndex2d)202 c $ zx_tmp_2d,iim*jjmp1,ndex2d) 237 203 c 238 204 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) 239 c z4_tmp_2d=zx_tmp_2d 240 c CALL histwrite(nid_ins,"cdrm",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 205 c CALL histwrite(nid_ins,"cdrm",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 241 206 c 242 207 c CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragh,zx_tmp_2d) 243 c z4_tmp_2d=zx_tmp_2d 244 c CALL histwrite(nid_ins,"cdrh",itau_w,z4_tmp_2d,iim*jjmp1,ndex2d) 208 c CALL histwrite(nid_ins,"cdrh",itau_w,zx_tmp_2d,iim*jjmp1,ndex2d) 245 209 c 246 210 ENDIF !lev_histday.GE.5
Note: See TracChangeset
for help on using the changeset viewer.
