- Timestamp:
- Feb 22, 2021, 12:44:07 PM (4 years ago)
- Location:
- LMDZ6/branches/LMDZ-tracers
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/LMDZ-tracers
- Property svn:mergeinfo changed
-
LMDZ6/branches/LMDZ-tracers/libf/phylmd/cospv2/lmdz_cosp_output_write_mod.F90
r3723 r3851 273 273 endif 274 274 if (cfg%LlidarBetaMol532) then 275 where(cospOUT%calipso_ tau_tot == R_UNDEF) cospOUT%calipso_tau_tot= missing_val275 where(cospOUT%calipso_beta_mol == R_UNDEF) cospOUT%calipso_beta_mol = missing_val 276 276 CALL histwrite3d_cosp(o_lidarBetaMol532,cospOUT%calipso_beta_mol,nvertmcosp) 277 277 endif … … 282 282 !!!! Sorties Ground Lidar 283 283 if (cfg%LgrLidar532) then 284 285 where(cospOUT%grLidar532_beta_tot == R_UNDEF) cospOUT%grLidar532_beta_tot = missing_val 286 where(cospOUT%grLidar532_cfad_sr == R_UNDEF) cospOUT%grLidar532_cfad_sr = missing_val 287 where(cospOUT%grLidar532_lidarcld == R_UNDEF) cospOUT%grLidar532_lidarcld = missing_val 288 where(cospOUT%grLidar532_cldlayer == R_UNDEF) cospOUT%grLidar532_cldlayer = missing_val 289 where(cospOUT%grLidar532_beta_mol == R_UNDEF) cospOUT%grLidar532_beta_mol = missing_val 290 where(cospOUT%grLidar532_srbval == R_UNDEF) cospOUT%grLidar532_srbval = missing_val 291 292 if (cfg%LcllgrLidar532) CALL histwrite2d_cosp(o_cllgrLidar532,cospOUT%grLidar532_cldlayer(:,1)) 293 if (cfg%LclmgrLidar532) CALL histwrite2d_cosp(o_clmgrLidar532,cospOUT%grLidar532_cldlayer(:,2)) 294 if (cfg%LclhgrLidar532) CALL histwrite2d_cosp(o_clhgrLidar532,cospOUT%grLidar532_cldlayer(:,3)) 295 if (cfg%LcltgrLidar532) CALL histwrite2d_cosp(o_cltgrLidar532,cospOUT%grLidar532_cldlayer(:,4)) 296 297 if (cfg%LclgrLidar532) CALL histwrite3d_cosp(o_clgrLidar532,cospOUT%grLidar532_lidarcld,nvert) 298 if (cfg%LlidarBetaMol532gr) CALL histwrite3d_cosp(o_lidarBetaMol532gr,cospOUT%grLidar532_beta_mol,nvertmcosp) 299 300 do icl=1,SR_BINS 284 ! AI juin 2020 Voir a quoi correspond ce champs 285 ! where(cospOUT%grLidar532_srbval == R_UNDEF) cospOUT%grLidar532_srbval = missing_val 286 287 if (cfg%LcllgrLidar532) then 288 where(cospOUT%grLidar532_cldlayer(:,1) == R_UNDEF) cospOUT%grLidar532_cldlayer(:,1) = missing_val 289 CALL histwrite2d_cosp(o_cllgrLidar532,cospOUT%grLidar532_cldlayer(:,1)) 290 endif 291 if (cfg%LclmgrLidar532) then 292 where(cospOUT%grLidar532_cldlayer(:,2) == R_UNDEF) cospOUT%grLidar532_cldlayer(:,2) = missing_val 293 CALL histwrite2d_cosp(o_clmgrLidar532,cospOUT%grLidar532_cldlayer(:,2)) 294 endif 295 if (cfg%LclhgrLidar532) then 296 where(cospOUT%grLidar532_cldlayer(:,3) == R_UNDEF) cospOUT%grLidar532_cldlayer(:,3) = missing_val 297 CALL histwrite2d_cosp(o_clhgrLidar532,cospOUT%grLidar532_cldlayer(:,3)) 298 endif 299 if (cfg%LcltgrLidar532) then 300 where(cospOUT%grLidar532_cldlayer(:,4) == R_UNDEF) cospOUT%grLidar532_cldlayer(:,4) = missing_val 301 CALL histwrite2d_cosp(o_cltgrLidar532,cospOUT%grLidar532_cldlayer(:,4)) 302 endif 303 304 if (cfg%LclgrLidar532) then 305 where(cospOUT%grLidar532_lidarcld == R_UNDEF) cospOUT%grLidar532_lidarcld = missing_val 306 CALL histwrite3d_cosp(o_clgrLidar532,cospOUT%grLidar532_lidarcld,nvert) 307 endif 308 if (cfg%LlidarBetaMol532gr) then 309 where(cospOUT%grLidar532_beta_mol == R_UNDEF) cospOUT%grLidar532_beta_mol = missing_val 310 CALL histwrite3d_cosp(o_lidarBetaMol532gr,cospOUT%grLidar532_beta_mol,nvertmcosp) 311 endif 312 if (cfg%LcfadLidarsr532gr) then 313 where(cospOUT%grLidar532_cfad_sr == R_UNDEF) cospOUT%grLidar532_cfad_sr = missing_val 314 do icl=1,SR_BINS 301 315 do k=1,Nlvgrid 302 316 do ip=1,Npoints … … 304 318 enddo 305 319 enddo 306 enddo 307 if (cfg%LcfadLidarsr532gr) CALL histwrite4d_cosp(o_cfadLidarsr532gr,tmp_fi4da_cfadLgr) 308 309 if (cfg%Latb532gr) CALL histwrite4d_cosp(o_atb532gr,cospOUT%grLidar532_beta_tot) 310 320 enddo 321 CALL histwrite4d_cosp(o_cfadLidarsr532gr,tmp_fi4da_cfadLgr) 322 endif 323 324 if (cfg%Latb532gr) then 325 where(cospOUT%grLidar532_beta_tot == R_UNDEF) cospOUT%grLidar532_beta_tot = missing_val 326 CALL histwrite4d_cosp(o_atb532gr,cospOUT%grLidar532_beta_tot) 327 endif 311 328 endif ! Ground Lidar 532 nm 312 329 … … 314 331 !!!! Sorties Atlid 315 332 if (cfg%Latlid) then 316 317 where(cospOUT%atlid_beta_tot == R_UNDEF) cospOUT%atlid_beta_tot = missing_val 318 where(cospOUT%atlid_cfad_sr == R_UNDEF) cospOUT%atlid_cfad_sr = missing_val 319 where(cospOUT%atlid_lidarcld == R_UNDEF) cospOUT%atlid_lidarcld = missing_val 320 where(cospOUT%atlid_cldlayer == R_UNDEF) cospOUT%atlid_cldlayer = missing_val 321 where(cospOUT%atlid_beta_mol == R_UNDEF) cospOUT%atlid_beta_mol = missing_val 322 where(cospOUT%atlid_srbval == R_UNDEF) cospOUT%atlid_srbval = missing_val 323 324 if (cfg%Lcllatlid) CALL histwrite2d_cosp(o_cllatlid,cospOUT%atlid_cldlayer(:,1)) 325 if (cfg%Lclmatlid) CALL histwrite2d_cosp(o_clmatlid,cospOUT%atlid_cldlayer(:,2)) 326 if (cfg%Lclhatlid) CALL histwrite2d_cosp(o_clhatlid,cospOUT%atlid_cldlayer(:,3)) 327 if (cfg%Lcltatlid) CALL histwrite2d_cosp(o_cltatlid,cospOUT%atlid_cldlayer(:,4)) 328 329 if (cfg%Lclatlid) CALL histwrite3d_cosp(o_clatlid,cospOUT%atlid_lidarcld,nvert) 330 if (cfg%LlidarBetaMol355) CALL histwrite3d_cosp(o_lidarBetaMol355,cospOUT%atlid_beta_mol,nvertmcosp) 331 332 do icl=1,SR_BINS 333 ! AI juin 2020 Voir a quoi correspond ce champs 334 ! where(cospOUT%atlid_srbval == R_UNDEF) cospOUT%atlid_srbval = missing_val 335 336 if (cfg%Lcllatlid) then 337 where(cospOUT%atlid_cldlayer(:,1) == R_UNDEF) cospOUT%atlid_cldlayer(:,1) = missing_val 338 CALL histwrite2d_cosp(o_cllatlid,cospOUT%atlid_cldlayer(:,1)) 339 endif 340 if (cfg%Lclmatlid) then 341 where(cospOUT%atlid_cldlayer(:,2) == R_UNDEF) cospOUT%atlid_cldlayer(:,2) = missing_val 342 CALL histwrite2d_cosp(o_clmatlid,cospOUT%atlid_cldlayer(:,2)) 343 endif 344 if (cfg%Lclhatlid) then 345 where(cospOUT%atlid_cldlayer(:,3) == R_UNDEF) cospOUT%atlid_cldlayer(:,3) = missing_val 346 CALL histwrite2d_cosp(o_clhatlid,cospOUT%atlid_cldlayer(:,3)) 347 endif 348 if (cfg%Lcltatlid) then 349 where(cospOUT%atlid_cldlayer(:,4) == R_UNDEF) cospOUT%atlid_cldlayer(:,4) = missing_val 350 CALL histwrite2d_cosp(o_cltatlid,cospOUT%atlid_cldlayer(:,4)) 351 endif 352 if (cfg%Lclatlid) then 353 where(cospOUT%atlid_lidarcld == R_UNDEF) cospOUT%atlid_lidarcld = missing_val 354 CALL histwrite3d_cosp(o_clatlid,cospOUT%atlid_lidarcld,nvert) 355 endif 356 if (cfg%LlidarBetaMol355) then 357 where(cospOUT%atlid_beta_mol == R_UNDEF) cospOUT%atlid_beta_mol = missing_val 358 CALL histwrite3d_cosp(o_lidarBetaMol355,cospOUT%atlid_beta_mol,nvertmcosp) 359 endif 360 if (cfg%LcfadLidarsr355) then 361 where(cospOUT%atlid_cfad_sr == R_UNDEF) cospOUT%atlid_cfad_sr = missing_val 362 do icl=1,SR_BINS 333 363 do k=1,Nlvgrid 334 364 do ip=1,Npoints … … 336 366 enddo 337 367 enddo 338 enddo 339 if (cfg%LcfadLidarsr355) CALL histwrite4d_cosp(o_cfadlidarsr355,tmp_fi4da_cfadLatlid) 340 341 if (cfg%Latb355) CALL histwrite4d_cosp(o_atb355,cospOUT%atlid_beta_tot) 342 368 enddo 369 CALL histwrite4d_cosp(o_cfadlidarsr355,tmp_fi4da_cfadLatlid) 370 endif 371 372 if (cfg%Latb355) then 373 where(cospOUT%atlid_beta_tot == R_UNDEF) cospOUT%atlid_beta_tot = missing_val 374 CALL histwrite4d_cosp(o_atb355,cospOUT%atlid_beta_tot) 375 endif 343 376 endif ! Atlid 344 377 … … 346 379 if (cfg%Lparasol) then 347 380 if (cfg%LparasolRefl) then 348 ! Ces 2 diagnostics sont controles par la clef logique "LparasolRefl" 349 350 !!! if (cfg%LparasolRefl) CALL histwrite3d_cosp(o_parasolrefl,cospOUT%parasolrefl,nvertp) 381 where(cospOUT%parasolGrid_refl == R_UNDEF) cospOUT%parasolGrid_refl = missing_val 382 where(cospOUT%parasolPix_refl == R_UNDEF) cospOUT%parasolPix_refl = missing_val 351 383 CALL histwrite3d_cosp(o_parasolGrid_refl,cospOUT%parasolGrid_refl,nvertp) 352 353 384 CALL histwrite4d_cosp(o_parasolPix_refl,cospOUT%parasolPix_refl) 354 355 endif ! LparasolRefl 356 endif ! Parasol 385 endif ! LparasolRefl 386 endif ! Parasol 357 387 358 388 ! if (cfg%LparasolRefl) then … … 376 406 !!! Sorties CloudSat 377 407 if (cfg%Lcloudsat) then 378 379 where(cospOUT%cloudsat_Ze_tot == R_UNDEF) cospOUT%cloudsat_Ze_tot = missing_val 380 where(cospOUT%cloudsat_cfad_ze == R_UNDEF) cospOUT%cloudsat_cfad_ze = missing_val 381 where(cospOUT%cloudsat_precip_cover == R_UNDEF) cospOUT%cloudsat_precip_cover = missing_val 382 where(cospOUT%cloudsat_pia == R_UNDEF) cospOUT%cloudsat_pia = missing_val 383 384 if (cfg%Lptradarflag0) CALL histwrite2d_cosp(o_ptradarflag0,cospOUT%cloudsat_precip_cover(:,1)) 385 if (cfg%Lptradarflag1) CALL histwrite2d_cosp(o_ptradarflag1,cospOUT%cloudsat_precip_cover(:,2)) 386 if (cfg%Lptradarflag2) CALL histwrite2d_cosp(o_ptradarflag2,cospOUT%cloudsat_precip_cover(:,3)) 387 if (cfg%Lptradarflag3) CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,4)) 388 if (cfg%Lptradarflag4) CALL histwrite2d_cosp(o_ptradarflag4,cospOUT%cloudsat_precip_cover(:,5)) 389 if (cfg%Lptradarflag5) CALL histwrite2d_cosp(o_ptradarflag5,cospOUT%cloudsat_precip_cover(:,6)) 390 if (cfg%Lptradarflag6) CALL histwrite2d_cosp(o_ptradarflag6,cospOUT%cloudsat_precip_cover(:,7)) 391 if (cfg%Lptradarflag7) CALL histwrite2d_cosp(o_ptradarflag7,cospOUT%cloudsat_precip_cover(:,8)) 392 if (cfg%Lptradarflag8) CALL histwrite2d_cosp(o_ptradarflag8,cospOUT%cloudsat_precip_cover(:,9)) 393 if (cfg%Lptradarflag9) CALL histwrite2d_cosp(o_ptradarflag9,cospOUT%cloudsat_precip_cover(:,10)) 394 if (cfg%Lradarpia) CALL histwrite2d_cosp(o_radarpia,cospOUT%cloudsat_pia) 395 396 do icl=1,CLOUDSAT_DBZE_BINS 408 ! AI juin 2020 voir a quoi correspond ce champs 409 ! where(cospOUT%cloudsat_pia == R_UNDEF) cospOUT%cloudsat_pia = missing_val 410 411 if (cfg%Lptradarflag0) then 412 where(cospOUT%cloudsat_precip_cover(:,1) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,1) = missing_val 413 CALL histwrite2d_cosp(o_ptradarflag0,cospOUT%cloudsat_precip_cover(:,1)) 414 endif 415 if (cfg%Lptradarflag1) then 416 where(cospOUT%cloudsat_precip_cover(:,2) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,2) = missing_val 417 CALL histwrite2d_cosp(o_ptradarflag1,cospOUT%cloudsat_precip_cover(:,2)) 418 endif 419 if (cfg%Lptradarflag2) then 420 where(cospOUT%cloudsat_precip_cover(:,3) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,3) = missing_val 421 CALL histwrite2d_cosp(o_ptradarflag2,cospOUT%cloudsat_precip_cover(:,3)) 422 endif 423 if (cfg%Lptradarflag3) then 424 where(cospOUT%cloudsat_precip_cover(:,4) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,4) = missing_val 425 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,4)) 426 endif 427 if (cfg%Lptradarflag4) then 428 where(cospOUT%cloudsat_precip_cover(:,5) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,5) = missing_val 429 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,5)) 430 endif 431 if (cfg%Lptradarflag5) then 432 where(cospOUT%cloudsat_precip_cover(:,6) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,6) = missing_val 433 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,6)) 434 endif 435 if (cfg%Lptradarflag6) then 436 where(cospOUT%cloudsat_precip_cover(:,7) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,7) = missing_val 437 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,7)) 438 endif 439 if (cfg%Lptradarflag7) then 440 where(cospOUT%cloudsat_precip_cover(:,8) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,8) = missing_val 441 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,8)) 442 endif 443 if (cfg%Lptradarflag8) then 444 where(cospOUT%cloudsat_precip_cover(:,9) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,9) = missing_val 445 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,9)) 446 endif 447 if (cfg%Lptradarflag9) then 448 where(cospOUT%cloudsat_precip_cover(:,10) == R_UNDEF) cospOUT%cloudsat_precip_cover(:,10) = missing_val 449 CALL histwrite2d_cosp(o_ptradarflag3,cospOUT%cloudsat_precip_cover(:,10)) 450 endif 451 452 if (cfg%Ldbze94) then 453 where(cospOUT%cloudsat_Ze_tot == R_UNDEF) cospOUT%cloudsat_Ze_tot = missing_val 454 CALL histwrite4d_cosp(o_dbze94,cospOUT%cloudsat_Ze_tot) 455 endif 456 if (cfg%LcfadDbze94) then 457 where(cospOUT%cloudsat_cfad_ze == R_UNDEF) cospOUT%cloudsat_cfad_ze = missing_val 458 do icl=1,CLOUDSAT_DBZE_BINS 397 459 do k=1,Nlvgrid 398 460 do ip=1,Npoints … … 400 462 enddo 401 463 enddo 402 enddo 403 if (cfg%Ldbze94) CALL histwrite4d_cosp(o_dbze94,cospOUT%cloudsat_Ze_tot) 404 ! if (cfg%LcfadDbze94) CALL histwrite4d_cosp(o_cfadDbze94,stradar%cfad_ze) 405 if (cfg%LcfadDbze94) CALL histwrite4d_cosp(o_cfadDbze94,tmp_fi4da_cfadR) 464 enddo 465 CALL histwrite4d_cosp(o_cfadDbze94,tmp_fi4da_cfadR) 466 endif 406 467 endif 407 468 ! endif pour CloudSat … … 410 471 !!! Sorties combinees Cloudsat et Calipso 411 472 if (cfg%Lcalipso .and. cfg%Lcloudsat) then 412 where(cospOUT%lidar_only_freq_cloud == R_UNDEF) & 473 474 if (cfg%Lclcalipso2) then 475 where(cospOUT%lidar_only_freq_cloud == R_UNDEF) & 413 476 cospOUT%lidar_only_freq_cloud = missing_val 414 where(cospOUT%cloudsat_tcc == R_UNDEF) & 477 CALL histwrite3d_cosp(o_clcalipso2,cospOUT%lidar_only_freq_cloud,nvert) 478 endif 479 if (cfg%Lcloudsat_tcc) then 480 where(cospOUT%cloudsat_tcc == R_UNDEF) & 415 481 cospOUT%cloudsat_tcc = missing_val 416 where(cospOUT%cloudsat_tcc2 == R_UNDEF) & 482 CALL histwrite2d_cosp(o_cloudsat_tcc,cospOUT%cloudsat_tcc) 483 endif 484 if (cfg%Lcloudsat_tcc2) then 485 where(cospOUT%cloudsat_tcc2 == R_UNDEF) & 417 486 cospOUT%cloudsat_tcc2 = missing_val 418 where(cospOUT%radar_lidar_tcc == R_UNDEF) &419 cospOUT%radar_lidar_tcc = missing_val420 421 if (cfg%Lclcalipso2) CALL histwrite3d_cosp(o_clcalipso2,cospOUT%lidar_only_freq_cloud,nvert)422 if (cfg%Lcloudsat_tcc) CALL histwrite2d_cosp(o_cloudsat_tcc,cospOUT%cloudsat_tcc)423 if (cfg%Lcloudsat_tcc2) CALL histwrite2d_cosp(o_cloudsat_tcc2,cospOUT%cloudsat_tcc2)424 if (cfg%Lcltlidarradar) CALL histwrite2d_cosp(o_cltlidarradar,cospOUT%radar_lidar_tcc)487 CALL histwrite2d_cosp(o_cloudsat_tcc2,cospOUT%cloudsat_tcc2) 488 endif 489 if (cfg%Lcltlidarradar) then 490 where(cospOUT%radar_lidar_tcc == R_UNDEF) & 491 cospOUT%radar_lidar_tcc = missing_val 492 CALL histwrite2d_cosp(o_cltlidarradar,cospOUT%radar_lidar_tcc) 493 endif 425 494 endif 426 495
Note: See TracChangeset
for help on using the changeset viewer.