Changeset 1664 for LMDZ5/branches/testing/libf/phylmd/radiation_AR4.F
- Timestamp:
- Oct 9, 2012, 3:29:15 PM (12 years ago)
- Location:
- LMDZ5/branches/testing
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ5/branches/testing
- Property svn:mergeinfo changed
/LMDZ5/trunk (added) merged: 1555-1557,1559-1566,1568-1575
- Property svn:mergeinfo changed
-
LMDZ5/branches/testing/libf/phylmd/radiation_AR4.F
r1279 r1664 19 19 cym#include "raddim.h" 20 20 #include "YOMCST.h" 21 #include "iniprint.h" 21 22 C 22 23 C ------------------------------------------------------------------ … … 163 164 allocate(ZFSUPAI(KDLON,KFLEV+1)) 164 165 allocate(ZFSDNAI(KDLON,KFLEV+1)) 165 DO JK = 1 , KDLON*(KFLEV+1) 166 ZFSUPAD(JK,1) = 0.0 ! ZFSUPAD(:,:)=0. 167 ZFSDNAD(JK,1) = 0.0 ! ZFSDNAD(:,:)=0. 168 ZFSUPAI(JK,1) = 0.0 ! ZFSUPAI(:,:)=0. 169 ZFSDNAI(JK,1) = 0.0 ! ZFSDNAI(:,:)=0. 170 END DO 166 167 ZFSUPAD(:,:)=0. 168 ZFSDNAD(:,:)=0. 169 ZFSUPAI(:,:)=0. 170 ZFSDNAI(:,:)=0. 171 171 endif 172 !rv 173 174 c 172 175 173 IF (appel1er) THEN 176 PRINT*,'SW calling frequency : ', swpas177 PRINT*," In general, it should be 1"174 WRITE(lunout,*) 'SW calling frequency : ', swpas 175 WRITE(lunout,*) " In general, it should be 1" 178 176 appel1er = .FALSE. 179 177 ENDIF … … 526 524 USE dimphy 527 525 USE radiation_AR4_param, only : RSUN, RRAY 526 USE infotrac, ONLY : type_trac 527 #ifdef REPROBUS 528 USE CHEM_REP, ONLY : RSUNTIME, ok_SUNTIME 529 #endif 530 528 531 IMPLICIT none 529 532 cym#include "dimensions.h" 530 533 cym#include "dimphy.h" 531 cym#include "raddim.h" 534 cym#include "raddim.h"i 535 #include "iniprint.h" 532 536 C 533 537 C ------------------------------------------------------------------ … … 613 617 INTEGER jl, jk, k, jaj, ikm1, ikl 614 618 619 C If running with Reporbus, overwrite default values of RSUN. 620 C Otherwise keep default values from radiation_AR4_param module. 621 IF (type_trac == 'repr') THEN 622 #ifdef REPROBUS 623 IF (ok_SUNTIME) THEN 624 RSUN(1) = RSUNTIME(1) 625 RSUN(2) = RSUNTIME(2) 626 ENDIF 627 WRITE(lunout,*)'RSUN(1): ',RSUN(1) 628 #endif 629 END IF 630 615 631 C ------------------------------------------------------------------ 616 632 C … … 754 770 USE dimphy 755 771 USE radiation_AR4_param, only : RSUN, RRAY 772 USE infotrac, ONLY : type_trac 773 #ifdef REPROBUS 774 use CHEM_REP, only : RSUNTIME, ok_SUNTIME 775 #endif 776 756 777 IMPLICIT none 757 778 cym#include "dimensions.h" … … 873 894 INTEGER jref, jkl, jklp1, jajp, jkki, jkkp4, jn2j, iabs 874 895 REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11 896 897 C If running with Reporbus, overwrite default values of RSUN. 898 C Otherwise keep default values from radiation_AR4_param module. 899 IF (type_trac == 'repr') THEN 900 #ifdef REPROBUS 901 IF (ok_SUNTIME) THEN 902 RSUN(1)=RSUNTIME(1) 903 RSUN(2)=RSUNTIME(2) 904 END IF 905 #endif 906 END IF 907 875 908 C 876 909 … … 2307 2340 #include "raddimlw.h" 2308 2341 #include "YOMCST.h" 2342 #include "iniprint.h" 2309 2343 C 2310 2344 C----------------------------------------------------------------------- … … 2412 2446 C ------------------------------------------------------------------ 2413 2447 IF (appel1er) THEN 2414 PRINT*,"LW clear-sky calling frequency: ", lw0pas2415 PRINT*,"LW cloudy-sky calling frequency: ", lwpas2416 PRINT*," In general, they should be 1"2448 WRITE(lunout,*) "LW clear-sky calling frequency: ", lw0pas 2449 WRITE(lunout,*) "LW cloudy-sky calling frequency: ", lwpas 2450 WRITE(lunout,*) " In general, they should be 1" 2417 2451 cym 2418 2452 allocate(ZFLUX(KDLON,2,KFLEV+1) ) … … 2485 2519 USE dimphy 2486 2520 USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT 2521 USE infotrac, ONLY : type_trac 2522 #ifdef REPROBUS 2523 USE CHEM_REP, ONLY: RCH42D, 2524 $ RN2O2D, 2525 $ RCFC112D, 2526 $ RCFC122D, 2527 $ ok_Rtime2D 2528 #endif 2529 2487 2530 IMPLICIT none 2488 2531 cym#include "dimensions.h" … … 2800 2843 S +ZUAER(JL,5) *ZDUC(JL,JC)*ZDIFF 2801 2844 C 2802 PABCU(JL,19,JC)=PABCU(JL,19,JCP1) 2803 S +ZABLY(JL,8,JC)*RCH4/RCO2*ZPHM6(JL)*ZDIFF 2804 PABCU(JL,20,JC)=PABCU(JL,20,JCP1) 2845 C 2846 2847 IF (type_trac == 'repr') THEN 2848 #ifdef REPROBUS 2849 IF (ok_Rtime2D) THEN 2850 PABCU(JL,19,JC)=PABCU(JL,19,JCP1) 2851 S +ZABLY(JL,8,JC)*RCH42D(JL,JC)/RCO2*ZPHM6(JL)*ZDIFF 2852 PABCU(JL,20,JC)=PABCU(JL,20,JCP1) 2853 S +ZABLY(JL,9,JC)*RCH42D(JL,JC)/RCO2*ZPSM6(JL)*ZDIFF 2854 PABCU(JL,21,JC)=PABCU(JL,21,JCP1) 2855 S +ZABLY(JL,8,JC)*RN2O2D(JL,JC)/RCO2*ZPHN6(JL)*ZDIFF 2856 PABCU(JL,22,JC)=PABCU(JL,22,JCP1) 2857 S +ZABLY(JL,9,JC)*RN2O2D(JL,JC)/RCO2*ZPSN6(JL)*ZDIFF 2858 C 2859 PABCU(JL,23,JC)=PABCU(JL,23,JCP1) 2860 S +ZABLY(JL,8,JC)*RCFC112D(JL,JC)/RCO2 *ZDIFF 2861 PABCU(JL,24,JC)=PABCU(JL,24,JCP1) 2862 S +ZABLY(JL,8,JC)*RCFC122D(JL,JC)/RCO2 *ZDIFF 2863 ELSE 2864 ! Same calculation as for type_trac /= repr 2865 PABCU(JL,19,JC)=PABCU(JL,19,JCP1) 2866 S +ZABLY(JL,8,JC)*RCH4/RCO2*ZPHM6(JL)*ZDIFF 2867 PABCU(JL,20,JC)=PABCU(JL,20,JCP1) 2868 S +ZABLY(JL,9,JC)*RCH4/RCO2*ZPSM6(JL)*ZDIFF 2869 PABCU(JL,21,JC)=PABCU(JL,21,JCP1) 2870 S +ZABLY(JL,8,JC)*RN2O/RCO2*ZPHN6(JL)*ZDIFF 2871 PABCU(JL,22,JC)=PABCU(JL,22,JCP1) 2872 S +ZABLY(JL,9,JC)*RN2O/RCO2*ZPSN6(JL)*ZDIFF 2873 C 2874 PABCU(JL,23,JC)=PABCU(JL,23,JCP1) 2875 S +ZABLY(JL,8,JC)*RCFC11/RCO2 *ZDIFF 2876 PABCU(JL,24,JC)=PABCU(JL,24,JCP1) 2877 S +ZABLY(JL,8,JC)*RCFC12/RCO2 *ZDIFF 2878 END IF 2879 #endif 2880 ELSE 2881 PABCU(JL,19,JC)=PABCU(JL,19,JCP1) 2882 S +ZABLY(JL,8,JC)*RCH4/RCO2*ZPHM6(JL)*ZDIFF 2883 PABCU(JL,20,JC)=PABCU(JL,20,JCP1) 2805 2884 S +ZABLY(JL,9,JC)*RCH4/RCO2*ZPSM6(JL)*ZDIFF 2806 PABCU(JL,21,JC)=PABCU(JL,21,JCP1) 2807 S +ZABLY(JL,8,JC)*RN2O/RCO2*ZPHN6(JL)*ZDIFF 2808 PABCU(JL,22,JC)=PABCU(JL,22,JCP1) 2809 S +ZABLY(JL,9,JC)*RN2O/RCO2*ZPSN6(JL)*ZDIFF 2810 C 2811 PABCU(JL,23,JC)=PABCU(JL,23,JCP1) 2812 S +ZABLY(JL,8,JC)*RCFC11/RCO2 *ZDIFF 2813 PABCU(JL,24,JC)=PABCU(JL,24,JCP1) 2814 S +ZABLY(JL,8,JC)*RCFC12/RCO2 *ZDIFF 2885 PABCU(JL,21,JC)=PABCU(JL,21,JCP1) 2886 S +ZABLY(JL,8,JC)*RN2O/RCO2*ZPHN6(JL)*ZDIFF 2887 PABCU(JL,22,JC)=PABCU(JL,22,JCP1) 2888 S +ZABLY(JL,9,JC)*RN2O/RCO2*ZPSN6(JL)*ZDIFF 2889 C 2890 PABCU(JL,23,JC)=PABCU(JL,23,JCP1) 2891 S +ZABLY(JL,8,JC)*RCFC11/RCO2 *ZDIFF 2892 PABCU(JL,24,JC)=PABCU(JL,24,JCP1) 2893 S +ZABLY(JL,8,JC)*RCFC12/RCO2 *ZDIFF 2894 END IF 2895 2815 2896 523 CONTINUE 2816 2897 524 CONTINUE
Note: See TracChangeset
for help on using the changeset viewer.