- Timestamp:
- Jul 19, 2024, 5:41:58 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/cosp/cosp_output_write_mod.F90
r4619 r5082 185 185 do k=1,PARASOL_NREFL 186 186 do ip=1, Npoints 187 if (stlidar%cldlayer(ip,4) .gt.0.01.and.stlidar%parasolrefl(ip,k).ne.missing_val) then187 if (stlidar%cldlayer(ip,4)>0.01.and.stlidar%parasolrefl(ip,k)/=missing_val) then 188 188 parasolcrefl(ip,k)=(stlidar%parasolrefl(ip,k)-0.03*(1.-stlidar%cldlayer(ip,4)))/ & 189 189 stlidar%cldlayer(ip,4) … … 456 456 CHARACTER(LEN=20) :: typeecrit 457 457 458 ! ug On r écupère le type écrit de la structure:459 ! Assez moche, Ã| refaire si meilleure méthode...458 ! ug On récupère le type écrit de la structure: 459 ! Assez moche, �| refaire si meilleure méthode... 460 460 IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN 461 461 typeecrit = 'once' … … 523 523 524 524 ! Axe vertical 525 IF (nvertsave .eq.nvertp(iff)) THEN525 IF (nvertsave==nvertp(iff)) THEN 526 526 klevs=PARASOL_NREFL 527 527 nam_axvert="sza" 528 ELSE IF (nvertsave .eq.nvertisccp(iff)) THEN528 ELSE IF (nvertsave==nvertisccp(iff)) THEN 529 529 klevs=7 530 530 nam_axvert="pressure2" 531 ELSE IF (nvertsave .eq.nvertcol(iff)) THEN531 ELSE IF (nvertsave==nvertcol(iff)) THEN 532 532 klevs=Ncolout 533 533 nam_axvert="column" 534 ELSE IF (nvertsave .eq.nverttemp(iff)) THEN534 ELSE IF (nvertsave==nverttemp(iff)) THEN 535 535 klevs=LIDAR_NTEMP 536 536 nam_axvert="temp" 537 ELSE IF (nvertsave .eq.nvertmisr(iff)) THEN537 ELSE IF (nvertsave==nvertmisr(iff)) THEN 538 538 klevs=MISR_N_CTH 539 539 nam_axvert="cth16" 540 ELSE IF (nvertsave .eq.nvertReffIce(iff)) THEN540 ELSE IF (nvertsave==nvertReffIce(iff)) THEN 541 541 klevs= numMODISReffIceBins 542 542 nam_axvert="ReffIce" 543 ELSE IF (nvertsave .eq.nvertReffLiq(iff)) THEN543 ELSE IF (nvertsave==nvertReffLiq(iff)) THEN 544 544 klevs= numMODISReffLiqBins 545 545 nam_axvert="ReffLiq" … … 558 558 END IF 559 559 560 ! ug On r écupère le type écrit de la structure:561 ! Assez moche, Ã| refaire si meilleure méthode...560 ! ug On récupère le type écrit de la structure: 561 ! Assez moche, �| refaire si meilleure méthode... 562 562 IF (INDEX(var%cosp_typeecrit(iff), "once") > 0) THEN 563 563 typeecrit = 'once' … … 628 628 IF (prt_level >= 9) WRITE(lunout,*)'Begin histrwrite2d ',var%name 629 629 630 ! On regarde si on est dans la phase de d éfinition ou d'écriture:630 ! On regarde si on est dans la phase de définition ou d'écriture: 631 631 IF(.NOT.cosp_varsdefined) THEN 632 632 !$OMP MASTER 633 !Si phase de d éfinition.... on définit633 !Si phase de définition.... on définit 634 634 CALL conf_cospoutputs(var%name,var%cles) 635 635 DO iff=1, 3 … … 640 640 !$OMP END MASTER 641 641 ELSE 642 !Et sinon on.... écrit642 !Et sinon on.... écrit 643 643 IF (SIZE(field)/=klon) & 644 644 CALL abort_physic('iophy::histwrite2d_cosp','Field first DIMENSION not equal to klon',1) … … 725 725 nom=var%name 726 726 END IF 727 ! On regarde si on est dans la phase de d éfinition ou d'écriture:727 ! On regarde si on est dans la phase de définition ou d'écriture: 728 728 IF(.NOT.cosp_varsdefined) THEN 729 !Si phase de d éfinition.... on définit729 !Si phase de définition.... on définit 730 730 !$OMP MASTER 731 731 CALL conf_cospoutputs(var%name,var%cles) … … 737 737 !$OMP END MASTER 738 738 ELSE 739 !Et sinon on.... écrit739 !Et sinon on.... écrit 740 740 IF (SIZE(field,1)/=klon) & 741 741 CALL abort_physic('iophy::histwrite3d','Field first DIMENSION not equal to klon',1) … … 809 809 810 810 IF(cosp_varsdefined) THEN 811 !Et sinon on.... écrit811 !Et sinon on.... écrit 812 812 IF (SIZE(field,1)/=klon) & 813 813 CALL abort_physic('iophy::histwrite3d','Field first DIMENSION not equal to klon',1)
Note: See TracChangeset
for help on using the changeset viewer.