Changeset 282 for LMDZ.3.3/branches/rel-LF
- Timestamp:
- Oct 19, 2001, 12:32:41 PM (23 years ago)
- File:
-
- 1 edited
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (25 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F
r271 r282 244 244 REAL falbe(klon,nbsrf) 245 245 SAVE falbe ! albedo par type de surface 246 REAL falblw(klon,nbsrf) 247 SAVE falblw ! albedo par type de surface 248 246 249 c 247 250 c … … 366 369 REAL albsol(klon) 367 370 SAVE albsol ! albedo du sol total 371 REAL albsollw(klon) 372 SAVE albsollw ! albedo du sol total 373 368 374 REAL wo(klon,klev) 369 375 SAVE wo ! ozone … … 770 776 zout = dtime * FLOAT(ecrit_day) 771 777 C Essai writephys 772 773 774 775 776 777 778 779 780 781 778 c nom_fichier = 'histday1' 779 c call writephy_ini(fid_day,nom_fichier,klon,iim,jjmp1,klev, 780 c . rlon,rlat, presnivs, 781 c . zjulian, dtime) 782 c call writephy_def(prof2d_on, fid_day, "once", zsto, zout, 0) 783 c call writephy_def(prof3d_on, fid_day, "once", zsto, zout, 784 c . klev) 785 c call writephy_def(prof2d_av, fid_day, "ave(X)", zsto, zout, 0) 786 c call writephy_def(prof3d_av, fid_day, "ave(X)", zsto, zout, 787 c . klev) 782 788 783 789 c … … 1162 1168 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1163 1169 . "ave(X)", zsto,zout) 1170 CALL histdef(nid_mth, "albslw", "Surface albedo LW", "-", 1171 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1172 . "ave(X)", zsto,zout) 1164 1173 c 1165 1174 CALL histdef(nid_mth, "cdrm", "Momentum drag coef.", "-", … … 1589 1598 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1590 1599 . "inst(X)", zsto,zout) 1600 CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", 1601 . iim,jjmp1,nhori, 1,1,1, -99, 32, 1602 . "inst(X)", zsto,zout) 1591 1603 c 1592 1604 c … … 1795 1807 e ok_veget, ocean, npas, nexca, ftsol, 1796 1808 $ soil_model,ftsoil, 1797 $ paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,fluxlat, 1809 $ paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,falblw, 1810 $ fluxlat, 1798 1811 e rain_fall, snow_fall, solsw, sollw, sollwdown, fder, 1799 1812 e rlon, rlat, cufi, cvfi, frugs, … … 2173 2186 . + falbe(i,is_ter) * pctsrf(i,is_ter) 2174 2187 . + falbe(i,is_sic) * pctsrf(i,is_sic) 2188 albsollw(i) = falblw(i,is_oce) * pctsrf(i,is_oce) 2189 . + falblw(i,is_lic) * pctsrf(i,is_lic) 2190 . + falblw(i,is_ter) * pctsrf(i,is_ter) 2191 . + falblw(i,is_sic) * pctsrf(i,is_sic) 2175 2192 ENDDO 2176 2193 CALL radlwsw ! nouveau rayonnement (compatible Arpege-IFS) 2177 2194 e (dist, rmu0, fract, co2_ppm, solaire, 2178 e paprs, pplay,zxtsol,albsol, t_seri,q_seri,wo, 2195 e paprs, pplay,zxtsol,albsol, albsollw, t_seri,q_seri, 2196 e wo, 2179 2197 e cldfra, cldemi, cldtau, 2180 2198 s heat,heat0,cool,cool0,radsol,albpla, … … 2374 2392 vartitle= 'Surface geop. height' 2375 2393 varunits= '-' 2376 2377 2394 c call writephy(fid_day,prof2d_on,varname,pphis,vartitle, 2395 c . varunits) 2378 2396 c 2379 2397 i = NINT(zout/zsto) … … 2383 2401 vartitle= 'Grid area' 2384 2402 varunits= '-' 2385 2386 2403 c call writephy(fid_day,prof2d_on,varname,paire,vartitle, 2404 c . varunits) 2387 2405 C 2388 2406 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2389 2407 CALL histwrite(nid_day,"tsol",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2390 2391 2408 c call writephy(fid_day,prof2d_av,'tsol',zxtsol, 2409 c . 'Surface Temperature','K') 2392 2410 c 2393 2411 C … … 2395 2413 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d ,zx_tmp_2d) 2396 2414 CALL histwrite(nid_day,"tter",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2397 2398 2415 c call writephy(fid_day,prof2d_av,'tter',ftsol(1 : klon, is_ter), 2416 c . 'Surface Temperature','K') 2399 2417 C 2400 2418 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_lic) 2401 2419 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2402 2420 CALL histwrite(nid_day,"tlic",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2403 2404 2421 c call writephy(fid_day,prof2d_av,'tlic',ftsol(1 : klon, is_lic), 2422 c . 'Surface Temperature','K') 2405 2423 C 2406 2424 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_oce) 2407 2425 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2408 2426 CALL histwrite(nid_day,"toce",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2409 2410 2427 c call writephy(fid_day,prof2d_av,'toce',ftsol(1 : klon, is_oce), 2428 c . 'Surface Temperature','K') 2411 2429 C 2412 2430 zx_tmp_fi2d(1 : klon) = ftsol(1 : klon, is_sic) 2413 2431 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxtsol,zx_tmp_2d) 2414 2432 CALL histwrite(nid_day,"tsic",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2415 2416 2433 c call writephy(fid_day,prof2d_av,'tsic',ftsol(1 : klon, is_sic), 2434 c . 'Surface Temperature','K') 2417 2435 C 2418 2436 DO i = 1, klon … … 2425 2443 vartitle= 'pression au sol' 2426 2444 varunits= 'hPa' 2427 2428 2445 c call writephy(fid_day,prof2d_av,varname,zx_tmp_fi2d,vartitle, 2446 c . varunits) 2429 2447 c 2430 2448 DO i = 1, klon … … 2433 2451 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zx_tmp_fi2d,zx_tmp_2d) 2434 2452 CALL histwrite(nid_day,"rain",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2435 2436 2453 c call writephy(fid_day,prof2d_av,'rain',zx_tmp_fi2d, 2454 c . 'Precipitation','mm/day') 2437 2455 2438 2456 … … 2440 2458 CALL gr_fi_ecrit(1, klon,iim,jjmp1, snow_fall,zx_tmp_2d) 2441 2459 CALL histwrite(nid_day,"snow",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2442 2443 2460 c call writephy(fid_day,prof2d_av,'snow',snow_fall, 2461 c . 'Snow','mm/day') 2444 2462 c 2445 2463 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d) 2446 2464 CALL histwrite(nid_day,"snow_cov",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2447 2448 2465 c call writephy(fid_day,prof2d_av,'snow_cov',zxsnow, 2466 c . 'Snow cover','mm') 2449 2467 c 2450 2468 CALL gr_fi_ecrit(1, klon,iim,jjmp1, evap,zx_tmp_2d) 2451 2469 CALL histwrite(nid_day,"evap",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2452 2453 2470 c call writephy(fid_day,prof2d_av,'evap',evap, 2471 c . 'Evaporation','mm/day') 2454 2472 c 2455 2473 CALL gr_fi_ecrit(1, klon,iim,jjmp1, topsw,zx_tmp_2d) 2456 2474 CALL histwrite(nid_day,"tops",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2457 2458 2475 c call writephy(fid_day,prof2d_av,'tops',topsw, 2476 c . 'Solar rad. at TOA','W/m2') 2459 2477 c 2460 2478 CALL gr_fi_ecrit(1, klon,iim,jjmp1, toplw,zx_tmp_2d) 2461 2479 CALL histwrite(nid_day,"topl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2462 2463 2480 c call writephy(fid_day,prof2d_av,'topl',toplw, 2481 c . 'IR rad. at TOA','W/m2') 2464 2482 c 2465 2483 CALL gr_fi_ecrit(1, klon,iim,jjmp1, solsw,zx_tmp_2d) 2466 2484 CALL histwrite(nid_day,"sols",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2467 2468 2485 c call writephy(fid_day,prof2d_av,'sols',solsw, 2486 c . 'Solar rad. at surf.','W/m2') 2469 2487 c 2470 2488 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollw,zx_tmp_2d) 2471 2489 CALL histwrite(nid_day,"soll",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2472 2473 2490 c call writephy(fid_day,prof2d_av,'soll',sollw, 2491 c . 'IR rad. at surface','W/m2') 2474 2492 c 2475 2493 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sollwdown,zx_tmp_2d) 2476 2494 CALL histwrite(nid_day,"solldown",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2477 2478 2495 c call writephy(fid_day,prof2d_av,'solldown',sollwdown, 2496 c . 'Down. IR rad. at surface','W/m2') 2479 2497 c 2480 2498 CALL gr_fi_ecrit(1, klon,iim,jjmp1, bils,zx_tmp_2d) 2481 2499 CALL histwrite(nid_day,"bils",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2482 2483 2500 c call writephy(fid_day,prof2d_av,'bils',bils, 2501 c . 'Surf. total heat flux','W/m2') 2484 2502 c 2485 2503 CALL gr_fi_ecrit(1, klon,iim,jjmp1, sens,zx_tmp_2d) 2486 2504 CALL histwrite(nid_day,"sens",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2487 2488 2505 c call writephy(fid_day,prof2d_av,'sens',sens, 2506 c . 'Sensible heat flux','W/m2') 2489 2507 c 2490 2508 CALL gr_fi_ecrit(1, klon,iim,jjmp1, fder,zx_tmp_2d) 2491 2509 CALL histwrite(nid_day,"fder",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2492 2493 2510 c call writephy(fid_day,prof2d_av,'fder',fder, 2511 c . 'Heat flux derivation','W/m2') 2494 2512 c 2495 2513 c … … 2500 2518 CALL histwrite(nid_day,"pourc_"//clnsurf(nsrf),itap, 2501 2519 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2502 2503 2504 2520 c call writephy(fid_day,prof2d_av,'pourc_'//clnsurf(nsrf), 2521 c . pctsrf( 1 : klon, nsrf), 2522 c . 'Fraction'//clnsurf(nsrf),'-') 2505 2523 C 2506 2524 zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf) … … 2508 2526 CALL histwrite(nid_day,"tsol_"//clnsurf(nsrf),itap, 2509 2527 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2510 2511 2512 2528 c call writephy(fid_day,prof2d_av,'tsol_'//clnsurf(nsrf), 2529 c . ftsol( 1 : klon, nsrf), 2530 c . 'Surf. Temp'//clnsurf(nsrf),'K') 2513 2531 C 2514 2532 zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf) … … 2516 2534 CALL histwrite(nid_day,"sens_"//clnsurf(nsrf),itap, 2517 2535 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2518 2519 2520 2536 c call writephy(fid_day,prof2d_av,'sens_'//clnsurf(nsrf), 2537 c . fluxt( 1 : klon, 1, nsrf), 2538 c . 'Sensible heat flux '//clnsurf(nsrf),'W/m2') 2521 2539 C 2522 2540 zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf) … … 2524 2542 CALL histwrite(nid_day,"lat_"//clnsurf(nsrf),itap, 2525 2543 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2526 2527 2528 2544 c call writephy(fid_day,prof2d_av,'lat_'//clnsurf(nsrf), 2545 c . fluxlat( 1 : klon, nsrf), 2546 c . 'Latent heat flux '//clnsurf(nsrf),'W/m2') 2529 2547 C 2530 2548 zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf) … … 2532 2550 CALL histwrite(nid_day,"taux_"//clnsurf(nsrf),itap, 2533 2551 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2534 2535 2536 2552 c call writephy(fid_day,prof2d_av,'taux_'//clnsurf(nsrf), 2553 c . fluxu( 1 : klon, 1, nsrf), 2554 c . 'Zonal wind stress '//clnsurf(nsrf),'Pa') 2537 2555 C 2538 2556 zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf) … … 2540 2558 CALL histwrite(nid_day,"tauy_"//clnsurf(nsrf),itap, 2541 2559 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2542 2543 2544 2560 c call writephy(fid_day,prof2d_av,'tauy_'//clnsurf(nsrf), 2561 c . fluxv( 1 : klon, 1, nsrf), 2562 c . 'Meridional wind stress '//clnsurf(nsrf),'Pa') 2545 2563 C 2546 2564 zx_tmp_fi2d(1 : klon) = falbe( 1 : klon, nsrf) … … 2548 2566 CALL histwrite(nid_day,"albe_"//clnsurf(nsrf),itap, 2549 2567 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2550 call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf), 2551 . falbe( 1 : klon, nsrf), 2552 . 'Albedo surf.'//clnsurf(nsrf),'-') 2568 c call writephy(fid_day,prof2d_av,'albe_'//clnsurf(nsrf), 2569 c . falbe( 1 : klon, nsrf), 2570 c . 'Albedo surf. SW'//clnsurf(nsrf),'-') 2571 c call writephy(fid_day,prof2d_av,'alblw_'//clnsurf(nsrf), 2572 c . falblw( 1 : klon, nsrf), 2573 c . 'Albedo surf. LW'//clnsurf(nsrf),'-') 2553 2574 C 2554 2575 zx_tmp_fi2d(1 : klon) = frugs( 1 : klon, nsrf) … … 2556 2577 CALL histwrite(nid_day,"rugs_"//clnsurf(nsrf),itap, 2557 2578 $ zx_tmp_2d,iim*jjmp1,ndex2d) 2558 2559 2560 2579 c call writephy(fid_day,prof2d_av,'rugs_'//clnsurf(nsrf), 2580 c . frugs( 1 : klon, nsrf), 2581 c . 'Rugosity '//clnsurf(nsrf),' - ') 2561 2582 C 2562 2583 END DO … … 2570 2591 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldl,zx_tmp_2d) 2571 2592 CALL histwrite(nid_day,"cldl",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2572 2573 2593 c call writephy(fid_day,prof2d_av,'cldl',cldl, 2594 c . 'Low-level cloudiness','-') 2574 2595 c 2575 2596 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldm,zx_tmp_2d) 2576 2597 CALL histwrite(nid_day,"cldm",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2577 2578 2598 c call writephy(fid_day,prof2d_av,'cldm',cldm, 2599 c . 'Mid-level cloudiness','-') 2579 2600 c 2580 2601 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldh,zx_tmp_2d) 2581 2602 CALL histwrite(nid_day,"cldh",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2582 2583 2603 c call writephy(fid_day,prof2d_av,'cldh',cldh, 2604 c . 'High-level cloudiness','-') 2584 2605 c 2585 2606 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldt,zx_tmp_2d) 2586 2607 CALL histwrite(nid_day,"cldt",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2587 2588 2608 c call writephy(fid_day,prof2d_av,'cldt',cldt, 2609 c . 'Total cloudiness','-') 2589 2610 c 2590 2611 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cldq,zx_tmp_2d) 2591 2612 CALL histwrite(nid_day,"cldq",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2592 2593 2613 c call writephy(fid_day,prof2d_av,'cldq',cldq, 2614 c . 'Cloud liquid water path','-') 2594 2615 c 2595 2616 c Champs 3D: … … 2602 2623 vartitle= 'temperature 3D' 2603 2624 varunits= 'K' 2604 2625 c call writephy(fid_day,prof3d_av,varname,t_seri,vartitle,varunits) 2605 2626 c 2606 2627 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, qx(1,1,ivap), zx_tmp_3d) 2607 2628 CALL histwrite(nid_day,"ovap",itap,zx_tmp_3d, 2608 2629 . iim*jjmp1*klev,ndex3d) 2609 2610 2630 c call writephy(fid_day,prof3d_av,'ovap',qx(1,1,ivap), 2631 c . 'Specific humidity','Kg/Kg') 2611 2632 c 2612 2633 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, zphi, zx_tmp_3d) 2613 2634 CALL histwrite(nid_day,"geop",itap,zx_tmp_3d, 2614 2635 . iim*jjmp1*klev,ndex3d) 2615 2616 2636 c call writephy(fid_day,prof3d_av,'geop',zphi, 2637 c . 'Geopotential height','m') 2617 2638 c 2618 2639 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, u_seri, zx_tmp_3d) 2619 2640 CALL histwrite(nid_day,"vitu",itap,zx_tmp_3d, 2620 2641 . iim*jjmp1*klev,ndex3d) 2621 2622 2642 c call writephy(fid_day,prof3d_av,'vitu',u_seri, 2643 c . 'Zonal wind','m/s') 2623 2644 c 2624 2645 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, v_seri, zx_tmp_3d) 2625 2646 CALL histwrite(nid_day,"vitv",itap,zx_tmp_3d, 2626 2647 . iim*jjmp1*klev,ndex3d) 2627 2628 2648 c call writephy(fid_day,prof3d_av,'vitv',v_seri, 2649 c . 'Meridional wind','m/s') 2629 2650 c 2630 2651 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, omega, zx_tmp_3d) 2631 2652 CALL histwrite(nid_day,"vitw",itap,zx_tmp_3d, 2632 2653 . iim*jjmp1*klev,ndex3d) 2633 2634 2654 c call writephy(fid_day,prof3d_av,'vitw',omega, 2655 c . 'Vertical wind','m/s') 2635 2656 c 2636 2657 CALL gr_fi_ecrit(klev,klon,iim,jjmp1, pplay, zx_tmp_3d) 2637 2658 CALL histwrite(nid_day,"pres",itap,zx_tmp_3d, 2638 2659 . iim*jjmp1*klev,ndex3d) 2639 2640 2660 c call writephy(fid_day,prof3d_av,'pres',pplay, 2661 c . 'Air pressure','Pa') 2641 2662 2642 2663 c 2643 2664 if (ok_sync) then 2644 2665 c call writephy_sync(fid_day) 2645 2666 call histsync(nid_day) 2646 2667 endif … … 2809 2830 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 2810 2831 CALL histwrite(nid_mth,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2832 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) 2833 CALL histwrite(nid_mth,"albslw",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 2811 2834 c 2812 2835 CALL gr_fi_ecrit(1, klon,iim,jjmp1, cdragm,zx_tmp_2d) … … 3179 3202 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsol,zx_tmp_2d) 3180 3203 CALL histwrite(nid_ins,"albs",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3204 CALL gr_fi_ecrit(1, klon,iim,jjmp1, albsollw,zx_tmp_2d) 3205 CALL histwrite(nid_ins,"albslw",itap,zx_tmp_2d,iim*jjmp1,ndex2d) 3181 3206 c 3182 3207 CALL gr_fi_ecrit(1, klon,iim,jjmp1, zxsnow,zx_tmp_2d)
Note: See TracChangeset
for help on using the changeset viewer.
