Changeset 634 for LMDZ4/branches/LMDZ4_par_0/libf/phylmd/radlwsw.F
- Timestamp:
- May 4, 2005, 5:11:29 PM (19 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4_par_0/libf/phylmd/radlwsw.F
r633 r634 16 16 . cldtaupi, topswai, solswai) 17 17 c 18 USE dimphy 18 19 IMPLICIT none 19 20 c====================================================================== … … 69 70 70 71 c====================================================================== 71 #include "dimensions.h"72 #include "dimphy.h"73 #include "raddim.h"72 cym#include "dimensions.h" 73 cym#include "dimphy.h" 74 cym#include "raddim.h" 74 75 #include "YOETHF.h" 75 76 c … … 428 429 S PTOPSWAD,PSOLSWAD,PTOPSWAI,PSOLSWAI, 429 430 J ok_ade, ok_aie ) 430 431 USE dimphy 431 432 IMPLICIT none 432 433 433 #include "dimensions.h"434 #include "dimphy.h"435 #include "raddim.h"434 cym#include "dimensions.h" 435 cym#include "dimphy.h" 436 cym#include "raddim.h" 436 437 #include "YOMCST.h" 437 438 C … … 547 548 REAL*8 PSOLSWAI(KDLON) ! SHORTWAVE FLUX AT SURFACE(+AEROSOL IND) 548 549 cjq - Fluxes including aerosol effects 549 REAL*8 ZFSUPAD(KDLON,KFLEV+1)550 REAL*8 ZFSDNAD(KDLON,KFLEV+1)551 REAL*8 ZFSUPAI(KDLON,KFLEV+1)552 REAL*8 ZFSDNAI(KDLON,KFLEV+1)550 REAL*8,allocatable,save :: ZFSUPAD(:,:) 551 REAL*8,allocatable,save :: ZFSDNAD(:,:) 552 REAL*8,allocatable,save :: ZFSUPAI(:,:) 553 REAL*8,allocatable,save :: ZFSDNAI(:,:) 553 554 logical initialized 554 SAVE ZFSUPAD, ZFSDNAD, ZFSUPAI, ZFSDNAI ! aerosol fluxes555 cym SAVE ZFSUPAD, ZFSDNAD, ZFSUPAI, ZFSDNAI ! aerosol fluxes 555 556 !rv 556 557 save flag_aer … … 560 561 flag_aer=0. 561 562 initialized=.TRUE. 563 allocate(ZFSUPAD(KDLON,KFLEV+1)) 564 allocate(ZFSDNAD(KDLON,KFLEV+1)) 565 allocate(ZFSUPAI(KDLON,KFLEV+1)) 566 allocate(ZFSDNAI(KDLON,KFLEV+1)) 567 ZFSUPAD(:,:)=0. 568 ZFSDNAD(:,:)=0. 569 ZFSUPAI(:,:)=0. 570 ZFSDNAI(:,:)=0. 571 562 572 endif 563 573 !rv … … 735 745 S PTAVE,PWV,PAKI,PCLD,PCLEAR,PDSIG,PFACT, 736 746 S PRMU,PSEC,PUD) 747 USE dimphy 737 748 IMPLICIT none 738 #include "dimensions.h"739 #include "dimphy.h"740 #include "raddim.h"749 cym#include "dimensions.h" 750 cym#include "dimphy.h" 751 cym#include "raddim.h" 741 752 #include "radepsi.h" 742 753 #include "radopt.h" … … 925 936 S , PDSIG , POMEGA, POZ , PRMU , PSEC , PTAU , PUD 926 937 S , PFD , PFU) 938 USE dimphy 927 939 IMPLICIT none 928 #include "dimensions.h"929 #include "dimphy.h"930 #include "raddim.h"940 cym#include "dimensions.h" 941 cym#include "dimphy.h" 942 cym#include "raddim.h" 931 943 C 932 944 C ------------------------------------------------------------------ … … 1165 1177 S , PUD ,PWV , PQS 1166 1178 S , PFDOWN,PFUP ) 1179 USE dimphy 1167 1180 IMPLICIT none 1168 #include "dimensions.h"1169 #include "dimphy.h"1170 #include "raddim.h"1181 cym#include "dimensions.h" 1182 cym#include "dimphy.h" 1183 cym#include "raddim.h" 1171 1184 #include "radepsi.h" 1172 1185 C … … 1711 1724 S , PCGAZ , PPIZAZ, PRAY1 , PRAY2 , PREFZ , PRJ 1712 1725 S , PRK , PRMU0 , PTAUAZ, PTRA1 , PTRA2 ) 1726 USE dimphy 1713 1727 IMPLICIT none 1714 #include "dimensions.h"1715 #include "dimphy.h"1716 #include "raddim.h"1728 cym#include "dimensions.h" 1729 cym#include "dimphy.h" 1730 cym#include "raddim.h" 1717 1731 #include "radepsi.h" 1718 1732 #include "radopt.h" … … 2071 2085 S , PCGAZ , PPIZAZ, PRAY1, PRAY2, PREFZ , PRJ , PRK , PRMUE 2072 2086 S , PTAUAZ, PTRA1 , PTRA2 ) 2087 USE dimphy 2073 2088 IMPLICIT none 2074 #include "dimensions.h"2075 #include "dimphy.h"2076 #include "raddim.h"2089 cym#include "dimensions.h" 2090 cym#include "dimphy.h" 2091 cym#include "raddim.h" 2077 2092 #include "radepsi.h" 2078 2093 #include "radopt.h" … … 2427 2442 SUBROUTINE SWDE (PGG,PREF,PRMUZ,PTO1,PW, 2428 2443 S PRE1,PRE2,PTR1,PTR2) 2444 USE dimphy 2429 2445 IMPLICIT none 2430 #include "dimensions.h"2431 #include "dimphy.h"2432 #include "raddim.h"2446 cym#include "dimensions.h" 2447 cym#include "dimphy.h" 2448 cym#include "raddim.h" 2433 2449 C 2434 2450 C ------------------------------------------------------------------ … … 2558 2574 END 2559 2575 SUBROUTINE SWTT (KNU,KA,PU,PTR) 2576 USE dimphy 2560 2577 IMPLICIT none 2561 #include "dimensions.h"2562 #include "dimphy.h"2563 #include "raddim.h"2578 cym#include "dimensions.h" 2579 cym#include "dimphy.h" 2580 cym#include "raddim.h" 2564 2581 C 2565 2582 C----------------------------------------------------------------------- … … 2674 2691 END 2675 2692 SUBROUTINE SWTT1(KNU,KABS,KIND, PU, PTR) 2693 USE dimphy 2676 2694 IMPLICIT none 2677 #include "dimensions.h"2678 #include "dimphy.h"2679 #include "raddim.h"2695 cym#include "dimensions.h" 2696 cym#include "dimphy.h" 2697 cym#include "raddim.h" 2680 2698 C 2681 2699 C----------------------------------------------------------------------- … … 2806 2824 cIM . ptoplwdown,ptoplwdownclr) 2807 2825 . plwup, plwdn, plwup0, plwdn0) 2826 USE dimphy 2808 2827 IMPLICIT none 2809 #include "dimensions.h"2810 #include "dimphy.h"2811 #include "raddim.h"2828 cym#include "dimensions.h" 2829 cym#include "dimphy.h" 2830 cym#include "raddim.h" 2812 2831 #include "raddimlw.h" 2813 2832 #include "YOMCST.h" … … 2884 2903 REAL*8 ZOZ(KDLON,KFLEV) 2885 2904 c 2886 REAL*8 ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down) 2887 REAL*8 ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES 2888 REAL*8 ZBINT(KDLON,KFLEV+1) ! Intermediate variable 2889 REAL*8 ZBSUI(KDLON) ! Intermediate variable 2890 REAL*8 ZCTS(KDLON,KFLEV) ! Intermediate variable 2891 REAL*8 ZCNTRB(KDLON,KFLEV+1,KFLEV+1) ! Intermediate variable 2892 SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB 2905 cym REAL*8 ZFLUX(KDLON,2,KFLEV+1) ! RADIATIVE FLUXES (1:up; 2:down) 2906 cym REAL*8 ZFLUC(KDLON,2,KFLEV+1) ! CLEAR-SKY RADIATIVE FLUXES 2907 cym REAL*8 ZBINT(KDLON,KFLEV+1) ! Intermediate variable 2908 cym REAL*8 ZBSUI(KDLON) ! Intermediate variable 2909 cym REAL*8,ZCTS(KDLON,KFLEV) ! Intermediate variable 2910 cym REAL*8 ZCNTRB(KDLON,KFLEV+1,KFLEV+1) ! Intermediate variable 2911 cym SAVE ZFLUX, ZFLUC, ZBINT, ZBSUI, ZCTS, ZCNTRB 2912 REAL*8,allocatable,save :: ZFLUX(:,:,:) ! RADIATIVE FLUXES (1:up; 2:down) 2913 REAL*8,allocatable,save :: ZFLUC(:,:,:) ! CLEAR-SKY RADIATIVE FLUXES 2914 REAL*8,allocatable,save :: ZBINT(:,:) ! Intermediate variable 2915 REAL*8,allocatable,save :: ZBSUI(:) ! Intermediate variable 2916 REAL*8,allocatable,save :: ZCTS(:,:) ! Intermediate variable 2917 REAL*8,allocatable,save :: ZCNTRB(:,:,:) ! Intermediate variable 2893 2918 c 2894 2919 INTEGER ilim, i, k, kpl1 … … 2909 2934 PRINT*, "LW cloudy-sky calling frequency: ", lwpas 2910 2935 PRINT*, " In general, they should be 1" 2936 cym 2937 allocate(ZFLUX(KDLON,2,KFLEV+1) ) 2938 allocate(ZFLUC(KDLON,2,KFLEV+1) ) 2939 allocate(ZBINT(KDLON,KFLEV+1)) 2940 allocate(ZBSUI(KDLON)) 2941 allocate(ZCTS(KDLON,KFLEV)) 2942 allocate(ZCNTRB(KDLON,KFLEV+1,KFLEV+1)) 2911 2943 appel1er=.FALSE. 2912 2944 ENDIF … … 2973 3005 S PAER,PDP,PPMB,PPSOL,POZ,PTAVE,PVIEW,PWV, 2974 3006 S PABCU) 3007 USE dimphy 2975 3008 IMPLICIT none 2976 #include "dimensions.h"2977 #include "dimphy.h"2978 #include "raddim.h"3009 cym#include "dimensions.h" 3010 cym#include "dimphy.h" 3011 cym#include "raddim.h" 2979 3012 #include "raddimlw.h" 2980 3013 #include "YOMCST.h" … … 3366 3399 SUBROUTINE LWBV(KLIM,PDP,PDT0,PEMIS,PPMB,PTL,PTAVE,PABCU, 3367 3400 S PFLUC,PBINT,PBSUI,PCTS,PCNTRB) 3401 USE dimphy 3368 3402 IMPLICIT none 3369 #include "dimensions.h"3370 #include "dimphy.h"3371 #include "raddim.h"3403 cym#include "dimensions.h" 3404 cym#include "dimphy.h" 3405 cym#include "raddim.h" 3372 3406 #include "raddimlw.h" 3373 3407 #include "YOMCST.h" … … 3457 3491 R PBINT,PBSUIN,PCTS,PCNTRB, 3458 3492 S PFLUX) 3493 USE dimphy 3459 3494 IMPLICIT none 3460 #include "dimensions.h"3461 #include "dimphy.h"3462 #include "raddim.h"3495 cym#include "dimensions.h" 3496 cym#include "dimphy.h" 3497 cym#include "raddim.h" 3463 3498 #include "radepsi.h" 3464 3499 #include "radopt.h" … … 3862 3897 S , PB,PBINT,PBSUIN,PBSUR,PBTOP,PDBSL 3863 3898 S , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP) 3899 USE dimphy 3864 3900 IMPLICIT none 3865 #include "dimensions.h"3866 #include "dimphy.h"3867 #include "raddim.h"3901 cym#include "dimensions.h" 3902 cym#include "dimphy.h" 3903 cym#include "raddim.h" 3868 3904 #include "raddimlw.h" 3869 3905 C … … 5273 5309 R , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP 5274 5310 S , PCNTRB,PCTS,PFLUC) 5311 USE dimphy 5275 5312 IMPLICIT none 5276 #include "dimensions.h"5277 #include "dimphy.h"5278 #include "raddim.h"5313 cym#include "dimensions.h" 5314 cym#include "dimphy.h" 5315 cym#include "raddim.h" 5279 5316 #include "raddimlw.h" 5280 5317 #include "YOMCST.h" … … 5386 5423 R , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP 5387 5424 S , PCTS,PFLUC) 5425 USE dimphy 5388 5426 IMPLICIT none 5389 #include "dimensions.h"5390 #include "dimphy.h"5391 #include "raddim.h"5427 cym#include "dimensions.h" 5428 cym#include "dimphy.h" 5429 cym#include "raddim.h" 5392 5430 #include "raddimlw.h" 5393 5431 #include "radopt.h" … … 5718 5756 R , PGA,PGB 5719 5757 S , PCNTRB,PDISD,PDISU) 5758 USE dimphy 5720 5759 IMPLICIT none 5721 #include "dimensions.h"5722 #include "dimphy.h"5723 #include "raddim.h"5760 cym#include "dimensions.h" 5761 cym#include "dimphy.h" 5762 cym#include "raddim.h" 5724 5763 #include "raddimlw.h" 5725 5764 C … … 5970 6009 R , PABCU,PDBSL,PGA,PGB 5971 6010 S , PADJD,PADJU,PCNTRB,PDBDT) 6011 USE dimphy 5972 6012 IMPLICIT none 5973 #include "dimensions.h"5974 #include "dimphy.h"5975 #include "raddim.h"6013 cym#include "dimensions.h" 6014 cym#include "dimphy.h" 6015 cym#include "raddim.h" 5976 6016 #include "raddimlw.h" 5977 6017 C … … 6177 6217 END 6178 6218 SUBROUTINE LWTT(PGA,PGB,PUU, PTT) 6219 USE dimphy 6179 6220 IMPLICIT none 6180 #include "dimensions.h"6181 #include "dimphy.h"6182 #include "raddim.h"6221 cym#include "dimensions.h" 6222 cym#include "dimphy.h" 6223 cym#include "raddim.h" 6183 6224 #include "raddimlw.h" 6184 6225 C … … 6352 6393 END 6353 6394 SUBROUTINE LWTTM(PGA,PGB,PUU1,PUU2, PTT) 6395 USE dimphy 6354 6396 IMPLICIT none 6355 #include "dimensions.h"6356 #include "dimphy.h"6357 #include "raddim.h"6397 cym#include "dimensions.h" 6398 cym#include "dimphy.h" 6399 cym#include "raddim.h" 6358 6400 #include "raddimlw.h" 6359 6401 C
Note: See TracChangeset
for help on using the changeset viewer.