Changeset 918 for LMDZ4/branches/LMDZ4_V3_patches
- Timestamp:
- Feb 27, 2008, 4:24:26 PM (16 years ago)
- Location:
- LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd
- Files:
-
- 12 edited
-
ini_histISCCP.h (modified) (2 diffs)
-
ini_histday.h (modified) (7 diffs)
-
ini_histhf.h (modified) (5 diffs)
-
ini_histmth.h (modified) (11 diffs)
-
isccp_cloud_types.F (modified) (2 diffs)
-
pbl_surface_mod.F90 (modified) (1 diff)
-
physiq.F (modified) (4 diffs)
-
phytrac.F (modified) (2 diffs)
-
surf_landice_mod.F90 (modified) (1 diff)
-
write_histISCCP.h (modified) (3 diffs)
-
write_histday.h (modified) (2 diffs)
-
write_histhf.h (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histISCCP.h
r845 r918 146 146 . "ave(X)", zcals(n),zout_isccp(n)) 147 147 c 148 CALL histdef(nid_isccp, "ptop", 149 . "Cld top pressure (cf ISCCP simulator)","Pa", 150 . iim,jj_nb,nhori, 1,1,1,-99,32, 151 . "ave(X)", zcals(n),zout_isccp(n)) 152 c 153 CALL histdef(nid_isccp,"meantaucld", 154 . "ISCCP mean cloud optical thickness","1", 155 . iim,jj_nb,nhori, 1,1,1, -99, 32, 156 . "ave(X)", zcals(n),zout_isccp(n)) 148 157 ENDDO 149 158 c … … 169 178 . "Nb of calls with sunlit ", "%", 170 179 . iim, jj_nb,nhori,1,1,1,-99,32, 180 . "ave(X)", zcals(n),zout_isccp(n)) 181 c 182 CALL histdef(nid_isccp, "ptop", 183 . "Cld top pressure (cf ISCCP simulator)","Pa", 184 . iim,jj_nb,nhori, 1,1,1,-99,32, 185 . "ave(X)", zcals(n),zout_isccp(n)) 186 c 187 CALL histdef(nid_isccp,"meantaucld", 188 . "ISCCP mean cloud optical thickness","1", 189 . iim,jj_nb,nhori, 1,1,1, -99, 32, 171 190 . "ave(X)", zcals(n),zout_isccp(n)) 172 191 c -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histday.h
r897 r918 39 39 . "once", zstophy,zout) 40 40 c 41 CALL histdef(nid_day, "aire", "Grid area", " -",41 CALL histdef(nid_day, "aire", "Grid area", "m2", 42 42 . iim,jj_nb,nhori, 1,1,1, -99, 32, 43 43 . "once", zstophy,zout) 44 44 c 45 CALL histdef(nid_day, "contfracATM","% sfce ter+lic "," -",45 CALL histdef(nid_day, "contfracATM","% sfce ter+lic ","1", 46 46 . iim,jj_nb,nhori, 1,1,1, -99, 32, 47 47 . "once", zstophy,zout) 48 48 c 49 CALL histdef(nid_day, "contfracOR","% sfce terre OR", " -",49 CALL histdef(nid_day, "contfracOR","% sfce terre OR", "1", 50 50 . iim,jj_nb,nhori, 1,1,1, -99, 32, 51 51 . "once", zstophy,zout) … … 100 100 . "ave(X)", zstophy,zout) 101 101 c 102 CALL histdef(nid_day, "sicf", "Sea-ice fraction", " -",102 CALL histdef(nid_day, "sicf", "Sea-ice fraction", "1", 103 103 . iim,jj_nb,nhori, 1,1,1, -99, 32, 104 104 . "inst(X)", zstoday,zout) … … 167 167 . "ave(X)", zstophy,zout) 168 168 c 169 CALL histdef(nid_day, "cldl", "Low-level cloudiness", " -",170 . iim,jj_nb,nhori, 1,1,1, -99, 32, 171 . "ave(X)", zstophy,zout) 172 c 173 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", " -",174 . iim,jj_nb,nhori, 1,1,1, -99, 32, 175 . "ave(X)", zstophy,zout) 176 c 177 CALL histdef(nid_day, "cldh", "High-level cloudiness", " -",178 . iim,jj_nb,nhori, 1,1,1, -99, 32, 179 . "ave(X)", zstophy,zout) 180 c 181 CALL histdef(nid_day, "cldt", "Total cloudiness", " -",169 CALL histdef(nid_day, "cldl", "Low-level cloudiness", "1", 170 . iim,jj_nb,nhori, 1,1,1, -99, 32, 171 . "ave(X)", zstophy,zout) 172 c 173 CALL histdef(nid_day, "cldm", "Mid-level cloudiness", "1", 174 . iim,jj_nb,nhori, 1,1,1, -99, 32, 175 . "ave(X)", zstophy,zout) 176 c 177 CALL histdef(nid_day, "cldh", "High-level cloudiness", "1", 178 . iim,jj_nb,nhori, 1,1,1, -99, 32, 179 . "ave(X)", zstophy,zout) 180 c 181 CALL histdef(nid_day, "cldt", "Total cloudiness", "1", 182 182 . iim,jj_nb,nhori, 1,1,1, -99, 32, 183 183 . "ave(X)", zstophy,zout) … … 375 375 . "ave(X)", zstophy,zout) 376 376 c 377 cIM IF (iflag_con.GE.3) THEN378 IF (ok_isccp) THEN379 CALL histdef(nid_day, "ptop",380 c . "Cld top pressure","mb",381 . "Cld top pressure (cf ISCCP simulator)","Pa",382 . iim,jj_nb,nhori, 1,1,1,-99,383 . 32, "ave(X)", zstophy,zout)384 ENDIF385 cIM ENDIF !(iflag_con.GE.3) THEN386 c387 377 ENDIF !lev_histday.GE.2 388 378 c================================================================= … … 390 380 c================================================================= 391 381 c INITIALISATION DES CHAMPS 3D 392 c=================================================================393 c=================================================================394 382 c================================================================= 395 383 c Champs 3D: … … 586 574 . "ave(X)", zstophy,zout) 587 575 c 588 CALL histdef(nid_day, "rneb", "Cloud fraction", " -",576 CALL histdef(nid_day, "rneb", "Cloud fraction", "1", 589 577 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 590 578 . "ave(X)", zstophy,zout) 591 579 c 592 580 CALL histdef(nid_day, "rnebcon", "Convective Cloud Fraction" 593 . , " -",581 . , "1", 594 582 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 595 583 . "ave(X)", zstophy,zout) … … 604 592 . "ave(X)", zstophy,zout) 605 593 c 606 IF (ok_isccp) THEN607 CALL histdef(nid_day,"meantaucld",608 . "ISCCP mean cloud optical thickness","1",609 . iim,jj_nb,nhori, 1,1,1, -99, 32,610 . "ave(X)", zstophy,zout)611 ENDIF612 c613 594 CALL histdef(nid_day,"cldtau","Cloud optical thickness","1", 614 595 . iim,jj_nb,nhori, klev,1,klev, nvert, 32, -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histhf.h
r897 r918 39 39 IF(lev_histhf.GE.1) THEN 40 40 c 41 CALL histdef(nid_hf, "aireTER","Grid area CONT"," -",41 CALL histdef(nid_hf, "aireTER","Grid area CONT","m2", 42 42 . iim,jj_nb,nhori, 1,1,1, -99, 32, 43 43 . "once", zstohf,zout) 44 44 c 45 CALL histdef(nid_hf, "contfracATM","% sfce ter+lic "," -",45 CALL histdef(nid_hf, "contfracATM","% sfce ter+lic ","1", 46 46 . iim,jj_nb,nhori, 1,1,1, -99, 32, 47 47 . "once", zstohf,zout) 48 48 c 49 CALL histdef(nid_hf, "contfracOR","% sfce terre OR", " -",49 CALL histdef(nid_hf, "contfracOR","% sfce terre OR", "1", 50 50 . iim,jj_nb,nhori, 1,1,1, -99, 32, 51 51 . "once", zstohf,zout) … … 93 93 c 94 94 bb2=clevSTD(k) 95 c 96 c IF(k.GE.2) THEN 97 c aa=clevSTD(k) 98 c bb=aa(1:lnblnk1(aa)) 99 c ENDIF 100 c 101 bb2=clevSTD(k) 95 102 96 IF(bb2.EQ."500") THEN 103 97 CALL histdef(nid_hf, "phi"//bb2, … … 112 106 IF(lev_histhf.GE.2) THEN 113 107 c 114 CALL histdef(nid_hf, "cldt", "Total cloudiness", " %",108 CALL histdef(nid_hf, "cldt", "Total cloudiness", "1", 115 109 . iim,jj_nb,nhori, 1,1,1, -99, 32, 116 110 . "ave(X)", zstophy,zout) … … 137 131 c 138 132 DO k=1, nlevSTD 139 c140 c bb=clevSTD(k)141 c142 c IF(k.GE.2) THEN143 c IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k)144 c IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k)145 c aa=clevSTD(k)146 c bb=aa(1:lnblnk1(aa))147 c ENDIF148 c149 133 bb2=clevSTD(k) 150 134 IF(bb2.EQ."1000".OR.bb2.EQ."925".OR. … … 335 319 . "inst(X)", zstohf,zout) 336 320 c 337 CALL histdef(nid_hf, "cldt", "Total cloudiness", " %",321 CALL histdef(nid_hf, "cldt", "Total cloudiness", "1", 338 322 . iim,jj_nb,nhori, 1,1,1, -99, 32, 339 323 . "ave(X)", zstohf,zout) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/ini_histmth.h
r906 r918 43 43 . "once", zstophy,zout) 44 44 c 45 CALL histdef(nid_mth, "aire", "Grid area", " -",45 CALL histdef(nid_mth, "aire", "Grid area", "m2", 46 46 . iim,jj_nb,nhori, 1,1,1, -99, 32, 47 47 . "once", zstophy,zout) … … 91 91 . "ave(X)", zstophy,zout) 92 92 c 93 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", " -",93 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "1", 94 94 . iim,jj_nb,nhori, 1,1,1, -99, 32, 95 95 . "ave(X)", zstophy,zout) … … 123 123 c 124 124 CALL histdef(nid_mth, "ndayrain", 125 . "Number of day with rain (liq+sol)", " -",125 . "Number of day with rain (liq+sol)", "days", 126 126 . iim,jj_nb,nhori, 1,1,1, -99, 32, 127 127 . "inst(X)", zstomth,zout) … … 294 294 c Effets des aerosols 295 295 c 296 cIF (ok_ade.OR.ok_aie) THEN297 CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2",298 . iim,jj_nb,nhori, 1,1,1, -99, 32,299 . "ave(X)", zstophy,zout)300 c 301 CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2",302 . iim,jj_nb,nhori, 1,1,1, -99, 32,303 . "ave(X)", zstophy,zout)304 c 305 CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2",306 . iim,jj_nb,nhori, 1,1,1, -99, 32,307 . "ave(X)", zstophy,zout)308 c 309 CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2",310 . iim,jj_nb,nhori, 1,1,1, -99, 32,311 . "ave(X)", zstophy,zout)312 c endif 296 IF (ok_ade.OR.ok_aie) THEN 297 CALL histdef(nid_mth, "topsad", "ADE at TOA", "W/m2", 298 . iim,jj_nb,nhori, 1,1,1, -99, 32, 299 . "ave(X)", zstophy,zout) 300 c 301 CALL histdef(nid_mth, "solsad", "ADE at sfc", "W/m2", 302 . iim,jj_nb,nhori, 1,1,1, -99, 32, 303 . "ave(X)", zstophy,zout) 304 c 305 CALL histdef(nid_mth, "topsai", "AIE at TOA", "W/m2", 306 . iim,jj_nb,nhori, 1,1,1, -99, 32, 307 . "ave(X)", zstophy,zout) 308 c 309 CALL histdef(nid_mth, "solsai", "AIE at sfc", "W/m2", 310 . iim,jj_nb,nhori, 1,1,1, -99, 32, 311 . "ave(X)", zstophy,zout) 312 ENDIF 313 313 c 314 314 c CALL histdef(nid_mth, "frtu", "Zonal wind stress", "Pa", … … 419 419 . "ave(X)", zstophy,zout) 420 420 c 421 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", " -",422 . iim,jj_nb,nhori, 1,1,1, -99, 32, 423 . "ave(X)", zstophy,zout) 424 c 425 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", " -",426 . iim,jj_nb,nhori, 1,1,1, -99, 32, 427 . "ave(X)", zstophy,zout) 428 c 429 CALL histdef(nid_mth, "cldh", "High-level cloudiness", " -",430 . iim,jj_nb,nhori, 1,1,1, -99, 32, 431 . "ave(X)", zstophy,zout) 432 c 433 CALL histdef(nid_mth, "cldt", "Total cloudiness", " -",421 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "1", 422 . iim,jj_nb,nhori, 1,1,1, -99, 32, 423 . "ave(X)", zstophy,zout) 424 c 425 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "1", 426 . iim,jj_nb,nhori, 1,1,1, -99, 32, 427 . "ave(X)", zstophy,zout) 428 c 429 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "1", 430 . iim,jj_nb,nhori, 1,1,1, -99, 32, 431 . "ave(X)", zstophy,zout) 432 c 433 CALL histdef(nid_mth, "cldt", "Total cloudiness", "1", 434 434 . iim,jj_nb,nhori, 1,1,1, -99, 32, 435 435 . "ave(X)", zstophy,zout) … … 658 658 . "ave(X)", zstophy,zout) 659 659 c 660 CALL histdef(nid_mth, "rneb", "Cloud fraction", " -",660 CALL histdef(nid_mth, "rneb", "Cloud fraction", "1", 661 661 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 662 662 . "ave(X)", zstophy,zout) 663 663 c 664 664 CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction" 665 . , " -",666 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 667 . "ave(X)", zstophy,zout) 668 c 669 CALL histdef(nid_mth, "rhum", "Relative humidity", " -",665 . , "1", 666 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 667 . "ave(X)", zstophy,zout) 668 c 669 CALL histdef(nid_mth, "rhum", "Relative humidity", "1", 670 670 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 671 671 . "ave(X)", zstophy,zout) … … 917 917 . "once", zstophy,zout) 918 918 c 919 CALL histdef(nid_mth, "aire", "Grid area", " -",919 CALL histdef(nid_mth, "aire", "Grid area", "m2", 920 920 . iim,jj_nb,nhori, 1,1,1, -99, 32, 921 921 . "once", zstophy,zout) … … 1083 1083 . "ave(X)", zstorad,zout) 1084 1084 c 1085 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", " -",1085 CALL histdef(nid_mth, "sicf", "Sea-ice fraction", "1", 1086 1086 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1087 1087 . "ave(X)", zstophy,zout) … … 1115 1115 c 1116 1116 CALL histdef(nid_mth, "ndayrain", 1117 . "Number of day with rain (liq+sol)", " -",1117 . "Number of day with rain (liq+sol)", "days", 1118 1118 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1119 1119 . "inst(X)", zstomth,zout) … … 1406 1406 . "ave(X)", zstophy,zout) 1407 1407 c 1408 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", " -",1409 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1410 . "ave(X)", zstophy,zout) 1411 c 1412 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", " -",1413 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1414 . "ave(X)", zstophy,zout) 1415 c 1416 CALL histdef(nid_mth, "cldh", "High-level cloudiness", " -",1417 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1418 . "ave(X)", zstophy,zout) 1419 c 1420 CALL histdef(nid_mth, "cldt", "Total cloudiness", " %",1408 CALL histdef(nid_mth, "cldl", "Low-level cloudiness", "1", 1409 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1410 . "ave(X)", zstophy,zout) 1411 c 1412 CALL histdef(nid_mth, "cldm", "Mid-level cloudiness", "1", 1413 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1414 . "ave(X)", zstophy,zout) 1415 c 1416 CALL histdef(nid_mth, "cldh", "High-level cloudiness", "1", 1417 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1418 . "ave(X)", zstophy,zout) 1419 c 1420 CALL histdef(nid_mth, "cldt", "Total cloudiness", "1", 1421 1421 . iim,jj_nb,nhori, 1,1,1, -99, 32, 1422 1422 . "ave(X)", zstophy,zout) … … 1650 1650 . "ave(X)", zstophy,zout) 1651 1651 c 1652 CALL histdef(nid_mth, "rneb", "Cloud fraction", " %",1652 CALL histdef(nid_mth, "rneb", "Cloud fraction", "1", 1653 1653 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 1654 1654 . "ave(X)", zstophy,zout) 1655 1655 c 1656 1656 CALL histdef(nid_mth, "rnebcon", "Convective Cloud Fraction" 1657 . , " -",1658 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 1659 . "ave(X)", zstophy,zout) 1660 c 1661 CALL histdef(nid_mth, "rhum", "Relative humidity", " -",1657 . , "1", 1658 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 1659 . "ave(X)", zstophy,zout) 1660 c 1661 CALL histdef(nid_mth, "rhum", "Relative humidity", "1", 1662 1662 . iim,jj_nb,nhori, klev,1,klev,nvert, 32, 1663 1663 . "ave(X)", zstophy,zout) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/isccp_cloud_types.F
r776 r918 530 530 ! Initialised frac_out to zero 531 531 532 do i box=1,ncol533 do i lev=1,nlev532 do ilev=1,nlev 533 do ibox=1,ncol 534 534 do j=1,npoints 535 535 frac_out(j,ibox,ilev)=0.0 … … 1219 1219 enddo 1220 1220 do 29 ilev=1,nlev-1 1221 1221 !cdir nodep 1222 1222 do j=1,npoints 1223 1223 if ((at(j,ilev) .ge. tb(j,ibox) .and. -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/pbl_surface_mod.F90
r803 r918 962 962 patm(j)=ypplay(j,1) 963 963 qairsol(j) = yqsurf(j) 964 END DO 965 964 ! Blindage provisoire 965 rugo1(j)=MAX(rugo1(j),1.E-5) 966 END DO 966 967 967 968 ! Calculate the temperature et relative humidity at 2m and the wind at 10m -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/physiq.F
r898 r918 47 47 c CLEFS CPP POUR LES IO 48 48 c ===================== 49 c#define histhf49 #define histhf 50 50 #define histday 51 51 #define histmth 52 52 c#define histins 53 c#define histmthNMC54 c#define histISCCP53 #define histmthNMC 54 #define histISCCP 55 55 c====================================================================== 56 56 c modif ( P. Le Van , 12/10/98 ) … … 1487 1487 cym Attention pbase pas initialise dans concvl !!!! 1488 1488 pbase(:)=0 1489 1490 1489 first=.false. 1491 1490 endif … … 1675 1674 rugoro(i) = MAX(1.0e-05, zstd(i)*zsig(i)/2.0) 1676 1675 ENDDO 1676 1677 cFRED WRITE(lunout,*)'FRED rugoro ...' 1678 cFRED WRITE(lunout,*) rugoro (416), rugoro (417), rugoro(418) 1679 cFRED . , rugoro(419) 1680 cFRED WRITE(lunout,*) rugoro (421), rugoro (422), rugoro(423) 1681 cFRED WRITE(lunout,*) rugoro (424), rugoro (425), rugoro(426) 1682 cFRED WRITE(lunout,*) rugoro (427), rugoro (428), rugoro(429) 1683 cFRED WRITE(lunout,*) rugoro (430), rugoro (431), rugoro(567) 1684 cFRED . , rugoro(568) 1685 cFRED WRITE(lunout,*) rugoro (574), rugoro (575), rugoro(576) 1686 cFRED WRITE(lunout,*) rugoro (577), rugoro (698), rugoro(699) 1687 cFRED WRITE(lunout,*) rugoro (711), rugoro (712), rugoro(713) 1688 1689 WRITE(lunout,*)'FRED rugoro =0.' 1690 rugoro=0. 1691 1677 1692 CALL SUGWD(klon,klev,paprs,pplay) 1678 1693 DO i=1,klon … … 2594 2609 #ifdef histISCCP 2595 2610 IF (ok_isccp) THEN 2611 c SD initialisation 2612 fq_isccp(:,:,:,:)= 0. 2613 totalcldarea(:,:)= 0. 2614 meanptop(:,:) = 0. 2615 meantaucld(:,:) = 0. 2616 nbsunlit(:,:,:) = 0. 2617 c 2596 2618 cIM appel simulateur toutes les NINT(freq_ISCCP/dtime) heures 2597 2619 IF (MOD(itap,NINT(freq_ISCCP/dtime)).EQ.0) THEN 2598 2620 #include "calcul_simulISCCP.h" 2599 ENDIF !(MOD(itap,NINT(freq_ISCCP/dtime)) 2621 ENDIF !(MOD(itap,NINT(freq_ISCCP/dtime))) 2600 2622 ENDIF !ok_isccp 2601 2623 #endif -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/phytrac.F
r862 r918 445 445 c============================================================= 446 446 447 #ifdef CPP_IOIPSL448 #include "ini_histrac.h"449 #endif447 c#ifdef CPP_IOIPSL 448 c#include "ini_histrac.h" 449 c#endif 450 450 451 451 c====================================================================== … … 929 929 c============================================================= 930 930 931 #ifdef CPP_IOIPSL932 #include "write_histrac.h"933 #endif931 c#ifdef CPP_IOIPSL 932 c#include "write_histrac.h" 933 c#endif 934 934 935 935 c============================================================= -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/surf_landice_mod.F90
r793 r918 141 141 ! 142 142 !**************************************************************************************** 143 z0_new(:) = rugoro(:) 143 ! z0_new(:) = rugoro(:) 144 z0_new(:) = MAX(1.E-3,rugoro(:)) 144 145 145 146 -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histISCCP.h
r848 r918 12 12 c 13 13 DO n=1, napisccp 14 c 15 DO k=1,kmaxm1 16 zx_tmp_fi3d(1:klon, 1:lmaxm1)=fq_isccp(1:klon,k,1:lmaxm1,n)*100. 17 cym CALL gr_fi_ecrit(lmaxm1,klon,iim,jjmp1,zx_tmp_fi3d, 18 cym . zx_tmp_3d) 19 c 20 cIM: champ 3d : (lon,lat,pres) pour un tau fixe 21 c 14 DO k=1, kmaxm1 15 DO l=1, lmaxm1 16 IF(top_height.LE.2) THEN 17 DO i=1, klon 18 IF(nbsunlit(1,i,n).NE.0.) THEN 19 zx_tmp_fi3d(i,l)= 20 $ fq_isccp(i,k,l,n)*100./nbsunlit(1,i,n) 21 ELSE 22 zx_tmp_fi3d(i,l)=0. 23 ENDIF 24 ENDDO 25 ELSE IF(top_height.EQ.3) THEN 26 DO i=1, klon 27 zx_tmp_fi3d(i,l) = fq_isccp(i,k,l,n)*100. 28 ENDDO ! klon 29 ENDIF 30 ENDDO !l 22 31 CALL histwrite_phy(nid_isccp,"cldISCCP_"//taulev(k)//verticaxe(n) 23 32 . ,itau_w,zx_tmp_fi3d) 24 ENDDO !k 25 c 26 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1,nbsunlit(1,:,n),zx_tmp_2d) 27 CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),itau_w, 28 . nbsunlit(1,:,n)) 29 c 33 ENDDO !k 34 CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n),itau_w, 35 . nbsunlit(1,:,n)) 30 36 ENDDO ! n=1, napisccp 37 c 38 CALL histwrite_phy(nid_isccp,"ptop",itau_w,meanptop(:,1)) 39 c 40 CALL histwrite_phy(nid_isccp,"meantaucld",itau_w,meantaucld) 41 c 31 42 ELSE IF(type_run.EQ."AMIP".OR.type_run.EQ."CFMI") THEN 32 43 c … … 42 53 $ fq_isccp(i,k,l,n)*100./nbsunlit(1,i,n) 43 54 ELSE 44 fq_is_true(i,k,l,n)=0 55 fq_is_true(i,k,l,n)=0. 45 56 ENDIF 46 57 ENDDO … … 64 75 CALL histwrite_phy(nid_isccp,"nsunlit"//verticaxe(n), 65 76 . itau_w,nbsunlit(1,:,n)) 77 c 78 CALL histwrite_phy(nid_isccp,"ptop",itau_w,meanptop(:,1)) 79 c 80 CALL histwrite_phy(nid_isccp,"meantaucld",itau_w,meantaucld) 66 81 c 67 82 zx_tmp_fi2d(1:klon)=float(seed(1:klon,n)) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histday.h
r897 r918 342 342 CALL histwrite_phy(nid_day,"tauy",itau_w,zx_tmp_fi2d) 343 343 c 344 cIM IF (iflag_con.GE.3) THEN345 cIM ema_pct= 0. CALL gr_fi_ecrit(1, klon,iim,jjmp1, ema_pct,zx_tmp_2d)346 cym CALL gr_fi_ecrit(1, klon,iim,jjmp1, meanptop(:,1),zx_tmp_2d)347 IF (ok_isccp) THEN348 CALL histwrite_phy(nid_day,"ptop",itau_w,meanptop(:,1))349 ENDIF350 cIM ENDIF351 c352 344 ENDIF !lev_histday.GE.2 353 345 c … … 542 534 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, fiwc,zx_tmp_3d) 543 535 CALL histwrite_phy(nid_day,"iwcon",itau_w,fiwc) 544 c545 cym CALL gr_fi_ecrit(1,klon,iim,jjmp1, meantaucld,zx_tmp_2d)546 IF (ok_isccp) THEN547 CALL histwrite_phy(nid_day,"meantaucld",itau_w,meantaucld)548 ENDIF549 536 c 550 537 cym CALL gr_fi_ecrit(klev,klon,iim,jjmp1, cldtau,zx_tmp_3d) -
LMDZ4/branches/LMDZ4_V3_patches/libf/phylmd/write_histhf.h
r897 r918 67 67 DO k=1, nlevSTD 68 68 c 69 c bb=clevSTD(k) 70 c 71 c IF(k.GE.2) THEN 72 IF(k.GE.2.AND.k.LE.12) bb2=clevSTD(k) 73 IF(k.GE.13.AND.k.LE.17) bb3=clevSTD(k) 74 c aa=clevSTD(k) 75 c bb=aa(1:lnblnk1(aa)) 76 c ENDIF 77 c 69 bb2=clevSTD(k) 78 70 IF(bb2.EQ."500") THEN 79 71 c … … 113 105 c 114 106 DO k=1, nlevSTD 115 c116 c bb=clevSTD(k)117 c118 c IF(k.GE.2) THEN119 c IF(k.GE.1.AND.k.LE.10) bb2=clevSTD(k)120 c IF(k.GE.11.AND.k.LE.17) bb3=clevSTD(k)121 c aa=clevSTD(k)122 c bb=aa(1:lnblnk1(aa))123 c ENDIF124 107 c 125 108 bb2=clevSTD(k)
Note: See TracChangeset
for help on using the changeset viewer.
