Changeset 218 for LMDZ.3.3/trunk/libf/dyn3d
- Timestamp:
- May 14, 2001, 5:01:36 PM (24 years ago)
- Location:
- LMDZ.3.3/trunk/libf/dyn3d
- Files:
-
- 5 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/trunk/libf/dyn3d/grilles_gcm_netcdf.F
r207 r218 168 168 c rajoute l'ecriture de la grille 169 169 170 status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg) 171 status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg) 172 status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg) 173 status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg) 170 #ifdef NC_DOUBLE 171 status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonuid,1,iim+1,rlonudeg) 172 status=NF_PUT_VARA_DOUBLE(ncid_out,out_lonvid,1,iim+1,rlonvdeg) 173 status=NF_PUT_VARA_DOUBLE(ncid_out,out_latuid,1,jjm+1,rlatudeg) 174 status=NF_PUT_VARA_DOUBLE(ncid_out,out_latvid,1,jjm,rlatvdeg) 175 #else 176 status=NF_PUT_VARA_REAL(ncid_out,out_lonuid,1,iim+1,rlonudeg) 177 status=NF_PUT_VARA_REAL(ncid_out,out_lonvid,1,iim+1,rlonvdeg) 178 status=NF_PUT_VARA_REAL(ncid_out,out_latuid,1,jjm+1,rlatudeg) 179 status=NF_PUT_VARA_REAL(ncid_out,out_latvid,1,jjm,rlatvdeg) 180 #endif 174 181 175 182 start(1)=1 … … 189 196 enddo 190 197 198 #ifdef NC_DOUBLE 199 status=NF_PUT_VARA_DOUBLE(ncid_out,out_varid,start, 200 s count,temp) 201 #else 191 202 status=NF_PUT_VARA_REAL(ncid_out,out_varid,start, 192 203 s count,temp) 204 #endif 205 193 206 194 207 c fermeture du fichier netcdf -
LMDZ.3.3/trunk/libf/dyn3d/read_dstoke.F
r207 r218 42 42 debut(2) = 1 43 43 44 #ifdef NC_DOUBLE 45 status=NF_GET_VARA_DOUBLE(ncidfd,variddt,debut,epais,dtv) 46 #else 44 47 status=NF_GET_VARA_REAL(ncidfd,variddt,debut,epais,dtv) 48 #endif 45 49 zdtvr=dtv(1,1,1) 46 50 51 #ifdef NC_DOUBLE 52 status=NF_GET_VARA_DOUBLE(ncidfd,varididvt,debut,epais,adv) 53 #else 47 54 status=NF_GET_VARA_REAL(ncidfd,varididvt,debut,epais,adv) 55 #endif 48 56 ziadvtrac= adv(1,1,1) 49 57 58 #ifdef NC_DOUBLE 59 status=NF_GET_VARA_DOUBLE(ncidfd,varididvp,debut,epais,adv2) 60 #else 50 61 status=NF_GET_VARA_REAL(ncidfd,varididvp,debut,epais,adv2) 62 #endif 51 63 ziadvtrac2= adv2(1,1,1) 52 64 -
LMDZ.3.3/trunk/libf/dyn3d/read_fstoke.F
r207 r218 28 28 . rlatu_dy(imo1,jmo1),rlatv_dy(imo1,jmo) 29 29 integer*4 ncrec,ncim,ncjm,nclm 30 realzrec,zim,zjm,zlm30 integer*4 zrec,zim,zjm,zlm 31 31 integer*4 xid,yid,zid,tid 32 32 real*4 zdtvr,ziadvtrac … … 126 126 c niveaux de pression 127 127 128 #ifdef NC_DOUBLE 129 status=NF_GET_VARA_DOUBLE(ncidf,varidpl,1,zlm,pl) 130 #else 128 131 status=NF_GET_VARA_REAL(ncidf,varidpl,1,zlm,pl) 132 #endif 129 133 130 134 c rlonu,rlonv,rlatu,rlatv … … 140 144 count(4)=0 141 145 142 status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy) 143 status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy) 146 #ifdef NC_DOUBLE 147 status=NF_GET_VARA_DOUBLE(ncidf,varidnlo,start,count,rlonu_dy) 148 status=NF_GET_VARA_DOUBLE(ncidf,varidnla,start,count,rlatu_dy) 149 #else 150 status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy) 151 status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy) 152 #endif 144 153 145 154 start(1)=1 … … 153 162 count(4)=0 154 163 155 status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy) 156 status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy) 164 #ifdef NC_DOUBLE 165 status=NF_GET_VARA_DOUBLE(ncidfv,varidnlov,start,count,rlonv_dy) 166 status=NF_GET_VARA_DOUBLE(ncidfv,varidnlav,start,count,rlatv_dy) 167 #else 168 status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy) 169 status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy) 170 #endif 157 171 158 172 c Lecture de phis et aire … … 169 183 170 184 c phis 185 #ifdef NC_DOUBLE 186 status=NF_GET_VARA_DOUBLE(ncidf,varidps,start,count,phis) 187 #else 171 188 status=NF_GET_VARA_REAL(ncidf,varidps,start,count,phis) 172 c print*,'WARNING!!! Correction bidon pour palier a un ' 173 c print*,'probleme dans la creation des fichiers nc' 174 c call correctbid(iim,jmo1*1,phis) 175 call dump2d(imo1,jmo1,phis,'PHISNC') 189 #endif 176 190 177 191 c aire 192 #ifdef NC_DOUBLE 193 status=NF_GET_VARA_DOUBLE(ncidf,varidai,start,count,airedy) 194 #else 178 195 status=NF_GET_VARA_REAL(ncidf,varidai,start,count,airedy) 179 c call correctbid(iim,jmo1*1,airedy) 180 c call dump2d(imo1,jmo1,airedy,'AIRENC') 181 182 c status= NF_CLOSE(ncidf) 196 #endif 183 197 184 198 else 185 199 186 200 print*,'ok1' 187 c stop'TEST lectfluxnc'188 201 189 202 c --------------------- … … 205 218 206 219 c masse 220 #ifdef NC_DOUBLE 221 status=NF_GET_VARA_DOUBLE(ncidf,varidm,start,count,masse) 222 #else 207 223 status=NF_GET_VARA_REAL(ncidf,varidm,start,count,masse) 224 #endif 208 225 c print*,'WARNING!!! Correction bidon pour palier a un ' 209 226 c print*,'probleme dans la creation des fichiers nc' … … 212 229 call dump2d(imo1,jmo1,masse,'MASSE COUCHE 1') 213 230 c pbaru 231 #ifdef NC_DOUBLE 232 status=NF_GET_VARA_DOUBLE(ncidf,varidpu,start,count,pbaru) 233 #else 214 234 status=NF_GET_VARA_REAL(ncidf,varidpu,start,count,pbaru) 235 #endif 215 236 c print*,'WARNING!!! Correction bidon pour palier a un ' 216 237 c print*,'probleme dans la creation des fichiers nc' … … 219 240 220 241 c w 242 #ifdef NC_DOUBLE 243 status=NF_GET_VARA_DOUBLE(ncidf,varidw,start,count,w) 244 #else 221 245 status=NF_GET_VARA_REAL(ncidf,varidw,start,count,w) 246 #endif 222 247 c print*,'WARNING!!! Correction bidon pour palier a un ' 223 248 c print*,'probleme dans la creation des fichiers nc' … … 226 251 227 252 c teta 253 #ifdef NC_DOUBLE 254 status=NF_GET_VARA_DOUBLE(ncidf,varidt,start,count,teta) 255 #else 228 256 status=NF_GET_VARA_REAL(ncidf,varidt,start,count,teta) 257 #endif 229 258 c print*,'WARNING!!! Correction bidon pour palier a un ' 230 259 c print*,'probleme dans la creation des fichiers nc' … … 233 262 234 263 c phi 264 #ifdef NC_DOUBLE 265 status=NF_GET_VARA_DOUBLE(ncidf,varidph,start,count,phi) 266 #else 235 267 status=NF_GET_VARA_REAL(ncidf,varidph,start,count,phi) 268 #endif 236 269 c print*,'WARNING!!! Correction bidon pour palier a un ' 237 270 c print*,'probleme dans la creation des fichiers nc' … … 243 276 c pbarv 244 277 print*,'Lecture de pbarv irec=',irec 278 #ifdef NC_DOUBLE 279 status=NF_GET_VARA_DOUBLE(ncidfv,varidpv,start,count,pbarv) 280 #else 245 281 status=NF_GET_VARA_REAL(ncidfv,varidpv,start,count,pbarv) 282 #endif 246 283 c call correctbid(iim,jjm*nlevnc,pbarv) 247 284 PRINT*,'*******LECTURE DE PBARV******' -
LMDZ.3.3/trunk/libf/dyn3d/read_fstoke0.F
r207 r218 24 24 . rlatu_dy(iip1,jjp1),rlatv_dy(iip1,jjm) 25 25 integer*4 ncrec,ncim,ncjm,nclm 26 realzrec,zim,zjm,zlm26 integer*4 zrec,zim,zjm,zlm 27 27 integer*4 xid,yid,zid,tid 28 28 real*4 zdtvr,ziadvtrac … … 122 122 c niveaux de pression 123 123 124 #ifdef NC_DOUBLE 125 status=NF_GET_VARA_DOUBLE(ncidf,varidpl,1,zlm,pl) 126 #else 124 127 status=NF_GET_VARA_REAL(ncidf,varidpl,1,zlm,pl) 128 #endif 125 129 126 130 c rlonu,rlonv,rlatu,rlatv … … 136 140 count(4)=0 137 141 138 status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy) 139 status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy) 140 142 #ifdef NC_DOUBLE 143 status=NF_GET_VARA_DOUBLE(ncidf,varidnlo,start,count,rlonu_dy) 144 status=NF_GET_VARA_DOUBLE(ncidf,varidnla,start,count,rlatu_dy) 145 #else 146 status=NF_GET_VARA_REAL(ncidf,varidnlo,start,count,rlonu_dy) 147 status=NF_GET_VARA_REAL(ncidf,varidnla,start,count,rlatu_dy) 148 #endif 149 141 150 start(1)=1 142 151 start(2)=1 … … 149 158 count(4)=0 150 159 151 status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy) 152 status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy) 160 #ifdef NC_DOUBLE 161 status=NF_GET_VARA_DOUBLE(ncidfv,varidnlov,start,count,rlonv_dy) 162 status=NF_GET_VARA_DOUBLE(ncidfv,varidnlav,start,count,rlatv_dy) 163 #else 164 status=NF_GET_VARA_REAL(ncidfv,varidnlov,start,count,rlonv_dy) 165 status=NF_GET_VARA_REAL(ncidfv,varidnlav,start,count,rlatv_dy) 166 #endif 153 167 154 168 c Lecture de phis et aire … … 165 179 166 180 c phis 181 #ifdef NC_DOUBLE 182 status=NF_GET_VARA_DOUBLE(ncidf,varidps,start,count,phis) 183 #else 167 184 status=NF_GET_VARA_REAL(ncidf,varidps,start,count,phis) 185 #endif 168 186 169 187 c aire 188 #ifdef NC_DOUBLE 189 status=NF_GET_VARA_DOUBLE(ncidf,varidai,start,count,airedy) 190 #else 170 191 status=NF_GET_VARA_REAL(ncidf,varidai,start,count,airedy) 192 #endif 171 193 172 194 else … … 193 215 c masse 194 216 PRINT*,'LECTURE DE masse' 217 #ifdef NC_DOUBLE 218 status=NF_GET_VARA_DOUBLE(ncidf,varidm,start,count,masse) 219 #else 195 220 status=NF_GET_VARA_REAL(ncidf,varidm,start,count,masse) 221 #endif 196 222 197 223 c pbaru 224 #ifdef NC_DOUBLE 225 status=NF_GET_VARA_DOUBLE(ncidf,varidpu,start,count,pbaru) 226 #else 198 227 status=NF_GET_VARA_REAL(ncidf,varidpu,start,count,pbaru) 228 #endif 199 229 200 230 c w 231 #ifdef NC_DOUBLE 232 status=NF_GET_VARA_DOUBLE(ncidf,varidw,start,count,w) 233 #else 201 234 status=NF_GET_VARA_REAL(ncidf,varidw,start,count,w) 235 #endif 202 236 203 237 c teta 238 #ifdef NC_DOUBLE 239 status=NF_GET_VARA_DOUBLE(ncidf,varidt,start,count,teta) 240 #else 204 241 status=NF_GET_VARA_REAL(ncidf,varidt,start,count,teta) 242 #endif 205 243 206 244 c phi 245 #ifdef NC_DOUBLE 246 status=NF_GET_VARA_DOUBLE(ncidf,varidph,start,count,phi) 247 #else 207 248 status=NF_GET_VARA_REAL(ncidf,varidph,start,count,phi) 249 #endif 208 250 209 251 count(2) = zjm-1 210 252 211 253 c pbarv 254 #ifdef NC_DOUBLE 255 status=NF_GET_VARA_DOUBLE(ncidfv,varidpv,start,count,pbarv) 256 #else 212 257 status=NF_GET_VARA_REAL(ncidfv,varidpv,start,count,pbarv) 258 #endif 213 259 214 260 start(3)=irec … … 222 268 return 223 269 end 224 -
LMDZ.3.3/trunk/libf/dyn3d/read_pstoke.F
r185 r218 196 196 c niveaux de pression 197 197 198 #ifdef NC_DOUBLE 199 status=NF_GET_VARA_DOUBLE(ncidp,varidpl,1,kev,pl) 200 #else 198 201 status=NF_GET_VARA_REAL(ncidp,varidpl,1,kev,pl) 202 #endif 199 203 200 204 c lecture de aire et phis … … 212 216 c 213 217 c phis 218 #ifdef NC_DOUBLE 219 status=NF_GET_VARA_DOUBLE(ncidp,varidps,start,count,phisfi2) 220 #else 214 221 status=NF_GET_VARA_REAL(ncidp,varidps,start,count,phisfi2) 222 #endif 223 224 215 225 call gr_ecrit_fi(1,kon,iim,jjm+1,phisfi2,phisfi) 216 226 217 227 c aire 228 #ifdef NC_DOUBLE 229 status=NF_GET_VARA_DOUBLE(ncidp,varidai,start,count,airefi2) 230 #else 218 231 status=NF_GET_VARA_REAL(ncidp,varidai,start,count,airefi2) 232 #endif 233 219 234 call gr_ecrit_fi(1,kon,iim,jjm+1,airefi2,airefi) 220 235 else … … 241 256 c frac_impa 242 257 258 #ifdef NC_DOUBLE 259 status=NF_GET_VARA_DOUBLE(ncidp,varidfi,start,count,frac_impa2) 260 #else 243 261 status=NF_GET_VARA_REAL(ncidp,varidfi,start,count,frac_impa2) 262 #endif 244 263 call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_impa2,frac_impa) 245 264 246 265 c frac_nucl 247 266 267 #ifdef NC_DOUBLE 268 status=NF_GET_VARA_DOUBLE(ncidp,varidfn,start,count,frac_nucl2) 269 #else 248 270 status=NF_GET_VARA_REAL(ncidp,varidfn,start,count,frac_nucl2) 271 #endif 249 272 call gr_ecrit_fi(kev,kon,iim,jjm+1,frac_nucl2,frac_nucl) 250 273 251 274 c abder t 275 276 #ifdef NC_DOUBLE 277 status=NF_GET_VARA_DOUBLE(ncidp,varidt,start,count,t2) 278 #else 252 279 status=NF_GET_VARA_REAL(ncidp,varidt,start,count,t2) 280 #endif 281 253 282 call gr_ecrit_fi(kev,kon,iim,jjm+1,t2,t) 254 283 255 284 c mfu 285 286 #ifdef NC_DOUBLE 287 status=NF_GET_VARA_DOUBLE(ncidp,varidmfu,start,count,mfu2) 288 #else 256 289 status=NF_GET_VARA_REAL(ncidp,varidmfu,start,count,mfu2) 290 #endif 291 257 292 call gr_ecrit_fi(kev,kon,iim,jjm+1,mfu2,mfu) 258 293 259 294 c mfd 295 #ifdef NC_DOUBLE 296 status=NF_GET_VARA_DOUBLE(ncidp,varidmfd,start,count,mfd2) 297 #else 260 298 status=NF_GET_VARA_REAL(ncidp,varidmfd,start,count,mfd2) 299 #endif 261 300 call gr_ecrit_fi(kev,kon,iim,jjm+1,mfd2,mfd) 262 301 263 302 c en_u 303 #ifdef NC_DOUBLE 304 status=NF_GET_VARA_DOUBLE(ncidp,varidenu,start,count,en_u2) 305 #else 264 306 status=NF_GET_VARA_REAL(ncidp,varidenu,start,count,en_u2) 307 #endif 265 308 call gr_ecrit_fi(kev,kon,iim,jjm+1,en_u2,en_u) 266 309 267 310 c de_u 311 #ifdef NC_DOUBLE 312 status=NF_GET_VARA_DOUBLE(ncidp,variddeu,start,count,de_u2) 313 #else 268 314 status=NF_GET_VARA_REAL(ncidp,variddeu,start,count,de_u2) 315 #endif 269 316 call gr_ecrit_fi(kev,kon,iim,jjm+1,de_u2,de_u) 270 317 271 318 c en_d 319 #ifdef NC_DOUBLE 320 status=NF_GET_VARA_DOUBLE(ncidp,varidend,start,count,en_d2) 321 #else 272 322 status=NF_GET_VARA_REAL(ncidp,varidend,start,count,en_d2) 323 #endif 273 324 call gr_ecrit_fi(kev,kon,iim,jjm+1,en_d2,en_d) 274 325 275 326 c de_d 327 #ifdef NC_DOUBLE 328 status=NF_GET_VARA_DOUBLE(ncidp,varidded,start,count,de_d2) 329 #else 276 330 status=NF_GET_VARA_REAL(ncidp,varidded,start,count,de_d2) 331 #endif 277 332 call gr_ecrit_fi(kev,kon,iim,jjm+1,de_d2,de_d) 278 333 279 334 c coefh 280 335 print*,'LECTURE de coefh a irec =',irec 336 #ifdef NC_DOUBLE 337 status=NF_GET_VARA_DOUBLE(ncidp,varidch,start,count,coefh2) 338 #else 281 339 status=NF_GET_VARA_REAL(ncidp,varidch,start,count,coefh2) 340 #endif 341 282 342 call gr_ecrit_fi(kev,kon,iim,jjm+1,coefh2,coefh) 283 343 … … 289 349 c pyu1 290 350 print*,'LECTURE de yu1 a irec =',irec 351 #ifdef NC_DOUBLE 352 status=NF_GET_VARA_DOUBLE(ncidp,varidyu1,start,count,pyu12) 353 #else 291 354 status=NF_GET_VARA_REAL(ncidp,varidyu1,start,count,pyu12) 355 #endif 356 292 357 call gr_ecrit_fi(1,kon,iim,jjm+1,pyu12,pyu1) 293 358 294 359 c pyv1 295 360 print*,'LECTURE de yv1 a irec =',irec 361 #ifdef NC_DOUBLE 362 status=NF_GET_VARA_DOUBLE(ncidp,varidyv1,start,count,pyv12) 363 #else 296 364 status=NF_GET_VARA_REAL(ncidp,varidyv1,start,count,pyv12) 365 #endif 297 366 call gr_ecrit_fi(1,kon,iim,jjm+1,pyv12,pyv1) 298 367 299 368 c ftsol1 300 369 print*,'LECTURE de ftsol1 a irec =',irec 370 #ifdef NC_DOUBLE 371 status=NF_GET_VARA_DOUBLE(ncidp,varidfts1,start,count,ftsol12) 372 #else 301 373 status=NF_GET_VARA_REAL(ncidp,varidfts1,start,count,ftsol12) 374 #endif 302 375 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol12,ftsol1) 303 376 304 377 c ftsol2 305 378 print*,'LECTURE de ftsol2 a irec =',irec 379 #ifdef NC_DOUBLE 380 status=NF_GET_VARA_DOUBLE(ncidp,varidfts2,start,count,ftsol22) 381 #else 306 382 status=NF_GET_VARA_REAL(ncidp,varidfts2,start,count,ftsol22) 383 #endif 307 384 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol22,ftsol2) 308 385 309 386 c ftsol3 310 387 print*,'LECTURE de ftsol3 a irec =',irec 388 #ifdef NC_DOUBLE 389 status=NF_GET_VARA_DOUBLE(ncidp,varidfts3,start,count,ftsol32) 390 #else 311 391 status=NF_GET_VARA_REAL(ncidp,varidfts3,start,count,ftsol32) 392 #endif 312 393 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol32,ftsol3) 313 394 314 395 c ftsol4 396 #ifdef NC_DOUBLE 397 status=NF_GET_VARA_DOUBLE(ncidp,varidfts4,start,count,ftsol42) 398 #else 315 399 status=NF_GET_VARA_REAL(ncidp,varidfts4,start,count,ftsol42) 400 #endif 316 401 call gr_ecrit_fi(1,kon,iim,jjm+1,ftsol42,ftsol4) 317 402 318 403 c psrf1 404 #ifdef NC_DOUBLE 405 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr1,start,count,psrf12) 406 #else 319 407 status=NF_GET_VARA_REAL(ncidp,varidpsr1,start,count,psrf12) 408 #endif 320 409 c call dump2d(iip1-1,jjm+1,psrf12,'PSRF1NC') 321 410 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf12,psrf1) 322 411 323 412 c psrf2 413 #ifdef NC_DOUBLE 414 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr2,start,count,psrf22) 415 #else 324 416 status=NF_GET_VARA_REAL(ncidp,varidpsr2,start,count,psrf22) 417 #endif 325 418 c call dump2d(iip1-1,jjm+1,psrf22,'PSRF2NC') 326 419 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf22,psrf2) 327 420 328 421 c psrf3 422 #ifdef NC_DOUBLE 423 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr3,start,count,psrf32) 424 #else 329 425 status=NF_GET_VARA_REAL(ncidp,varidpsr3,start,count,psrf32) 426 #endif 330 427 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf32,psrf3) 331 428 332 429 c psrf4 430 #ifdef NC_DOUBLE 431 status=NF_GET_VARA_DOUBLE(ncidp,varidpsr4,start,count,psrf42) 432 #else 333 433 status=NF_GET_VARA_REAL(ncidp,varidpsr4,start,count,psrf42) 434 #endif 334 435 call gr_ecrit_fi(1,kon,iim,jjm+1,psrf42,psrf4) 335 436
Note: See TracChangeset
for help on using the changeset viewer.