Index: /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_BRODIE.opt =================================================================== --- /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_BRODIE.opt (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_BRODIE.opt (revision 1250) @@ -0,0 +1,12 @@ +%INLINE -pi auto exp=swtt1_lmdar4,swtt_lmdar4,swde_lmdar4,lwttm_lmdar4,lwtt_lmdar4,swr_lmdar4,swclr_lmdar4 noexp=SW_LMDAR4,SWU_LMDAR4,SW1S_LMDAR4,SW2S_LMDAR4,LW_LMDAR4,LWU_LMDAR4,LWBV_LMDAR4,LWC_LMDAR4,LWB_LMDAR4,LWV_LMDAR4,LWVB_LMDAR4,LWVD_LMDAR4,LWVN_LMDAR4 line=2000 + +bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto +bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR +bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -Wf,-O,extendreorder %INLINE +bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt %INLINE +bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder +bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt Index: /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_MERCURE.opt =================================================================== --- /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_MERCURE.opt (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/arch/arch-SX8_MERCURE.opt (revision 1250) @@ -0,0 +1,12 @@ +%INLINE -pi auto exp=swtt1_lmdar4,swtt_lmdar4,swde_lmdar4,lwttm_lmdar4,lwtt_lmdar4,swr_lmdar4,swclr_lmdar4 noexp=SW_LMDAR4,SWU_LMDAR4,SW1S_LMDAR4,SW2S_LMDAR4,LW_LMDAR4,LWU_LMDAR4,LWBV_LMDAR4,LWC_LMDAR4,LWB_LMDAR4,LWV_LMDAR4,LWVB_LMDAR4,LWVD_LMDAR4,LWVN_LMDAR4 line=2000 + +bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto +bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR +bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -Wf,-O,extendreorder %INLINE +bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt %INLINE +bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder +bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt Index: /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.fcm =================================================================== --- /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.fcm (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.fcm (revision 1250) @@ -0,0 +1,15 @@ +%COMPILER sxmpif90 +%LINK sxmpif90 +%AR sxar +%MAKE make +%FPP_FLAGS -P -traditional +%FPP_DEF NC_DOUBLE BLAS SGEMV=DGEMV SGEMM=DGEMM FFT_MATHKEISAN +%BASE_FFLAGS -P stack -Wf,-pvctl res=whole,-A dbl4,-ptr byte -EP -R2 -float0 -size_t64 -dw -Wf,"-pvctl loopcnt=999999 fullmsg noassume" +%PROD_FFLAGS -C vopt -pi expin=%SRC_PATH/%DYN/cray.F exp=ssum,scopy +%DEV_FFLAGS -C vsafe -gv -Wf,-init stack=nan,-init heap=nan +%DEBUG_FFLAGS -C debug -eC -Wf,-init stack=nan,-init heap=nan +%MPI_FFLAGS +%OMP_FFLAGS -P openmp +%BASE_LD -size_t64 -lblas -lfft +%MPI_LD +%OMP_LD -P openmp -Wl,"-ZL 3G" Index: /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.opt =================================================================== --- /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.opt (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.opt (revision 1250) @@ -0,0 +1,12 @@ +%INLINE -pi auto exp=swtt1_lmdar4,swtt_lmdar4,swde_lmdar4,lwttm_lmdar4,lwtt_lmdar4,swr_lmdar4,swclr_lmdar4 noexp=SW_LMDAR4,SWU_LMDAR4,SW1S_LMDAR4,SW2S_LMDAR4,LW_LMDAR4,LWU_LMDAR4,LWBV_LMDAR4,LWC_LMDAR4,LWB_LMDAR4,LWV_LMDAR4,LWVB_LMDAR4,LWVD_LMDAR4,LWVN_LMDAR4 line=2000 + +bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto +bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR +bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -Wf,-O,extendreorder %INLINE +bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt %INLINE +bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder +bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt Index: /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.path =================================================================== --- /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.path (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/arch/arch-SX9_MERCURE.path (revision 1250) @@ -0,0 +1,10 @@ +NETCDF_LIBDIR=/usr/local/SX8/soft/netcdf/lib +NETCDF_INCDIR=/usr/local/SX8/soft/netcdf/include +IOIPSL_INCDIR=$WORKDIR/IOIPSL_MPP/src +IOIPSL_LIBDIR=$WORKDIR/IOIPSL_MPP/src +ORCH_INCDIR=$LMDGCM/../../lib +ORCH_LIBDIR=$LMDGCM/../../lib +OASIS_INCDIR=$LMDGCM/../../prism/SX/build/lib/psmile.$couple +OASIS_LIBDIR=$LMDGCM/../../prism/SX/lib +INCA_LIBDIR=$LMDGCM/../INCA3/config/lib +INCA_INCDIR=$LMDGCM/../INCA3/config/lib Index: /LMDZ4/branches/LMDZ4-dev/bld.cfg =================================================================== --- /LMDZ4/branches/LMDZ4-dev/bld.cfg (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/bld.cfg (revision 1250) @@ -27,5 +27,4 @@ src::filtrez %SRC_PATH/filtrez src::bibio %SRC_PATH/bibio -src::ext_src %EXT_SRC bld::lib::dyn %DYN @@ -39,5 +38,5 @@ bld::target lib%{DYN}.a lib%{PHYS}.a libgrid.a libfiltrez.a libbibio.a bld::target %EXEC%SUFF_NAME.e -bld::exe_dep %{DYN} %{PHYS} grid filtrez bibio ext_src +bld::exe_dep %{DYN} %{PHYS} grid filtrez bibio @@ -60,4 +59,18 @@ bld::tool::fppkeys %CPP_KEY %FPP_DEF + +#bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto +#bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR +#bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -Wf,-O,extendreorder +#bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -f3 +#bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +#bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder +#bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +#bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +#bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt +#bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt + + +inc arch.opt # Pre-process code before analysing dependencies Index: /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F (revision 1250) @@ -158,6 +158,4 @@ REAL,ALLOCATABLE,SAVE :: pcvgt(:,:), pcvgq(:,:,:) c -c REAL,ALLOCATABLE,SAVE :: pvervel(:,:) -c REAL,ALLOCATABLE,SAVE :: zdufi(:,:),zdvfi(:,:) REAL,ALLOCATABLE,SAVE :: zdtfi(:,:),zdqfi(:,:,:) @@ -175,5 +173,4 @@ REAL,ALLOCATABLE,SAVE :: ztfi_omp(:,:) REAL,ALLOCATABLE,SAVE :: zqfi_omp(:,:,:) -c REAL,ALLOCATABLE,SAVE :: pvervel_omp(:,:) REAL,ALLOCATABLE,SAVE :: zdufi_omp(:,:) REAL,ALLOCATABLE,SAVE :: zdvfi_omp(:,:) @@ -252,5 +249,4 @@ ALLOCATE(pcvgu(klon,llm), pcvgv(klon,llm)) ALLOCATE(pcvgt(klon,llm), pcvgq(klon,llm,2)) -c ALLOCATE(pvervel(klon,llm)) ALLOCATE(zdufi(klon,llm),zdvfi(klon,llm)) ALLOCATE(zdtfi(klon,llm),zdqfi(klon,llm,nqtot)) @@ -278,4 +274,6 @@ c$OMP MASTER +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) do ig0=1,klon i=index_i(ig0) @@ -300,4 +298,6 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l = 1, llmp1 +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) do ig0=1,klon i=index_i(ig0) @@ -314,5 +314,6 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm - +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) do ig0=1,klon i=index_i(ig0) @@ -321,5 +322,4 @@ zplay(ig0,l) = preff * pksurcp ** unskap ztfi(ig0,l) = pteta(i,j,l) * pksurcp -c pcvgt(ig0,l) = pdteta(i,j,l) * pksurcp / pmasse(i,j,l) enddo @@ -335,4 +335,6 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) do ig0=1,klon i=index_i(ig0) @@ -344,20 +346,4 @@ ENDDO -c convergence dynamique pour les traceurs "EAU" -! Earth-specific treatment of first 2 tracers (water) - if (planet_type=="earth") then - DO iq=1,2 -c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) - DO l=1,llm - do ig0=1,klon - i=index_i(ig0) - j=index_j(ig0) -c pcvgq(ig0,l,iq) = pdq(i,j,l,iq) / pmasse(i,j,l) - enddo - ENDDO -c$OMP END DO NOWAIT - ENDDO - endif ! of if (planet_type=="earth") - c Geopotentiel calcule par rapport a la surface locale: @@ -378,19 +364,4 @@ c$OMP END DO NOWAIT -c .... Calcul de la vitesse verticale ( en Pa*m*s ou Kg/s ) .... -c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux -c de masse est calclue dans advtrac_p.F -c -cc$OMP DO SCHEDULE(STATIC,OMP_CHUNK) -c DO l=1,llm -c do ig0=1,klon -c i=index_i(ig0) -c j=index_j(ig0) -c pvervel(ig0,l) = pw(i,j,l)*g* unsaire(i,j) -c enddo -c if (is_north_pole) pvervel(1,l)=pw(1,1,l)*g /apoln -c if (is_south_pole) pvervel(klon,l)=pw(1,jjp1,l)*g/apols -c ENDDO -cc$OMP END DO NOWAIT c @@ -406,4 +377,7 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) +!CDIR SPARSE do ig0=kstart,kend i=index_i(ig0) @@ -412,19 +386,19 @@ zufi(ig0,l)= 0.5 *( pucov(iim,j,l)/cu(iim,j) $ + pucov(1,j,l)/cu(1,j) ) -c pcvgu(ig0,l)= 0.5*( pducov(iim,j,l)/cu(iim,j) -c $ + pducov(1,j,l)/cu(1,j) ) else zufi(ig0,l)= 0.5*( pucov(i-1,j,l)/cu(i-1,j) $ + pucov(i,j,l)/cu(i,j) ) -c pcvgu(ig0,l)= 0.5*( pducov(i-1,j,l)/cu(i-1,j) -c $ + pducov(i,j,l)/cu(i,j) ) endif enddo ENDDO c$OMP END DO NOWAIT + c 46.champ v: c ----------- + c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) DO ig0=kstart,kend i=index_i(ig0) @@ -433,6 +407,4 @@ $ + pvcov(i,j,l)/cv(i,j) ) -c pcvgv(ig0+i,l)= 0.5 * ( pdvcov(i,j-1,l)/cv(i,j-1) -c $ + pdvcov(i,j,l)/cv(i,j) ) ENDDO ENDDO @@ -449,21 +421,15 @@ z1(1) =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,1,l)/cv(1,1) -c z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,1,l)/cv(1,1) DO i=2,iim z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,1,l)/cv(i,1) -c z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,1,l)/cv(i,1) ENDDO DO i=1,iim zcos(i) = COS(rlonv(i))*z1(i) -c zcosbis(i)= COS(rlonv(i))*z1bis(i) zsin(i) = SIN(rlonv(i))*z1(i) -c zsinbis(i)= SIN(rlonv(i))*z1bis(i) ENDDO zufi(1,l) = SSUM(iim,zcos,1)/pi -c pcvgu(1,l) = SSUM(iim,zcosbis,1)/pi zvfi(1,l) = SSUM(iim,zsin,1)/pi -c pcvgv(1,l) = SSUM(iim,zsinbis,1)/pi ENDDO @@ -482,22 +448,15 @@ z1(1) =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,jjm,l)/cv(1,jjm) -c z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,jjm,l)/cv(1,jjm) DO i=2,iim - z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm) -c z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,jjm,l)/cv(i,jjm) + z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm) ENDDO DO i=1,iim zcos(i) = COS(rlonv(i))*z1(i) -c zcosbis(i) = COS(rlonv(i))*z1bis(i) zsin(i) = SIN(rlonv(i))*z1(i) -c zsinbis(i) = SIN(rlonv(i))*z1bis(i) ENDDO zufi(klon,l) = SSUM(iim,zcos,1)/pi -c pcvgu(klon,l) = SSUM(iim,zcosbis,1)/pi zvfi(klon,l) = SSUM(iim,zsin,1)/pi -c pcvgv(klon,l) = SSUM(iim,zsinbis,1)/pi - ENDDO c$OMP END DO NOWAIT @@ -506,6 +465,4 @@ IF (is_sequential) THEN - if (planet_type=="earth") then -#ifdef CPP_EARTH c cIM calcul PV a teta=350, 380, 405K @@ -514,6 +471,4 @@ $ ntetaSTD,rtetaSTD,PVteta) c -#endif - endif ENDIF @@ -525,15 +480,4 @@ c --------------------- -cc$OMP PARALLEL DEFAULT(NONE) -cc$OMP+ PRIVATE(i,l,offset,iq) -cc$OMP+ SHARED(klon_omp_nb,nqtot,klon_omp_begin, -cc$OMP+ debut,lafin,rdayvrai,heure,dtphys,zplev,zplay, -cc$OMP+ zphi,zphis,presnivs,clesphy0,zufi,zvfi,ztfi, -cc$OMP+ zqfi,pvervel,zdufi,zdvfi,zdtfi,zdqfi,zdpsrf) - -c PRIVATE(zplev_omp,zplay_omp,zphi_omp,zphis_omp, -c c$OMP+ presnivs_omp,zufi_omp,zvfi_omp,ztfi_omp, -c c$OMP+ zqfi_omp,pvervel_omp,zdufi_omp,zdvfi_omp, -c c$OMP+ zdtfi_omp,zdqfi_omp,zdpsrf_omp) c$OMP BARRIER @@ -550,5 +494,4 @@ allocate(ztfi_omp(klon,llm)) allocate(zqfi_omp(klon,llm,nqtot)) -c allocate(pvervel_omp(klon,llm)) allocate(zdufi_omp(klon,llm)) allocate(zdvfi_omp(klon,llm)) @@ -617,10 +560,4 @@ enddo -c do l=1,llm -c do i=1,klon -c pvervel_omp(i,l)=pvervel(offset+i,l) -c enddo -c enddo - do l=1,llm do i=1,klon @@ -660,5 +597,4 @@ c$OMP BARRIER -cym call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm) if (planet_type=="earth") then @@ -681,5 +617,4 @@ . ztfi_omp, . zqfi_omp, -c . pvervel_omp, c#ifdef INCA . flxwfi_omp, @@ -695,7 +630,4 @@ #endif endif !of if (planet_type=="earth") - -cym call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm) - c$OMP BARRIER @@ -753,10 +685,4 @@ enddo enddo - -c do l=1,llm -c do i=1,klon -c pvervel(offset+i,l)=pvervel_omp(i,l) -c enddo -c enddo do l=1,llm @@ -791,5 +717,4 @@ -cc$OMP END PARALLEL klon=klon_mpi 500 CONTINUE @@ -797,5 +722,4 @@ c$OMP MASTER -cym call WriteField_phy('zdtfi',zdtfi(:,:),llm) call stop_timer(timer_physic) c$OMP END MASTER @@ -913,5 +837,7 @@ DO l=1,llm -!!cdir NODEP +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) +!cdir NODEP do ig0=kstart,kend i=index_i(ig0) @@ -975,11 +901,12 @@ C - +!cdir NODEP DO iq=1,nqtot iiq=niadv(iq) c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm - -!!cdir NODEP +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) +!cdir NODEP DO ig0=kstart,kend i=index_i(ig0) @@ -1009,5 +936,7 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm -!!cdir NODEP +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) +!cdir NODEP do ig0=kstart,kend i=index_i(ig0) @@ -1052,5 +981,7 @@ c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) DO l=1,llm -!!cdir NODEP +!CDIR ON_ADB(index_i) +!CDIR ON_ADB(index_j) +!cdir NODEP do ig0=kstart,kend i=index_i(ig0) Index: /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F (revision 1250) @@ -13,11 +13,17 @@ real sx((n-1)*incx+1),sy((n-1)*incy+1) c + if (incx.eq.1.and.incy.eq.1) then + do 10 i=1,n + sy(i)=sx(i) +10 continue + else iy=1 ix=1 - do 10 i=1,n + do 11 i=1,n sy(iy)=sx(ix) ix=ix+incx iy=iy+incy -10 continue +11 continue + endif c return @@ -32,9 +38,15 @@ c ssum=0. + if (incx.eq.1) then + do 10 i=1,n + ssum=ssum+sx(i) +10 continue + else ix=1 - do 10 i=1,n + do 11 i=1,n ssum=ssum+sx(ix) ix=ix+incx -10 continue +11 continue + endif c return Index: /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/groupeun_p.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/groupeun_p.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/groupeun_p.F (revision 1250) @@ -1,4 +1,5 @@ SUBROUTINE groupeun_p(jjmax,llmax,jjb,jje,q) USE parallel + USE Write_Field_p IMPLICIT NONE @@ -17,5 +18,5 @@ REAL airecs,qs - INTEGER i,j,l,ig,j1,j2,i0,jd + INTEGER i,j,l,ig,ig2,j1,j2,i0,jd c--------------------------------------------------------------------c @@ -37,4 +38,7 @@ LOGICAL, SAVE :: first = .TRUE. !$OMP THREADPRIVATE(first) + INTEGER,SAVE :: i_index(iim,ngroup) + INTEGER :: offset + REAL :: qsum(iim/ngroup) IF (first) THEN @@ -54,32 +58,67 @@ j_start = MAX(jjb, j1-jd) j_finish = MIN(jje, j2-jd) - DO j=j_start, j_finish - DO i0=1,iim,2**(ngroup-ig+1) - qn=0. - DO i=i0,i0+2**(ngroup-ig+1)-1 - qn=qn+q(i,j,l) - ENDDO - DO i=i0,i0+2**(ngroup-ig+1)-1 - q(i,j,l)=qn*airen_tab(i,j,jd) - ENDDO + DO ig2=1,ngroup-ig+1 + offset=2**(ig2-1) + DO j=j_start, j_finish +!CDIR NODEP +!CDIR ON_ADB(q) + DO i0=1,iim,2**ig2 + q(i0,j,l)=q(i0,j,l)+q(i0+offset,j,l) + ENDDO + ENDDO + ENDDO + + DO j=j_start, j_finish +!CDIR NODEP +!CDIR ON_ADB(q) + DO i=1,iim + q(i,j,l)=q(i-MOD(i-1,2**(ngroup-ig+1)),j,l) + ENDDO + ENDDO + + DO j=j_start, j_finish +!CDIR ON_ADB(airen_tab) +!CDIR ON_ADB(q) + DO i=1,iim + q(i,j,l)=q(i,j,l)*airen_tab(i,j,jd) ENDDO q(iip1,j,l)=q(1,j,l) ENDDO - + !c Concerne le pole sud j_start = MAX(1+jjp1-jje-jd, j1-jd) j_finish = MIN(1+jjp1-jjb-jd, j2-jd) - DO j=j_start, j_finish - DO i0=1,iim,2**(ngroup-ig+1) - qs=0. - DO i=i0,i0+2**(ngroup-ig+1)-1 - qs=qs+q(i,jjp1-j+1-jd,l) - ENDDO - DO i=i0,i0+2**(ngroup-ig+1)-1 - q(i,jjp1-j+1-jd,l)=qs*aires_tab(i,jjp1-j+1,jd) - ENDDO + DO ig2=1,ngroup-ig+1 + offset=2**(ig2-1) + DO j=j_start, j_finish +!CDIR NODEP +!CDIR ON_ADB(q) + DO i0=1,iim,2**ig2 + q(i0,jjp1-j+1-jd,l)= q(i0,jjp1-j+1-jd,l) + & +q(i0+offset,jjp1-j+1-jd,l) + ENDDO + ENDDO + ENDDO + + + DO j=j_start, j_finish +!CDIR NODEP +!CDIR ON_ADB(q) + DO i=1,iim + q(i,jjp1-j+1-jd,l)=q(i-MOD(i-1,2**(ngroup-ig+1)), + & jjp1-j+1-jd,l) + ENDDO + ENDDO + + DO j=j_start, j_finish +!CDIR ON_ADB(aires_tab) +!CDIR ON_ADB(q) + DO i=1,iim + q(i,jjp1-j+1-jd,l)=q(i,jjp1-j+1-jd,l)* + & aires_tab(i,jjp1-j+1,jd) ENDDO q(iip1,jjp1-j+1-jd,l)=q(1,jjp1-j+1-jd,l) ENDDO + j1=j2+1 Index: /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F (revision 1250) @@ -327,5 +327,5 @@ c$OMP MASTER ItCount=ItCount+1 - if (MOD(ItCount,1)==1) then + if (MOD(ItCount,12)==0) then debug=.true. else Index: DZ4/branches/LMDZ4-dev/libf/dyn3dpar/pres2lev.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/pres2lev.F (revision 1249) +++ (revision ) @@ -1,84 +1,0 @@ -! $Id$ -! -c****************************************************** - SUBROUTINE pres2lev(varo,varn,lmo,lmn,po,pn, - % ni,nj,ok_invertp) -c -c interpolation lineaire pour passer -c a une nouvelle discretisation verticale pour -c les variables de GCM -c Francois Forget (01/1995) -c MOdif remy roca 12/97 pour passer de pres2sig -c Modif F.Codron 07/08 po en 3D -c********************************************************** - - IMPLICIT NONE - -c Declarations: -c ============== -c -c ARGUMENTS -c """"""""" - LOGICAL, INTENT(IN) :: ok_invertp - INTEGER, INTENT(IN) :: lmo ! dimensions ancienne couches - INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches - INTEGER lmomx ! dimensions ancienne couches - INTEGER lmnmx ! dimensions nouvelle couches - - parameter(lmomx=10000,lmnmx=10000) - - real, INTENT(IN) :: po(ni,nj,lmo) ! niveau de pression ancienne grille - real, INTENT(IN) :: pn(ni,nj,lmn) ! niveau de pression nouvelle grille - - INTEGER, INTENT(IN) :: ni,nj ! nombre de point horizontale - - REAL, INTENT(IN) :: varo(ni,nj,lmo) ! var dans l'ancienne grille - REAL, INTENT(OUT) :: varn(ni,nj,lmn) ! var dans la nouvelle grille - - real zvaro(lmomx),zpo(lmomx) - -c Autres variables -c """""""""""""""" - INTEGER n, ln ,lo, i, j, Nhoriz - REAL coef - -c run -c ==== - do i=1,ni - do j=1,nj -! Inversion de l'ordre des niveaux verticaux - IF (ok_invertp) THEN - do lo=1,lmo - zpo(lo)=po(i,j,lmo+1-lo) - zvaro(lo)=varo(i,j,lmo+1-lo) - enddo - ELSE - do lo=1,lmo - zpo(lo)=po(i,j,lo) - zvaro(lo)=varo(i,j,lo) - enddo - ENDIF - - do ln=1,lmn - if (pn(i,j,ln).ge.zpo(1))then - varn(i,j,ln) = zvaro(1) - else if (pn(i,j,ln).le.zpo(lmo)) then - varn(i,j,ln) = zvaro(lmo) - else - do lo=1,lmo-1 - if ( (pn(i,j,ln).le.zpo(lo)).and. - & (pn(i,j,ln).gt.zpo(lo+1)) )then - coef=(pn(i,j,ln)-zpo(lo)) - & /(zpo(lo+1)-zpo(lo)) - varn(i,j,ln)=zvaro(lo) - & +coef*(zvaro(lo+1)-zvaro(lo)) -c print*,'pn(',ln,')=',pn(i,j,ln),varn(i,j,ln) - end if - enddo - endif - enddo - - enddo - enddo - return - end Index: /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/pres2lev.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/pres2lev.F90 (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/pres2lev.F90 (revision 1250) @@ -0,0 +1,74 @@ +! $Id: pres2lev.F 1179 2009-06-11 14:18:47Z jghattas $ +! +!****************************************************** +SUBROUTINE pres2lev(varo,varn,lmo,lmn,po,pn,ni,nj,ok_invertp) +! +! interpolation lineaire pour passer +! a une nouvelle discretisation verticale pour +! les variables de GCM +! Francois Forget (01/1995) +! MOdif remy roca 12/97 pour passer de pres2sig +! Modif F.Codron 07/08 po en 3D +!********************************************************** + + IMPLICIT NONE + +! Declarations: +! ============== +! +! ARGUMENTS +! """"""""" + LOGICAL, INTENT(IN) :: ok_invertp + INTEGER, INTENT(IN) :: lmo ! dimensions ancienne couches + INTEGER, INTENT(IN) :: lmn ! dimensions nouvelle couches + + REAL, INTENT(IN) :: po(ni*nj,lmo) ! niveau de pression ancienne grille + REAL, INTENT(IN) :: pn(ni*nj,lmn) ! niveau de pression nouvelle grille + + INTEGER, INTENT(IN) :: ni,nj ! nombre de point horizontal + + REAL, INTENT(IN) :: varo(ni*nj,lmo) ! var dans l'ancienne grille + REAL, INTENT(OUT) :: varn(ni*nj,lmn) ! var dans la nouvelle grille + + REAL :: zvaro(ni*nj,lmo),zpo(ni*nj,lmn) + +! Autres variables +! """""""""""""""" + INTEGER :: ln ,lo, k + REAL :: coef + + +! Inversion de l'ordre des niveaux verticaux + IF (ok_invertp) THEN + DO lo=1,lmo + DO k=1,ni*nj + zpo(k,lo)=po(k,lmo+1-lo) + zvaro(k,lo)=varo(k,lmo+1-lo) + ENDDO + ENDDO + ELSE + DO lo=1,lmo + DO k=1,ni*nj + zpo(k,lo)=po(k,lo) + zvaro(k,lo)=varo(k,lo) + ENDDO + ENDDO + ENDIF + + DO ln=1,lmn + DO lo=1,lmo-1 + DO k=1,ni*nj + IF (pn(k,ln) >= zpo(k,1) ) THEN + varn(k,ln) = varo(k,1) + ELSE IF (pn(k,ln) <= zpo(k,lmo)) THEN + varn(k,ln) = zvaro(k,lmo) + ELSE IF ( pn(k,ln) <= zpo(k,lo) .AND. pn(k,ln) > zpo(k,lo+1) ) THEN + coef = (pn(k,ln)-zpo(k,lo)) / (zpo(k,lo+1)-zpo(k,lo)) + varn(k,ln) = zvaro(k,lo) + coef*(zvaro(k,lo+1)-zvaro(k,lo)) + ENDIF + + ENDDO + ENDDO + ENDDO + +END SUBROUTINE pres2lev Index: /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_wrapper.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_wrapper.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_wrapper.F90 (revision 1250) @@ -19,5 +19,5 @@ INTEGER,INTENT(IN) :: nb_vect REAL,INTENT(IN) :: vect(vsize+inc,nb_vect) - COMPLEX,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect) + COMPLEX*16,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect) STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" @@ -29,5 +29,5 @@ INTEGER,INTENT(IN) :: nb_vect REAL,INTENT(INOUT) :: vect(vsize+inc,nb_vect) - COMPLEX,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect) + COMPLEX*16,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect) STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" Index: /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90 (revision 1250) @@ -199,5 +199,5 @@ REAL :: vect(iim+inc,jj_end-jj_begin+1,nbniv) - COMPLEX :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) + COMPLEX*16 :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) INTEGER :: nb_vect INTEGER :: i,j,l @@ -260,6 +260,6 @@ REAL,INTENT(INOUT) :: vect_inout(iim+1,nlat,nbniv) - REAL :: vect(iim+inc,jj_end-jj_begin+1,nbniv) - COMPLEX :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) + REAL :: vect(iim+inc,jj_end-jj_begin+1,nbniv) + COMPLEX*16 :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) INTEGER :: nb_vect INTEGER :: i,j,l Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90 (revision 1250) @@ -51,7 +51,7 @@ REAL, DIMENSION(klon,klev,naero_tot,nbands) :: cg_ae LOGICAL :: soluble - INTEGER :: i, k, ierr, inu, m, mrfspecies + INTEGER :: i, k,n, ierr, inu, m, mrfspecies INTEGER :: spsol, spinsol, spss - INTEGER :: RH_num + INTEGER :: RH_num(klon,klev) INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA @@ -122,8 +122,9 @@ REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./) REAL, PARAMETER :: RH_MAX=95. - REAL:: DELTA, rh, H - REAL:: tau_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of epaisseur optique aerosol - REAL:: piz_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of Single scattering albedo - REAL:: cg_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of Assymetry parameter + REAL:: DELTA(klon,klev), rh(klon,klev), H + REAL:: tau_ae2b_int ! Intermediate computation of epaisseur optique aerosol + REAL:: piz_ae2b_int ! Intermediate computation of Single scattering albedo + REAL:: cg_ae2b_int ! Intermediate computation of Assymetry parameter + REAL :: Fact_RH(nbre_RH) REAL :: zrho REAL :: fac @@ -146,4 +147,7 @@ REAL:: piz_aeri_2bands(nbands,naero_insoluble) !-- unit + INTEGER :: id + LOGICAL :: used_aer(naero_tot) + REAL :: tmp_var, tmp_var_pi DATA presnivs_19/& @@ -586,25 +590,20 @@ DO k=1, klev - DO i=1, klon -! IF (t_seri(i,k).EQ.0.) THEN -! WRITE(lunout,*) 't_seri(i,k)=0 for i=',i,'k=',k -! CALL abort_gcm('aeropt_2bands','t_seri=0',1) -! END IF -! IF (pplay(i,k).EQ.0.) THEN -! WRITE(lunout,*) 'pplay(i,k)=0 for i=',i,'k=',k -! CALL abort_gcm('aeropt_2bands','pplay=0',1) -! END IF - zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 - mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 + DO i=1, klon + zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 +!CDIR UNROLL=naero_spc + mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 !RAF zrho - mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 - ENDDO +!CDIR UNROLL=naero_spc + mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 + zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] + ENDDO ENDDO IF (flag_aerosol .EQ. 1) THEN - nb_aer = 1 + nb_aer = 2 ALLOCATE (aerosol_name(nb_aer)) aerosol_name(1) = id_ASSO4M - + aerosol_name(2) = id_CSSO4M ELSEIF (flag_aerosol .EQ. 2) THEN nb_aer = 2 @@ -650,23 +649,30 @@ ! compute optical_thickness_at_gridpoint_per_species - tau_ae(:,:,:,:)=0. -!RAF - tau_ae_pi(:,:,:,:)=0. - piz_ae(:,:,:,:)=0. - cg_ae(:,:,:,:)=0. - tau_allaer(:,:,:,:)=0. - piz_allaer(:,:,:,:)=0. - cg_allaer(:,:,:,:)=0. - - ! - ! Calculations that need to be done since we are not in the subroutines INCA - ! - ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] - zdp1(:,:)=pdel(:,:)/(gravit*delt) - - + + +!!CDIR ON_ADB(RH_tab) +!CDIR ON_ADB(fact_RH) +!CDIR SHORTLOOP + DO n=1,nbre_RH-1 + fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n)) + ENDDO + + DO k=1, KLEV +!!CDIR ON_ADB(RH_tab) +!CDIR ON_ADB(fact_RH) + DO i=1, KLON + rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) + RH_num(i,k) = INT( rh(i,k)/10. + 1.) + IF (rh(i,k).GT.85.) RH_num(i,k)=10 + IF (rh(i,k).GT.90.) RH_num(i,k)=11 + + DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k)) + ENDDO + ENDDO + + used_aer(:)=.FALSE. + DO m=1,nb_aer ! tau is only computed for each mass - - fac=1.0 + fac=1.0 IF (aerosol_name(m).EQ.id_ASBCM) THEN soluble=.TRUE. @@ -715,130 +721,243 @@ ENDIF - - tau_ae2b_int(:,:,:)=0. - piz_ae2b_int(:,:,:)=0. - cg_ae2b_int(:,:,:)=0. - - DO inu=1,nbands + id=aerosol_name(m) + used_aer(id)=.TRUE. + + + IF (soluble) THEN + + IF (spss.NE.0) THEN + + IF (spss.EQ.1) THEN !accumulation mode + DO k=1, KLEV +!CDIR ON_ADB(A1_ASSSM_b1) +!CDIR ON_ADB(A2_ASSSM_b1) +!CDIR ON_ADB(A3_ASSSM_b1) +!CDIR ON_ADB(B1_ASSSM_b1) +!CDIR ON_ADB(B2_ASSSM_b1) +!CDIR ON_ADB(C1_ASSSM_b1) +!CDIR ON_ADB(C2_ASSSM_b2) +!CDIR ON_ADB(A1_ASSSM_b2) +!CDIR ON_ADB(A2_ASSSM_b2) +!CDIR ON_ADB(A3_ASSSM_b2) +!CDIR ON_ADB(B1_ASSSM_b2) +!CDIR ON_ADB(B2_ASSSM_b2) +!CDIR ON_ADB(C1_ASSSM_b2) +!CDIR ON_ADB(C2_ASSSM_b2) + DO i=1, KLON + H=rh(i,k)/100 + tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + + ! band 1 + tau_ae2b_int=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05) + piz_ae2b_int=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H + cg_ae2b_int=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H + + tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,1) = piz_ae2b_int + cg_ae(i,k,id,1)= cg_ae2b_int + + !band 2 + tau_ae2b_int=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05) + piz_ae2b_int=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H + cg_ae2b_int=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H + + tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,2) = piz_ae2b_int + cg_ae(i,k,id,2)= cg_ae2b_int + + ENDDO + ENDDO + ENDIF + + IF (spss.EQ.2) THEN !coarse mode + DO k=1, KLEV +!CDIR ON_ADB(A1_CSSSM_b1) +!CDIR ON_ADB(A2_CSSSM_b1) +!CDIR ON_ADB(A3_CSSSM_b1) +!CDIR ON_ADB(B1_CSSSM_b1) +!CDIR ON_ADB(B2_CSSSM_b1) +!CDIR ON_ADB(C1_CSSSM_b1) +!CDIR ON_ADB(C2_CSSSM_b2) +!CDIR ON_ADB(A1_CSSSM_b2) +!CDIR ON_ADB(A2_CSSSM_b2) +!CDIR ON_ADB(A3_CSSSM_b2) +!CDIR ON_ADB(B1_CSSSM_b2) +!CDIR ON_ADB(B2_CSSSM_b2) +!CDIR ON_ADB(C1_CSSSM_b2) +!CDIR ON_ADB(C2_CSSSM_b2) + DO i=1, KLON + H=rh(i,k)/100 + tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + ! band 1 + tau_ae2b_int=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05) + piz_ae2b_int=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H + cg_ae2b_int=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H + + tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,1) = piz_ae2b_int + cg_ae(i,k,id,1)= cg_ae2b_int + + ! band 2 + tau_ae2b_int=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05) + piz_ae2b_int=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H + cg_ae2b_int=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H + + tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,2) = piz_ae2b_int + cg_ae(i,k,id,2)= cg_ae2b_int + + ENDDO + ENDDO + ENDIF + + IF (spss.EQ.3) THEN !super coarse mode + DO k=1, KLEV +!CDIR ON_ADB(A1_SSSSM_b1) +!CDIR ON_ADB(A2_SSSSM_b1) +!CDIR ON_ADB(A3_SSSSM_b1) +!CDIR ON_ADB(B1_SSSSM_b1) +!CDIR ON_ADB(B2_SSSSM_b1) +!CDIR ON_ADB(C1_SSSSM_b1) +!CDIR ON_ADB(C2_SSSSM_b2) +!CDIR ON_ADB(A1_SSSSM_b2) +!CDIR ON_ADB(A2_SSSSM_b2) +!CDIR ON_ADB(A3_SSSSM_b2) +!CDIR ON_ADB(B1_SSSSM_b2) +!CDIR ON_ADB(B2_SSSSM_b2) +!CDIR ON_ADB(C1_SSSSM_b2) +!CDIR ON_ADB(C2_SSSSM_b2) + DO i=1, KLON + H=rh(i,k)/100 + tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + + ! band 1 + tau_ae2b_int=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05) + piz_ae2b_int=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H + cg_ae2b_int=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H + + tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,1) = piz_ae2b_int + cg_ae(i,k,id,1)= cg_ae2b_int + + ! band 2 + tau_ae2b_int=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05) + piz_ae2b_int=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H + cg_ae2b_int=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H + + tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,2) = piz_ae2b_int + cg_ae(i,k,id,2)= cg_ae2b_int + + ENDDO + ENDDO + ENDIF + + ELSE + +!CDIR ON_ADB(alpha_aers_2bands) +!CDIR ON_ADB(piz_aers_2bands) +!CDIR ON_ADB(cg_aers_2bands) + DO k=1, KLEV + DO i=1, KLON + tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac + tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac +!CDIR UNROLL=nbands + DO inu=1,nbands + + tau_ae2b_int= alpha_aers_2bands(RH_num(i,k),inu,spsol)+ & + DELTA(i,k)* (alpha_aers_2bands(RH_num(i,k)+1,inu,spsol) - & + alpha_aers_2bands(RH_num(i,k),inu,spsol)) + + piz_ae2b_int = piz_aers_2bands(RH_num(i,k),inu,spsol) + & + DELTA(i,k)* (piz_aers_2bands(RH_num(i,k)+1,inu,spsol) - & + piz_aers_2bands(RH_num(i,k),inu,spsol)) + + cg_ae2b_int = cg_aers_2bands(RH_num(i,k),inu,spsol) + & + DELTA(i,k)* (cg_aers_2bands(RH_num(i,k)+1,inu,spsol) - & + cg_aers_2bands(RH_num(i,k),inu,spsol)) + + tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int + piz_ae(i,k,id,inu) = piz_ae2b_int + cg_ae(i,k,id,inu)= cg_ae2b_int + + ENDDO + ENDDO + ENDDO + + ENDIF + + ELSE ! For all aerosol insoluble components + +!CDIR ON_ADB(alpha_aers_2bands) +!CDIR ON_ADB(piz_aers_2bands) +!CDIR ON_ADB(cg_aers_2bands) DO k=1, KLEV - DO i=1, KLON - - rh=MIN(RHcl(i,k)*100.,RH_MAX) - RH_num = INT( rh/10. + 1.) - - IF (rh.GT.85.) RH_num=10 - IF (rh.GT.90.) RH_num=11 - DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num)) - -! DO inu=1,nbands - IF (soluble) THEN - - ! First optical parameters are computed for seasalt - IF (spss.NE.0) THEN - H=rh/100 - IF (spss.EQ.1) THEN !accumulation mode - ! band 1 - tau_ae2b_int(i,k,1)=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05) - piz_ae2b_int(i,k,1)=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H - cg_ae2b_int(i,k,1)=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H - !band 2 - tau_ae2b_int(i,k,2)=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05) - piz_ae2b_int(i,k,2)=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H - cg_ae2b_int(i,k,2)=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H - ENDIF - IF (spss.EQ.2) THEN !coarse mode - ! band 1 - tau_ae2b_int(i,k,1)=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05) - piz_ae2b_int(i,k,1)=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H - cg_ae2b_int(i,k,1)=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H - ! band 2 - tau_ae2b_int(i,k,2)=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05) - piz_ae2b_int(i,k,2)=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H - cg_ae2b_int(i,k,2)=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H - ENDIF - IF (spss.EQ.3) THEN !super coarse mode - ! band 1 - tau_ae2b_int(i,k,1)=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05) - piz_ae2b_int(i,k,1)=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H - cg_ae2b_int(i,k,1)=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H - ! band 2 - tau_ae2b_int(i,k,2)=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05) - piz_ae2b_int(i,k,2)=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H - cg_ae2b_int(i,k,2)=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H - ENDIF - ELSE - tau_ae2b_int(i,k,inu)= & - alpha_aers_2bands(RH_num,inu,spsol)+ & - DELTA* (alpha_aers_2bands(RH_num+1,inu,spsol) - & - alpha_aers_2bands(RH_num,inu,spsol)) - - piz_ae2b_int(i,k,inu)= & - piz_aers_2bands(RH_num,inu,spsol) + & - DELTA* (piz_aers_2bands(RH_num+1,inu,spsol) - & - piz_aers_2bands(RH_num,inu,spsol)) - - cg_ae2b_int(i,k,inu)= & - cg_aers_2bands(RH_num,inu,spsol) + & - DELTA* (cg_aers_2bands(RH_num+1,inu,spsol) - & - cg_aers_2bands(RH_num,inu,spsol)) - ENDIF - - tau_ae(i,k,aerosol_name(m),inu) = & - mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt* & - tau_ae2b_int(i,k,inu)*fac -!RAF fac - tau_ae_pi(i,k,aerosol_name(m),inu) = & - mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt* & - tau_ae2b_int(i,k,inu)*fac - - ELSE ! For all aerosol insoluble components - tau_ae2b_int(i,k,inu) = alpha_aeri_2bands(inu,spinsol) - piz_ae2b_int(i,k,inu) = piz_aeri_2bands(inu,spinsol) - cg_ae2b_int(i,k,inu) = cg_aeri_2bands(inu,spinsol) - - tau_ae(i,k,aerosol_name(m),inu) = & - mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* & - delt*tau_ae2b_int(i,k,inu)*fac -!RAF fac - tau_ae_pi(i,k,aerosol_name(m),inu) = & - mass_temp_pi(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* & - delt*tau_ae2b_int(i,k,inu)*fac - ENDIF - - piz_ae(i,k,aerosol_name(m),inu) = piz_ae2b_int(i,k,inu) - - cg_ae(i,k,aerosol_name(m),inu)= cg_ae2b_int(i,k,inu) - - ENDDO ! nbands : boucle sur les bandes spectrale - ENDDO ! klon : Boucle sur les points geographiques (grille horizontale) - ENDDO ! klev : Boucle sur les niveaux verticaux - ENDDO ! nb_aer : Boucle sur les masses de traceurs - + DO i=1, KLON + tmp_var=mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)*delt*fac + tmp_var_pi=mass_temp_pi(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)*delt*fac +!CDIR UNROLL=nbands + DO inu=1,nbands + tau_ae2b_int = alpha_aeri_2bands(inu,spinsol) + piz_ae2b_int = piz_aeri_2bands(inu,spinsol) + cg_ae2b_int = cg_aeri_2bands(inu,spinsol) + + tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int + tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int + piz_ae(i,k,id,inu) = piz_ae2b_int + cg_ae(i,k,id,inu)= cg_ae2b_int + ENDDO + ENDDO + ENDDO + + ENDIF ! soluble + + ENDDO ! nb_aer + + DO m=1,nb_aer + IF (.NOT. used_aer(m)) THEN + tau_ae(:,:,:,:)=0. + tau_ae_pi(:,:,:,:)=0. + piz_ae(:,:,:,:)=0. + cg_ae(:,:,:,:)=0. + ENDIF + ENDDO DO inu=1, nbands - DO mrfspecies=1,naero_grp + DO mrfspecies=1,naero_grp + IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER DO k=1, KLEV - DO i=1, KLON - IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ & - tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ & - tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ & - tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)+ & - tau_ae(i,k,id_CIDUSTM,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ & + tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ & + tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ & + tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+ & + tau_ae(i,k,id_SSSSM,inu)+ tau_ae(i,k,id_CIDUSTM,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & - tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ & - tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ & - tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ & - tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ & - tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ & - tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ & - tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ & - tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ & - tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu))/tau_allaer(i,k,mrfspecies,inu) - piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) - - cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ & + piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & + tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ & + tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ & + tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ & + tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ & + tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ & + tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ & + tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ & + tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ & + tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) & + /tau_allaer(i,k,mrfspecies,inu) + piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) + + cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu)+ & tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ & tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ & @@ -851,6 +970,11 @@ tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ & (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - - ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT + + DO k=1, KLEV + DO i=1, KLON !RAF tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M,inu)+ & @@ -892,66 +1016,106 @@ cg_ae(i,k,id_CIDUSTM,inu))/ & (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - -!END RAF - ELSEIF (mrfspecies .EQ. 4) THEN ! = BC - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) & + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 4) THEN ! = BC + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) & +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ & tau_allaer(i,k,mrfspecies,inu) - piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) - cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) *cg_ae(i,k,id_ASBCM,inu)& + piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) + cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) *cg_ae(i,k,id_ASBCM,inu)& +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu))/ & (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - - ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4 - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) & + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4 + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) & +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ & tau_allaer(i,k,mrfspecies,inu) - piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) - cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) *cg_ae(i,k,id_CSSO4M,inu)& + piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) + cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) *cg_ae(i,k,id_CSSO4M,inu)& +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu))/ & (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - - ELSEIF (mrfspecies .EQ. 6) THEN ! = POM - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) & + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 6) THEN ! = POM + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) & +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ & tau_allaer(i,k,mrfspecies,inu) - piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) - cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) *cg_ae(i,k,id_ASPOMM,inu)& + piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) + cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) *cg_ae(i,k,id_ASPOMM,inu)& +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu))/ & (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - - ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu) - cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu) - - ELSEIF (mrfspecies .EQ. 8) THEN ! = SS - tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu) - tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) - piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) & - +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) & - +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ & - tau_allaer(i,k,mrfspecies,inu) - piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) - cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) *cg_ae(i,k,id_ASSSM,inu)& - +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) & - +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ & - (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) - - ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3 - tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary - piz_allaer(i,k,mrfspecies,inu)=0. - cg_allaer(i,k,mrfspecies,inu)=0. - ENDIF - ENDDO + ENDDO ENDDO - ENDDO + + ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu) + cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu) + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 8) THEN ! = SS + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu) + tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) + piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) & + +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) & + +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ & + tau_allaer(i,k,mrfspecies,inu) + piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) + cg_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) *cg_ae(i,k,id_ASSSM,inu)& + +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) & + +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ & + (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) + ENDDO + ENDDO + + ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3 + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary + piz_allaer(i,k,mrfspecies,inu)=0. + cg_allaer(i,k,mrfspecies,inu)=0. + ENDDO + ENDDO + + ELSE + + DO k=1, KLEV + DO i=1, KLON + tau_allaer(i,k,mrfspecies,inu)=0. + piz_allaer(i,k,mrfspecies,inu)=0. + cg_allaer(i,k,mrfspecies,inu)=0. + ENDDO + ENDDO + + ENDIF + + ENDDO ENDDO Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90 (revision 1250) @@ -126,8 +126,8 @@ REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./) - REAL :: DELTA, rh, H - REAL :: tau_ae5wv_int(KLON,KLEV,las) ! Intermediate computation of epaisseur optique aerosol - REAL :: piz_ae5wv_int(KLON,KLEV,las) ! Intermediate single scattering albedo aerosol - REAL :: cg_ae5wv_int(KLON,KLEV,las) ! Intermediate asymmetry parameter aerosol + REAL :: DELTA(klon,klev), rh(klon,klev), H + REAL :: tau_ae5wv_int ! Intermediate computation of epaisseur optique aerosol + REAL :: piz_ae5wv_int ! Intermediate single scattering albedo aerosol + REAL :: cg_ae5wv_int ! Intermediate asymmetry parameter aerosol REAL, PARAMETER :: RH_MAX=95. REAL :: taue670(KLON) ! epaisseur optique aerosol absorption 550 nm @@ -158,6 +158,9 @@ ! Proprietes optiques ! - REAL :: radry = 287.054 ! dry air mass constant - + REAL :: radry = 287.054 + REAL :: tau_tmp ! dry air mass constant + REAL :: fact_RH(nbre_RH) + LOGICAL :: used_tau(naero_spc) + DATA presnivs_19/& 100426.5, 98327.6, 95346.5, 90966.8, 84776.9, & @@ -477,9 +480,5 @@ ! Initialisations ai(:) = 0. - tau_ae5wv_int(:,:,:) = 0. - piz_ae5wv_int(:,:,:) = 0. - cg_ae5wv_int(:,:,:) = 0. tausum(:,:,:) = 0. - tau(:,:,:,:) = 0. @@ -489,5 +488,8 @@ ! IF (pplay(i,k).EQ.0) stop 'stop aeropt_5wv p ' zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 +!CDIR UNROLL=naero_spc mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 + zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] + ENDDO ENDDO @@ -545,7 +547,27 @@ ! Calculations that need to be done since we are not in the subroutines INCA ! - ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] - zdp1=pdel(:,:)/(gravit*delt) - + +!!CDIR ON_ADB(RH_tab) +!CDIR ON_ADB(fact_RH) +!CDIR NOVECTOR + DO RH_Num=1,nbre_RH-1 + fact_RH(RH_num)=1./(RH_tab(RH_num+1)-RH_tab(RH_num)) + ENDDO + + DO k=1, KLEV +!!CDIR ON_ADB(RH_tab) +!CDIR ON_ADB(fact_RH) + DO i=1, KLON + rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) + RH_num = INT( rh(i,k)/10. + 1.) + IF (rh(i,k).GT.85.) RH_num=10 + IF (rh(i,k).GT.90.) RH_num=11 + DELTA(i,k)=(rh(i,k)-RH_tab(RH_num))*fact_RH(RH_num) + ENDDO + ENDDO + +!CDIR SHORTLOOP + used_tau(:)=.FALSE. + DO m=1,nb_aer ! tau is only computed for each mass fac=1.0 @@ -595,96 +617,89 @@ CYCLE ENDIF - + + used_tau(spsol)=.TRUE. DO la=1,las - tau3d(:,:)=0. - piz3d(:,:)=0. - cg3d(:,:)=0. - abs3d(:,:)=0. - - DO k=1, KLEV - DO i=1, KLON - - rh=MIN(RHcl(i,k)*100.,RH_MAX) - RH_num = INT( rh/10. + 1.) - - IF (rh.GT.85.) RH_num=10 - IF (rh.GT.90.) RH_num=11 - DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num)) - - IF (soluble) THEN - - IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm - H=rh/100 - IF (spss.EQ.1) THEN !accumulation mode - tau_ae5wv_int(i,k,la)=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05) - piz_ae5wv_int(i,k,la)=1-B1_ASSSM(k)-B2_ASSSM(k)*H - cg_ae5wv_int(i,k,la)=C1_ASSSM(k)+C2_ASSSM(k)*H - ENDIF - IF (spss.EQ.2) THEN !coarse mode - tau_ae5wv_int(i,k,la)=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05) - piz_ae5wv_int(i,k,la)=1-B1_CSSSM(k)-B2_CSSSM(k)*H - cg_ae5wv_int(i,k,la)=C1_CSSSM(k)+C2_CSSSM(k)*H - ENDIF - IF (spss.EQ.3) THEN !super coarse mode - tau_ae5wv_int(i,k,la)=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05) - piz_ae5wv_int(i,k,la)=1-B1_SSSSM(k)-B2_SSSSM(k)*H - cg_ae5wv_int(i,k,la)=C1_SSSSM(k)+C2_SSSSM(k)*H - ENDIF - ELSE - - tau_ae5wv_int(i,k,la) = & - alpha_aers_5wv(RH_num,la,spsol)+DELTA* & - (alpha_aers_5wv(RH_num+1,la,spsol) - & - alpha_aers_5wv(RH_num,la,spsol)) - - piz_ae5wv_int(i,k,la) = & - piz_aers_5wv(RH_num,la,spsol)+DELTA* & - (piz_aers_5wv(RH_num+1,la,spsol) - & - piz_aers_5wv(RH_num,la,spsol)) - - cg_ae5wv_int(i,k,la) = & - cg_aers_5wv(RH_num,la,spsol)+DELTA* & - (cg_aers_5wv(RH_num+1,la,spsol) - & - cg_aers_5wv(RH_num,la,spsol)) - ENDIF - - tau3d(i,k) = & - mass_temp(i,k,spsol)*1000.*zdp1(i,k)*tau_ae5wv_int(i,k,la)*delt*fac - - ELSE ! For insoluble aerosol - tau_ae5wv_int(i,k,la) = alpha_aeri_5wv(la,spinsol) - piz_ae5wv_int(i,k,la) = piz_aeri_5wv(la,spinsol) - cg_ae5wv_int(i,k,la) = cg_aeri_5wv(la,spinsol) - - tau3d(i,k) = & - mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* & - tau_ae5wv_int(i,k,la)*delt*fac + + IF (soluble) THEN + + IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm + IF (spss.EQ.1) THEN !accumulation mode + DO k=1, KLEV +!CDIR ON_ADB(A1_ASSSM) +!CDIR ON_ADB(A2_ASSSM) +!CDIR ON_ADB(A3_ASSSM) + DO i=1, KLON + H=rh(i,k)/100 + tau_ae5wv_int=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05) + tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & + *tau_ae5wv_int*delt*fac + tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) + ENDDO + ENDDO ENDIF - - ENDDO ! Boucle sur les points géographiques (grille horizontale) - ENDDO ! Boucle sur les niveaux verticaux - - IF (soluble) THEN - - tau(:,:,la,spsol)=tau3d(:,:) - + + IF (spss.EQ.2) THEN !coarse mode + DO k=1, KLEV +!CDIR ON_ADB(A1_CSSSM) +!CDIR ON_ADB(A2_CSSSM) +!CDIR ON_ADB(A3_CSSSM) + DO i=1, KLON + H=rh(i,k)/100 + tau_ae5wv_int=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05) + tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & + *tau_ae5wv_int*delt*fac + tausum(i,la,spsol) = tausum(i,la,spsol)+tau(i,k,la,spsol) + ENDDO + ENDDO + ENDIF + + IF (spss.EQ.3) THEN !super coarse mode + DO k=1, KLEV +!CDIR ON_ADB(A1_SSSSM) +!CDIR ON_ADB(A2_SSSSM) +!CDIR ON_ADB(A3_SSSSM) + DO i=1, KLON + H=rh(i,k)/100 + tau_ae5wv_int=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05) + tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & + *tau_ae5wv_int*delt*fac + tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) + ENDDO + ENDDO + ENDIF + + ELSE DO k=1, KLEV - DO i=1,KLON - tausum(i,la,spsol)=tausum(i,la,spsol)+tau3d(i,k) +!CDIR ON_ADB(alpha_aers_5wv) + DO i=1, KLON + tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* & + (alpha_aers_5wv(RH_num+1,la,spsol) - & + alpha_aers_5wv(RH_num,la,spsol)) + + tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & + *tau_ae5wv_int*delt*fac + tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) ENDDO ENDDO - ELSE ! For insoluble aerosol - tau(:,:,la,naero_soluble+spinsol)=tau3d(:,:) - - DO k=1, KLEV - DO i=1,KLON - tausum(i,la,naero_soluble+spinsol)= & - tausum(i,la,naero_soluble+spinsol)+tau3d(i,k) - ENDDO ! Boucle sur les points géographiques (grille horizontale) - ENDDO ! Boucle sur les niveaux verticaux + ENDIF + + ELSE ! For insoluble aerosol + DO k=1, KLEV +!CDIR ON_ADB(alpha_aeri_5wv) + DO i=1, KLON + tau_ae5wv_int = alpha_aeri_5wv(la,spinsol) + tau(i,k,la,naero_soluble+spinsol) = mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* & + tau_ae5wv_int*delt*fac + tausum(i,la,naero_soluble+spinsol)= tausum(i,la,naero_soluble+spinsol) & + +tau(i,k,la,naero_soluble+spinsol) + ENDDO + ENDDO ENDIF - ENDDO ! boucle sur les longueurs d'onde ENDDO ! Boucle sur les masses de traceurs + + DO m=1,naero_spc + IF (.NOT.used_tau(m)) tau(:,:,:,m)=0. + ENDDO ! ! Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/carbon_cycle_mod.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/carbon_cycle_mod.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/carbon_cycle_mod.F90 (revision 1250) @@ -9,7 +9,7 @@ ! Variables read from parmeter file physiq.def - LOGICAL, SAVE, PUBLIC :: carbon_cycle_tr ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys + LOGICAL, PUBLIC :: carbon_cycle_tr ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys !$OMP THREADPRIVATE(carbon_cycle_tr) - LOGICAL, SAVE,PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) + LOGICAL, PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) !$OMP THREADPRIVATE(carbon_cycle_cpl) LOGICAL :: carbon_cycle_emis_comp=.FALSE. ! Calculation of emission compatible @@ -17,5 +17,7 @@ ! Scalare values when no transport, from physiq.def REAL :: fos_fuel_s ! carbon_cycle_fos_fuel dans physiq.def +!$OMP THREADPRIVATE(fos_fuel_s) REAL :: emis_land_s ! not yet implemented +!$OMP THREADPRIVATE(emis_land_s) INTEGER :: ntr_co2 ! Number of tracers concerning the carbon cycle @@ -25,17 +27,24 @@ INTEGER :: id_fco2_land_use ! - " - INTEGER :: id_fco2_fos_fuel ! - " - - +!$OMP THREADPRIVATE(ntr_co2, id_fco2_tot, id_fco2_ocn, id_fco2_land, id_fco2_land_use, id_fco2_fos_fuel) + REAL, DIMENSION(:), ALLOCATABLE :: fos_fuel ! CO2 fossil fuel emission from file [gC/m2/d] +!$OMP THREADPRIVATE(fos_fuel) REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocn_day ! flux CO2 from ocean for 1 day (cumulated) [gC/m2/d] +!$OMP THREADPRIVATE(fco2_ocn_day) REAL, DIMENSION(:), ALLOCATABLE :: fco2_land_day ! flux CO2 from land for 1 day (cumulated) [gC/m2/d] +!$OMP THREADPRIVATE(fco2_land_day) REAL, DIMENSION(:), ALLOCATABLE :: fco2_lu_day ! Emission from land use change for 1 day (cumulated) [gC/m2/d] +!$OMP THREADPRIVATE(fco2_lu_day) ! Following 2 fields will be initialized in surf_land_orchidee at each time step REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_inst ! flux CO2 from land at one time step +!$OMP THREADPRIVATE(fco2_land_inst) REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_lu_inst ! Emission from land use change at one time step +!$OMP THREADPRIVATE(fco2_lu_inst) ! Calculated co2 field to be send to the ocean via the coupler and to ORCHIDEE REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: co2_send - +!$OMP THREADPRIVATE(co2_send) CONTAINS Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/clouds_gno.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/clouds_gno.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/clouds_gno.F (revision 1250) @@ -47,25 +47,25 @@ INTEGER i,K, n, m - REAL mu(klon), qsat(klon), delta(klon), beta(klon) - real zu2(klon),zv2(klon) - REAL xx(klon), aux(klon), coeff(klon), block(klon) - REAL dist(klon), fprime(klon), det(klon) - REAL pi, u(klon), v(klon), erfcu(klon), erfcv(klon) - REAL xx1(klon), xx2(klon) - real erf,kkk + REAL mu(klon), qsat, delta(klon), beta(klon) + real zu2,zv2 + REAL xx(klon), aux(klon), coeff, block + REAL dist, fprime, det + REAL pi, u, v, erfcu, erfcv + REAL xx1, xx2 + real erf,hsqrtlog_2,v2 real sqrtpi,sqrt2,zx1,zx2,exdel c lconv = true si le calcul a converge (entre autre si qsub < min_q) LOGICAL lconv(klon) -cym - cldf(:,:)=0.0 +!cdir arraycomb + cldf (1:klon,1:ND)=0.0 ! cym + ratqsc(1:klon,1:ND)=0.0 + ptconv(1:klon,1:ND)=.false. +!cdir end arraycomb pi = ACOS(-1.) sqrtpi=sqrt(pi) sqrt2=sqrt(2.) - - ptconv=.false. - ratqsc=0. - + hsqrtlog_2=0.5*SQRT(log(2.)) DO 500 K = 1, ND @@ -74,8 +74,8 @@ mu(i) = R(i,K) mu(i) = MAX(mu(i),min_mu) - qsat(i) = RS(i,K) - qsat(i) = MAX(qsat(i),min_mu) - delta(i) = log(mu(i)/qsat(i)) - enddo ! vector + qsat = RS(i,K) + qsat = MAX(qsat,min_mu) + delta(i) = log(mu(i)/qsat) +c enddo ! vector C @@ -106,5 +106,5 @@ c suffisamment d'eau nuageuse. - do i=1,klon ! vector +c do i=1,klon ! vector IF ( QSUB(i,K) .lt. min_Q ) THEN @@ -124,19 +124,19 @@ c -- roots of equation v > vmax: - det(i) = delta(i) + vmax(i)**2. - if (det(i).LE.0.0) vmax(i) = vmax0 + 1.0 - det(i) = delta(i) + vmax(i)**2. - - if (det(i).LE.0.) then + det = delta(i) + vmax(i)*vmax(i) + if (det.LE.0.0) vmax(i) = vmax0 + 1.0 + det = delta(i) + vmax(i)*vmax(i) + + if (det.LE.0.) then xx(i) = -0.0001 else zx1=-sqrt2*vmax(i) - zx2=SQRT(1.0+delta(i)/(vmax(i)**2.)) - xx1(i)=zx1*(1.0-zx2) - xx2(i)=zx1*(1.0+zx2) - xx(i) = 1.01 * xx1(i) - if ( xx1(i) .GE. 0.0 ) xx(i) = 0.5*xx2(i) + zx2=SQRT(1.0+delta(i)/(vmax(i)*vmax(i))) + xx1=zx1*(1.0-zx2) + xx2=zx1*(1.0+zx2) + xx(i) = 1.01 * xx1 + if ( xx1 .GE. 0.0 ) xx(i) = 0.5*xx2 endif - if (delta(i).LT.0.) xx(i) = -0.5*SQRT(log(2.)) + if (delta(i).LT.0.) xx(i) = -hsqrtlog_2 ENDIF @@ -153,10 +153,11 @@ if (.not.lconv(i)) then - u(i) = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2) - v(i) = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2) - - IF ( v(i) .GT. vmax(i) ) THEN - - IF ( ABS(u(i)) .GT. vmax(i) + u = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2) + v = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2) + v2 = v*v + + IF ( v .GT. vmax(i) ) THEN + + IF ( ABS(u) .GT. vmax(i) : .AND. delta(i) .LT. 0. ) THEN @@ -171,7 +172,7 @@ endif xx(i) = -SQRT(aux(i)) - block(i) = EXP(-v(i)*v(i)) / v(i) / sqrtpi - dist(i) = 0.0 - fprime(i) = 1.0 + block = EXP(-v*v) / v / sqrtpi + dist = 0.0 + fprime = 1.0 ELSE @@ -179,13 +180,13 @@ c -- erfv -> 1.0, use an asymptotic expression of erfv for v large: - erfcu(i) = 1.0-ERF(u(i)) + erfcu = 1.0-ERF(u) c !!! ATTENTION : rajout d'un seuil pour l'exponentiel - aux(i) = sqrtpi*erfcu(i)*EXP(min(v(i)*v(i),100.)) - coeff(i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.) - block(i) = coeff(i) * EXP(-v(i)*v(i)) / v(i) / sqrtpi - dist(i) = v(i) * aux(i) / coeff(i) - beta(i) - fprime(i) = 2.0 / xx(i) * (v(i)**2.) - : * ( coeff(i)*EXP(-delta(i)) - u(i) * aux(i) ) - : / coeff(i) / coeff(i) + aux(i) = sqrtpi*erfcu*EXP(min(v2,100.)) + coeff = 1.0 - 0.5/(v2) + 0.75/(v2*v2) + block = coeff * EXP(-v2) / v / sqrtpi + dist = v * aux(i) / coeff - beta(i) + fprime = 2.0 / xx(i) * (v2) + : * ( EXP(-delta(i)) - u * aux(i) / coeff ) + : / coeff ENDIF ! ABS(u) @@ -195,11 +196,11 @@ c -- general case: - erfcu(i) = 1.0-ERF(u(i)) - erfcv(i) = 1.0-ERF(v(i)) - block(i) = erfcv(i) - dist(i) = erfcu(i) / erfcv(i) - beta(i) - zu2(i)=u(i)*u(i) - zv2(i)=v(i)*v(i) - if(zu2(i).gt.20..or. zv2(i).gt.20.) then + erfcu = 1.0-ERF(u) + erfcv = 1.0-ERF(v) + block = erfcv + dist = erfcu / erfcv - beta(i) + zu2=u*u + zv2=v2 + if(zu2.gt.20..or. zv2.gt.20.) then c print*,'ATTENTION !!! xx(',i,') =', xx(i) c print*,'ATTENTION !!! klon,ND,R,RS,QSUB,PTCONV,RATQSC,CLDF', @@ -207,11 +208,11 @@ c .CLDF(i,k) c print*,'ATTENTION !!! zu2 zv2 =',zu2(i),zv2(i) - zu2(i)=20. - zv2(i)=20. - fprime(i) = 0. + zu2=20. + zv2=20. + fprime = 0. else - fprime(i) = 2. /sqrtpi /xx(i) /erfcv(i)**2. - : * ( erfcv(i)*v(i)*EXP(-zu2(i)) - : - erfcu(i)*u(i)*EXP(-zv2(i)) ) + fprime = 2. /sqrtpi /xx(i) /(erfcv*erfcv) + : * ( erfcv*v*EXP(-zu2) + : - erfcu*u*EXP(-zv2) ) endif ENDIF ! x @@ -223,14 +224,14 @@ ! stop ! endif - if (abs(fprime(i)).lt.1.e-11) then + if (abs(fprime).lt.1.e-11) then ! print*,'avant test fprime<.e-11 ' ! s ,i,k,lconv(i),u(i),v(i),beta(i),fprime(i) ! print*,'klon,ND,R,RS,QSUB', ! s klon,ND,R(i,k),rs(i,k),qsub(i,k) - fprime(i)=sign(1.e-11,fprime(i)) + fprime=sign(1.e-11,fprime) endif - if ( ABS(dist(i)/beta(i)) .LT. epsilon ) then + if ( ABS(dist/beta(i)) .LT. epsilon ) then c print*,'v-u **2',(v(i)-u(i))**2 c print*,'exp v-u **2',exp((v(i)-u(i))**2) @@ -238,9 +239,9 @@ lconv(i)=.true. c borne pour l'exponentielle - ratqsc(i,k)=min(2.*(v(i)-u(i))**2,20.) + ratqsc(i,k)=min(2.*(v-u)*(v-u),20.) ratqsc(i,k)=sqrt(exp(ratqsc(i,k))-1.) - CLDF(i,K) = 0.5 * block(i) + CLDF(i,K) = 0.5 * block else - xx(i) = xx(i) - dist(i)/fprime(i) + xx(i) = xx(i) - dist/fprime endif c print*,'apres test ',i,k,lconv(i) Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/cvltr.F90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/cvltr.F90 (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/cvltr.F90 (revision 1250) @@ -28,5 +28,5 @@ ! Variables locales - REAL,DIMENSION(klon,klev) :: zed +! REAL,DIMENSION(klon,klev) :: zed REAL,DIMENSION(klon,klev,klev) :: zmd REAL,DIMENSION(klon,klev,klev) :: za @@ -34,28 +34,41 @@ REAL,DIMENSION(klon,klev) :: zmfp,zmfu INTEGER :: i,k,j + REAL :: pdtimeRG ! ========================================= ! calcul des tendances liees au downdraft ! ========================================= - zed(:,:)=0. - zmfd(:,:)=0. - zmfa(:,:)=0. - zmfu(:,:)=0. - zmfp(:,:)=0. - zmd(:,:,:)=0. - za(:,:,:)=0. +!cdir collapse + DO j=1,klev + DO i=1,klon +! zed(i,j)=0. + zmfd(i,j)=0. + zmfa(i,j)=0. + zmfu(i,j)=0. + zmfp(i,j)=0. + END DO + END DO +!cdir collapse + DO k=1,klev + DO j=1,klev + DO i=1,klon + zmd(i,j,k)=0. + za (i,j,k)=0. + END DO + END DO + END DO ! entrainement - DO k=1,klev-1 - DO i=1,klon - zed(i,k)=max(0.,mp(i,k)-mp(i,k+1)) - END DO - END DO +! DO k=1,klev-1 +! DO i=1,klon +! zed(i,k)=max(0.,mp(i,k)-mp(i,k+1)) +! END DO +! END DO ! calcul de la matrice d echange ! matrice de distribution de la masse entrainee en k - DO k=1,klev + DO k=1,klev-1 DO i=1,klon - zmd(i,k,k)=zed(i,k) + zmd(i,k,k)=max(0.,mp(i,k)-mp(i,k+1)) END DO END DO @@ -126,7 +139,13 @@ DO k=1, klev DO i=1, klon + dx(i,k)=paprs(i,k)-paprs(i,k+1) + ENDDO + ENDDO + pdtimeRG=pdtime*RG +!cdir collapse + DO k=1, klev + DO i=1, klon dx(i,k)=(zmfd(i,k)+zmfu(i,k) & - +zmfa(i,k)+zmfp(i,k))*pdtime & - *RG/(paprs(i,k)-paprs(i,k+1)) + +zmfa(i,k)+zmfp(i,k))*pdtimeRG/dx(i,k) ! print*,'dx',k,dx(i,k) ENDDO Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/fisrtilp.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/fisrtilp.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/fisrtilp.F (revision 1250) @@ -54,5 +54,5 @@ REAL frac_impa(klon,klev) REAL frac_nucl(klon,klev) - real zct(klon),zcl(klon) + real zct ,zcl cAA c @@ -87,6 +87,6 @@ REAL ztglace, zt(klon) INTEGER nexpo ! exponentiel pour glace/eau - REAL zdz(klon),zrho(klon),ztot(klon), zrhol(klon) - REAL zchau(klon),zfroi(klon),zfice(klon),zneb(klon) + REAL zdz(klon),zrho(klon),ztot , zrhol(klon) + REAL zchau ,zfroi ,zfice(klon),zneb(klon) c LOGICAL appel1er @@ -150,4 +150,5 @@ cAA Initialisation a 1 des coefs des fractions lessivees c +!cdir collapse DO k = 1, klev DO i = 1, klon @@ -161,7 +162,7 @@ c cMAf Initialisation a 0 de zoliq - DO i = 1, klon - zoliq(i)=0. - ENDDO +c DO i = 1, klon +c zoliq(i)=0. +c ENDDO c Determiner les nuages froids par leur temperature c nexpo regle la raideur de la transition eau liquide / eau glace. @@ -173,4 +174,5 @@ c Initialiser les sorties: c +!cdir collapse DO k = 1, klev+1 DO i = 1, klon @@ -180,4 +182,5 @@ ENDDO +!cdir collapse DO k = 1, klev DO i = 1, klon @@ -194,9 +197,10 @@ rain(i) = 0.0 snow(i) = 0.0 - ENDDO + zoliq(i)=0. +c ENDDO c c Initialiser le flux de precipitation a zero c - DO i = 1, klon +c DO i = 1, klon zrfl(i) = 0.0 zneb(i) = seuil_neb @@ -441,26 +445,27 @@ zrhol(i) = zrho(i) * zoliq(i) / zneb(i) - if (ptconv(i,k)) then - zcl(i)=cld_lc_con - zct(i)=1./cld_tau_con - else - zcl(i)=cld_lc_lsc - zct(i)=1./cld_tau_lsc - endif -c quantit�d'eau ��minier. - zchau(i) = zct(i)*dtime/FLOAT(ninter) * zoliq(i) - . *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl(i))**2)) *(1.-zfice(i)) -c meme chose pour la glace. - if (ptconv(i,k)) then - zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) + IF (zneb(i).EQ.seuil_neb) THEN + ztot = 0.0 + ELSE +c quantite d'eau a eliminer: zchau +c meme chose pour la glace: zfroi + if (ptconv(i,k)) then + zcl =cld_lc_con + zct =1./cld_tau_con + zfroi = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) . *fallvc(zrhol(i)) * zfice(i) - else - zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) + else + zcl =cld_lc_lsc + zct =1./cld_tau_lsc + zfroi = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) . *fallvs(zrhol(i)) * zfice(i) - endif - ztot(i) = zchau(i) + zfroi(i) - IF (zneb(i).EQ.seuil_neb) ztot(i) = 0.0 - ztot(i) = MIN(MAX(ztot(i),0.0),zoliq(i)) - zoliq(i) = MAX(zoliq(i)-ztot(i), 0.0) + endif + zchau = zct *dtime/FLOAT(ninter) * zoliq(i) + . *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl )**2)) *(1.-zfice(i)) + ztot = zchau + zfroi + ztot = MAX(ztot ,0.0) + ENDIF + ztot = MIN(ztot,zoliq(i)) + zoliq(i) = MAX(zoliq(i)-ztot , 0.0) radliq(i,k) = radliq(i,k) + zoliq(i)/FLOAT(ninter+1) ENDIF Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/hbtm.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/hbtm.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/hbtm.F (revision 1250) @@ -762,6 +762,6 @@ endif c + qsatbef(i) = qqsat ! bug dans la version orig ??? endif - qsatbef(i) = qqsat camn ???? cette ligne a deja ete faite normalement ? endif Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F (revision 1250) @@ -155,4 +155,6 @@ c$OMP THREADPRIVATE(initialized) cjq-end + INTEGER JK + REAL tmp_ if(.not.initialized) then flag_aer=0. @@ -162,9 +164,10 @@ allocate(ZFSUPAI(KDLON,KFLEV+1)) allocate(ZFSDNAI(KDLON,KFLEV+1)) - ZFSUPAD(:,:)=0. - ZFSDNAD(:,:)=0. - ZFSUPAI(:,:)=0. - ZFSDNAI(:,:)=0. - + DO JK = 1 , KDLON*(KFLEV+1) + ZFSUPAD(JK,1) = 0.0 ! ZFSUPAD(:,:)=0. + ZFSDNAD(JK,1) = 0.0 ! ZFSDNAD(:,:)=0. + ZFSUPAI(JK,1) = 0.0 ! ZFSUPAI(:,:)=0. + ZFSDNAI(JK,1) = 0.0 ! ZFSDNAI(:,:)=0. + END DO endif !rv @@ -179,9 +182,11 @@ IF (MOD(itapsw,swpas).EQ.0) THEN c + tmp_ = 1./( dobson_u * 1e3 * RG) +!cdir collapse DO JK = 1 , KFLEV - DO JL = 1, KDLON - ZCLDSW0(JL,JK) = 0.0 - ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG * PDP(JL,JK) - ENDDO + DO JL = 1, KDLON + ZCLDSW0(JL,JK) = 0.0 + ZOZ(JL,JK) = POZON(JL,JK)*tmp_*PDP(JL,JK) + ENDDO ENDDO C @@ -342,4 +347,6 @@ S PRMU,PSEC,PUD) USE dimphy + USE radiation_AR4_param, only : + S ZPDH2O,ZPDUMG,ZPRH2O,ZPRUMG,RTDH2O,RTDUMG,RTH2O,RTUMG IMPLICIT none cym#include "dimensions.h" @@ -389,22 +396,4 @@ INTEGER jl, jk, jkp1, jkl, jklp1, ja C -C* Prescribed Data: -c - REAL(KIND=8) ZPDH2O,ZPDUMG - SAVE ZPDH2O,ZPDUMG -c$OMP THREADPRIVATE(ZPDH2O,ZPDUMG) - REAL(KIND=8) ZPRH2O,ZPRUMG - SAVE ZPRH2O,ZPRUMG -c$OMP THREADPRIVATE(ZPRH2O,ZPRUMG) - REAL(KIND=8) RTDH2O,RTDUMG - SAVE RTDH2O,RTDUMG -c$OMP THREADPRIVATE(RTDH2O,RTDUMG) - REAL(KIND=8) RTH2O ,RTUMG - SAVE RTH2O ,RTUMG -c$OMP THREADPRIVATE(RTH2O ,RTUMG) - DATA ZPDH2O,ZPDUMG / 0.8 , 0.75 / - DATA ZPRH2O,ZPRUMG / 30000., 30000. / - DATA RTDH2O,RTDUMG / 0.40 , 0.375 / - DATA RTH2O ,RTUMG / 240. , 240. / C ------------------------------------------------------------------ C @@ -537,4 +526,5 @@ S , PFD , PFU) USE dimphy + USE radiation_AR4_param, only : RSUN, RRAY IMPLICIT none cym#include "dimensions.h" @@ -623,21 +613,5 @@ C INTEGER jl, jk, k, jaj, ikm1, ikl -c -c Prescribed Data: -c - REAL(KIND=8) RSUN(2) - SAVE RSUN -c$OMP THREADPRIVATE(RSUN) - REAL(KIND=8) RRAY(2,6) - SAVE RRAY -c$OMP THREADPRIVATE(RRAY) - DATA RSUN(1) / 0.441676 / - DATA RSUN(2) / 0.558324 / - DATA (RRAY(1,K),K=1,6) / - S .428937E-01, .890743E+00,-.288555E+01, - S .522744E+01,-.469173E+01, .161645E+01/ - DATA (RRAY(2,K),K=1,6) / - S .697200E-02, .173297E-01,-.850903E-01, - S .248261E+00,-.302031E+00, .129662E+00/ + C ------------------------------------------------------------------ C @@ -780,4 +754,5 @@ S , PFDOWN,PFUP ) USE dimphy + USE radiation_AR4_param, only : RSUN, RRAY IMPLICIT none cym#include "dimensions.h" @@ -900,20 +875,5 @@ REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11 C -C* Prescribed Data: -C - REAL(KIND=8) RSUN(2) - SAVE RSUN -c$OMP THREADPRIVATE(RSUN) - REAL(KIND=8) RRAY(2,6) - SAVE RRAY -c$OMP THREADPRIVATE(RRAY) - DATA RSUN(1) / 0.441676 / - DATA RSUN(2) / 0.558324 / - DATA (RRAY(1,K),K=1,6) / - S .428937E-01, .890743E+00,-.288555E+01, - S .522744E+01,-.469173E+01, .161645E+01/ - DATA (RRAY(2,K),K=1,6) / - S .697200E-02, .173297E-01,-.850903E-01, - S .248261E+00,-.302031E+00, .129662E+00/ + C C ------------------------------------------------------------------ @@ -1329,4 +1289,5 @@ S , PRK , PRMU0 , PTAUAZ, PTRA1 , PTRA2 ) USE dimphy + USE radiation_AR4_param, only : TAUA, RPIZA, RCGA IMPLICIT none cym#include "dimensions.h" @@ -1398,18 +1359,5 @@ REAL(KIND=8) ZBMU0, ZBMU1, ZRE11 C -C* Prescribed Data for Aerosols: -C - REAL(KIND=8) TAUA(2,5), RPIZA(2,5), RCGA(2,5) - SAVE TAUA, RPIZA, RCGA -c$OMP THREADPRIVATE(TAUA, RPIZA, RCGA) - DATA ((TAUA(IN,JA),JA=1,5),IN=1,2) / - S .730719, .912819, .725059, .745405, .682188 , - S .730719, .912819, .725059, .745405, .682188 / - DATA ((RPIZA(IN,JA),JA=1,5),IN=1,2) / - S .872212, .982545, .623143, .944887, .997975 , - S .872212, .982545, .623143, .944887, .997975 / - DATA ((RCGA (IN,JA),JA=1,5),IN=1,2) / - S .647596, .739002, .580845, .662657, .624246 , - S .647596, .739002, .580845, .662657, .624246 / + C ------------------------------------------------------------------ C @@ -1419,4 +1367,5 @@ 100 CONTINUE C +!cdir collapse DO 103 JK = 1 , KFLEV+1 DO 102 JA = 1 , 6 @@ -2179,4 +2128,5 @@ SUBROUTINE SWTT_LMDAR4 (KNU,KA,PU,PTR) USE dimphy + USE radiation_AR4_param, only : APAD, BPAD, D IMPLICIT none cym#include "dimensions.h" @@ -2225,45 +2175,5 @@ INTEGER jl, i,j C -C* Prescribed Data: -C - REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3) - SAVE APAD, BPAD, D -c$OMP THREADPRIVATE(APAD, BPAD, D) - DATA ((APAD(1,I,J),I=1,3),J=1,7) / - S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, - S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01, - S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00, - S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02, - S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02, - S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02, - S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 / - DATA ((APAD(2,I,J),I=1,3),J=1,7) / - S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, - S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02, - S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00, - S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00, - S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00, - S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00, - S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 / -C - DATA ((BPAD(1,I,J),I=1,3),J=1,7) / - S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, - S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01, - S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00, - S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02, - S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02, - S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02, - S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 / - DATA ((BPAD(2,I,J),I=1,3),J=1,7) / - S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, - S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02, - S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01, - S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00, - S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00, - S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00, - S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 / -c - DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 / - DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 / + C C----------------------------------------------------------------------- @@ -2297,4 +2207,5 @@ SUBROUTINE SWTT1_LMDAR4(KNU,KABS,KIND, PU, PTR) USE dimphy + USE radiation_AR4_param, only : APAD, BPAD, D IMPLICIT none cym#include "dimensions.h" @@ -2345,45 +2256,5 @@ INTEGER jl, ja, i, j, ia C -C* Prescribed Data: -C - REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3) - SAVE APAD, BPAD, D -c$OMP THREADPRIVATE(APAD, BPAD, D) - DATA ((APAD(1,I,J),I=1,3),J=1,7) / - S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, - S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01, - S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00, - S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02, - S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02, - S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02, - S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 / - DATA ((APAD(2,I,J),I=1,3),J=1,7) / - S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, - S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02, - S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00, - S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00, - S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00, - S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00, - S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 / -C - DATA ((BPAD(1,I,J),I=1,3),J=1,7) / - S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, - S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01, - S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00, - S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02, - S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02, - S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02, - S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 / - DATA ((BPAD(2,I,J),I=1,3),J=1,7) / - S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, - S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02, - S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01, - S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00, - S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00, - S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00, - S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 / -c - DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 / - DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 / + C----------------------------------------------------------------------- C @@ -2614,4 +2485,5 @@ S PABCU) USE dimphy + USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT IMPLICIT none cym#include "dimensions.h" @@ -2707,63 +2579,5 @@ REAL(KIND=8) R10E ! DECIMAL/NATURAL LOG.FACTOR PARAMETER (R10E=0.4342945) -c -c Used Data Block: -c - REAL(KIND=8) TREF - SAVE TREF -c$OMP THREADPRIVATE(TREF) - REAL(KIND=8) RT1(2) - SAVE RT1 -c$OMP THREADPRIVATE(RT1) - REAL(KIND=8) RAER(5,5) - SAVE RAER -c$OMP THREADPRIVATE(RAER) - REAL(KIND=8) AT(8,3), BT(8,3) - SAVE AT, BT -c$OMP THREADPRIVATE(AT, BT) - REAL(KIND=8) OCT(4) - SAVE OCT -c$OMP THREADPRIVATE(OCT) - DATA TREF /250.0/ - DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 / - DATA RAER / .038520, .037196, .040532, .054934, .038520 - 1 , .12613 , .18313 , .10357 , .064106, .126130 - 2 , .012579, .013649, .018652, .025181, .012579 - 3 , .011890, .016142, .021105, .028908, .011890 - 4 , .013792, .026810, .052203, .066338, .013792 / - DATA (AT(1,IR),IR=1,3) / - S 0.298199E-02,-.394023E-03,0.319566E-04 / - DATA (BT(1,IR),IR=1,3) / - S-0.106432E-04,0.660324E-06,0.174356E-06 / - DATA (AT(2,IR),IR=1,3) / - S 0.143676E-01,0.366501E-02,-.160822E-02 / - DATA (BT(2,IR),IR=1,3) / - S-0.553979E-04,-.101701E-04,0.920868E-05 / - DATA (AT(3,IR),IR=1,3) / - S 0.197861E-01,0.315541E-02,-.174547E-02 / - DATA (BT(3,IR),IR=1,3) / - S-0.877012E-04,0.513302E-04,0.523138E-06 / - DATA (AT(4,IR),IR=1,3) / - S 0.289560E-01,-.208807E-02,-.121943E-02 / - DATA (BT(4,IR),IR=1,3) / - S-0.165960E-03,0.157704E-03,-.146427E-04 / - DATA (AT(5,IR),IR=1,3) / - S 0.103800E-01,0.436296E-02,-.161431E-02 / - DATA (BT(5,IR),IR=1,3) / - S -.276744E-04,-.327381E-04,0.127646E-04 / - DATA (AT(6,IR),IR=1,3) / - S 0.868859E-02,-.972752E-03,0.000000E-00 / - DATA (BT(6,IR),IR=1,3) / - S -.278412E-04,-.713940E-06,0.117469E-05 / - DATA (AT(7,IR),IR=1,3) / - S 0.250073E-03,0.455875E-03,0.109242E-03 / - DATA (BT(7,IR),IR=1,3) / - S 0.199846E-05,-.216313E-05,0.175991E-06 / - DATA (AT(8,IR),IR=1,3) / - S 0.307423E-01,0.110879E-02,-.322172E-03 / - DATA (BT(8,IR),IR=1,3) / - S-0.108482E-03,0.258096E-05,-.814575E-06 / -c - DATA OCT /-.326E-03, -.102E-05, .137E-02, -.535E-05/ + C----------------------------------------------------------------------- c @@ -3510,4 +3324,5 @@ S , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP) USE dimphy + USE radiation_AR4_param, only : TINTP, XP, GA, GB IMPLICIT none cym#include "dimensions.h" @@ -3604,16 +3419,16 @@ C C* Used Data Block: - REAL(KIND=8) TINTP(11) - SAVE TINTP +c REAL*8 TINTP(11) +c SAVE TINTP c$OMP THREADPRIVATE(TINTP) - REAL(KIND=8) GA(11,16,3), GB(11,16,3) - SAVE GA, GB +c REAL*8 GA(11,16,3), GB(11,16,3) +c SAVE GA, GB c$OMP THREADPRIVATE(GA, GB) - REAL(KIND=8) XP(6,6) - SAVE XP +c REAL*8 XP(6,6) +c SAVE XP c$OMP THREADPRIVATE(XP) c - DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250., - S 262.5, 275., 287.5, 300., 312.5 / +c DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250., +c S 262.5, 275., 287.5, 300., 312.5 / C----------------------------------------------------------------------- C-- WATER VAPOR -- INT.1 -- 0- 500 CM-1 -- FROM ABS225 ---------------- @@ -3628,132 +3443,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 1, 1,IC),IC=1,3) / - S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/ - DATA (GB( 1, 1,IC),IC=1,3) / - S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/ - DATA (GA( 1, 2,IC),IC=1,3) / - S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/ - DATA (GB( 1, 2,IC),IC=1,3) / - S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/ +C DATA (GA( 1, 1,IC),IC=1,3) / +C S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/ +C DATA (GB( 1, 1,IC),IC=1,3) / +C S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/ +C DATA (GA( 1, 2,IC),IC=1,3) / +C S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/ +C DATA (GB( 1, 2,IC),IC=1,3) / +C S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 2, 1,IC),IC=1,3) / - S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/ - DATA (GB( 2, 1,IC),IC=1,3) / - S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/ - DATA (GA( 2, 2,IC),IC=1,3) / - S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/ - DATA (GB( 2, 2,IC),IC=1,3) / - S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/ +C DATA (GA( 2, 1,IC),IC=1,3) / +C S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/ +C DATA (GB( 2, 1,IC),IC=1,3) / +C S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/ +C DATA (GA( 2, 2,IC),IC=1,3) / +C S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/ +C DATA (GB( 2, 2,IC),IC=1,3) / +C S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 3, 1,IC),IC=1,3) / - S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/ - DATA (GB( 3, 1,IC),IC=1,3) / - S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/ - DATA (GA( 3, 2,IC),IC=1,3) / - S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/ - DATA (GB( 3, 2,IC),IC=1,3) / - S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/ +C DATA (GA( 3, 1,IC),IC=1,3) / +C S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/ +C DATA (GB( 3, 1,IC),IC=1,3) / +C S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/ +C DATA (GA( 3, 2,IC),IC=1,3) / +C S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/ +C DATA (GB( 3, 2,IC),IC=1,3) / +C S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 4, 1,IC),IC=1,3) / - S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/ - DATA (GB( 4, 1,IC),IC=1,3) / - S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/ - DATA (GA( 4, 2,IC),IC=1,3) / - S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/ - DATA (GB( 4, 2,IC),IC=1,3) / - S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/ +C DATA (GA( 4, 1,IC),IC=1,3) / +C S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/ +C DATA (GB( 4, 1,IC),IC=1,3) / +C S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/ +C DATA (GA( 4, 2,IC),IC=1,3) / +C S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/ +C DATA (GB( 4, 2,IC),IC=1,3) / +C S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 5, 1,IC),IC=1,3) / - S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/ - DATA (GB( 5, 1,IC),IC=1,3) / - S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/ - DATA (GA( 5, 2,IC),IC=1,3) / - S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/ - DATA (GB( 5, 2,IC),IC=1,3) / - S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/ +C DATA (GA( 5, 1,IC),IC=1,3) / +C S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/ +C DATA (GB( 5, 1,IC),IC=1,3) / +C S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/ +C DATA (GA( 5, 2,IC),IC=1,3) / +C S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/ +C DATA (GB( 5, 2,IC),IC=1,3) / +C S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 6, 1,IC),IC=1,3) / - S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/ - DATA (GB( 6, 1,IC),IC=1,3) / - S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/ - DATA (GA( 6, 2,IC),IC=1,3) / - S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/ - DATA (GB( 6, 2,IC),IC=1,3) / - S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/ +C DATA (GA( 6, 1,IC),IC=1,3) / +C S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/ +C DATA (GB( 6, 1,IC),IC=1,3) / +C S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/ +C DATA (GA( 6, 2,IC),IC=1,3) / +C S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/ +C DATA (GB( 6, 2,IC),IC=1,3) / +C S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 7, 1,IC),IC=1,3) / - S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/ - DATA (GB( 7, 1,IC),IC=1,3) / - S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/ - DATA (GA( 7, 2,IC),IC=1,3) / - S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/ - DATA (GB( 7, 2,IC),IC=1,3) / - S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/ +C DATA (GA( 7, 1,IC),IC=1,3) / +C S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/ +C DATA (GB( 7, 1,IC),IC=1,3) / +C S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/ +C DATA (GA( 7, 2,IC),IC=1,3) / +C S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/ +C DATA (GB( 7, 2,IC),IC=1,3) / +C S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 8, 1,IC),IC=1,3) / - S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/ - DATA (GB( 8, 1,IC),IC=1,3) / - S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/ - DATA (GA( 8, 2,IC),IC=1,3) / - S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/ - DATA (GB( 8, 2,IC),IC=1,3) / - S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/ +C DATA (GA( 8, 1,IC),IC=1,3) / +C S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/ +C DATA (GB( 8, 1,IC),IC=1,3) / +C S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/ +C DATA (GA( 8, 2,IC),IC=1,3) / +C S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/ +C DATA (GB( 8, 2,IC),IC=1,3) / +C S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 9, 1,IC),IC=1,3) / - S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/ - DATA (GB( 9, 1,IC),IC=1,3) / - S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/ - DATA (GA( 9, 2,IC),IC=1,3) / - S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/ - DATA (GB( 9, 2,IC),IC=1,3) / - S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/ +C DATA (GA( 9, 1,IC),IC=1,3) / +C S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/ +C DATA (GB( 9, 1,IC),IC=1,3) / +C S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/ +C DATA (GA( 9, 2,IC),IC=1,3) / +C S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/ +C DATA (GB( 9, 2,IC),IC=1,3) / +C S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA(10, 1,IC),IC=1,3) / - S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/ - DATA (GB(10, 1,IC),IC=1,3) / - S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/ - DATA (GA(10, 2,IC),IC=1,3) / - S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/ - DATA (GB(10, 2,IC),IC=1,3) / - S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/ +C DATA (GA(10, 1,IC),IC=1,3) / +C S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/ +C DATA (GB(10, 1,IC),IC=1,3) / +C S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/ +C DATA (GA(10, 2,IC),IC=1,3) / +C S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/ +C DATA (GB(10, 2,IC),IC=1,3) / +C S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/ C C----- INTERVAL = 1 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA(11, 1,IC),IC=1,3) / - S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/ - DATA (GB(11, 1,IC),IC=1,3) / - S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/ - DATA (GA(11, 2,IC),IC=1,3) / - S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/ - DATA (GB(11, 2,IC),IC=1,3) / - S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/ +C DATA (GA(11, 1,IC),IC=1,3) / +C S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/ +C DATA (GB(11, 1,IC),IC=1,3) / +C S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/ +C DATA (GA(11, 2,IC),IC=1,3) / +C S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/ +C DATA (GB(11, 2,IC),IC=1,3) / +C S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/ C C @@ -3770,132 +3585,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 1, 3,IC),IC=1,3) / - S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/ - DATA (GB( 1, 3,IC),IC=1,3) / - S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/ - DATA (GA( 1, 4,IC),IC=1,3) / - S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/ - DATA (GB( 1, 4,IC),IC=1,3) / - S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/ +C DATA (GA( 1, 3,IC),IC=1,3) / +C S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/ +C DATA (GB( 1, 3,IC),IC=1,3) / +C S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/ +C DATA (GA( 1, 4,IC),IC=1,3) / +C S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/ +C DATA (GB( 1, 4,IC),IC=1,3) / +C S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 2, 3,IC),IC=1,3) / - S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/ - DATA (GB( 2, 3,IC),IC=1,3) / - S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/ - DATA (GA( 2, 4,IC),IC=1,3) / - S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/ - DATA (GB( 2, 4,IC),IC=1,3) / - S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/ +C DATA (GA( 2, 3,IC),IC=1,3) / +C S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/ +C DATA (GB( 2, 3,IC),IC=1,3) / +C S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/ +C DATA (GA( 2, 4,IC),IC=1,3) / +C S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/ +C DATA (GB( 2, 4,IC),IC=1,3) / +C S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 3, 3,IC),IC=1,3) / - S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/ - DATA (GB( 3, 3,IC),IC=1,3) / - S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/ - DATA (GA( 3, 4,IC),IC=1,3) / - S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/ - DATA (GB( 3, 4,IC),IC=1,3) / - S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/ +C DATA (GA( 3, 3,IC),IC=1,3) / +C S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/ +C DATA (GB( 3, 3,IC),IC=1,3) / +C S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/ +C DATA (GA( 3, 4,IC),IC=1,3) / +C S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/ +C DATA (GB( 3, 4,IC),IC=1,3) / +C S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 4, 3,IC),IC=1,3) / - S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/ - DATA (GB( 4, 3,IC),IC=1,3) / - S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/ - DATA (GA( 4, 4,IC),IC=1,3) / - S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/ - DATA (GB( 4, 4,IC),IC=1,3) / - S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/ +C DATA (GA( 4, 3,IC),IC=1,3) / +C S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/ +C DATA (GB( 4, 3,IC),IC=1,3) / +C S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/ +C DATA (GA( 4, 4,IC),IC=1,3) / +C S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/ +C DATA (GB( 4, 4,IC),IC=1,3) / +C S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 5, 3,IC),IC=1,3) / - S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/ - DATA (GB( 5, 3,IC),IC=1,3) / - S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/ - DATA (GA( 5, 4,IC),IC=1,3) / - S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/ - DATA (GB( 5, 4,IC),IC=1,3) / - S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/ +C DATA (GA( 5, 3,IC),IC=1,3) / +C S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/ +C DATA (GB( 5, 3,IC),IC=1,3) / +C S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/ +C DATA (GA( 5, 4,IC),IC=1,3) / +C S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/ +C DATA (GB( 5, 4,IC),IC=1,3) / +C S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 6, 3,IC),IC=1,3) / - S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/ - DATA (GB( 6, 3,IC),IC=1,3) / - S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/ - DATA (GA( 6, 4,IC),IC=1,3) / - S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/ - DATA (GB( 6, 4,IC),IC=1,3) / - S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/ +C DATA (GA( 6, 3,IC),IC=1,3) / +C S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/ +C DATA (GB( 6, 3,IC),IC=1,3) / +C S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/ +C DATA (GA( 6, 4,IC),IC=1,3) / +C S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/ +C DATA (GB( 6, 4,IC),IC=1,3) / +C S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 7, 3,IC),IC=1,3) / - S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/ - DATA (GB( 7, 3,IC),IC=1,3) / - S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/ - DATA (GA( 7, 4,IC),IC=1,3) / - S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/ - DATA (GB( 7, 4,IC),IC=1,3) / - S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/ +C DATA (GA( 7, 3,IC),IC=1,3) / +C S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/ +C DATA (GB( 7, 3,IC),IC=1,3) / +C S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/ +C DATA (GA( 7, 4,IC),IC=1,3) / +C S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/ +C DATA (GB( 7, 4,IC),IC=1,3) / +C S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 8, 3,IC),IC=1,3) / - S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/ - DATA (GB( 8, 3,IC),IC=1,3) / - S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/ - DATA (GA( 8, 4,IC),IC=1,3) / - S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/ - DATA (GB( 8, 4,IC),IC=1,3) / - S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/ +C DATA (GA( 8, 3,IC),IC=1,3) / +C S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/ +C DATA (GB( 8, 3,IC),IC=1,3) / +C S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/ +C DATA (GA( 8, 4,IC),IC=1,3) / +C S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/ +C DATA (GB( 8, 4,IC),IC=1,3) / +C S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 9, 3,IC),IC=1,3) / - S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/ - DATA (GB( 9, 3,IC),IC=1,3) / - S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/ - DATA (GA( 9, 4,IC),IC=1,3) / - S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/ - DATA (GB( 9, 4,IC),IC=1,3) / - S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/ +C DATA (GA( 9, 3,IC),IC=1,3) / +C S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/ +C DATA (GB( 9, 3,IC),IC=1,3) / +C S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/ +C DATA (GA( 9, 4,IC),IC=1,3) / +C S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/ +C DATA (GB( 9, 4,IC),IC=1,3) / +C S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(10, 3,IC),IC=1,3) / - S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/ - DATA (GB(10, 3,IC),IC=1,3) / - S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/ - DATA (GA(10, 4,IC),IC=1,3) / - S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/ - DATA (GB(10, 4,IC),IC=1,3) / - S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/ +C DATA (GA(10, 3,IC),IC=1,3) / +C S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/ +C DATA (GB(10, 3,IC),IC=1,3) / +C S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/ +C DATA (GA(10, 4,IC),IC=1,3) / +C S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/ +C DATA (GB(10, 4,IC),IC=1,3) / +C S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(11, 3,IC),IC=1,3) / - S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/ - DATA (GB(11, 3,IC),IC=1,3) / - S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/ - DATA (GA(11, 4,IC),IC=1,3) / - S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/ - DATA (GB(11, 4,IC),IC=1,3) / - S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/ +C DATA (GA(11, 3,IC),IC=1,3) / +C S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/ +C DATA (GB(11, 3,IC),IC=1,3) / +C S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/ +C DATA (GA(11, 4,IC),IC=1,3) / +C S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/ +C DATA (GB(11, 4,IC),IC=1,3) / +C S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/ C C @@ -3916,132 +3731,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 1, 7,IC),IC=1,3) / - S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/ - DATA (GB( 1, 7,IC),IC=1,3) / - S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/ - DATA (GA( 1, 8,IC),IC=1,3) / - S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/ - DATA (GB( 1, 8,IC),IC=1,3) / - S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/ +C DATA (GA( 1, 7,IC),IC=1,3) / +C S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/ +C DATA (GB( 1, 7,IC),IC=1,3) / +C S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/ +C DATA (GA( 1, 8,IC),IC=1,3) / +C S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/ +C DATA (GB( 1, 8,IC),IC=1,3) / +C S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 2, 7,IC),IC=1,3) / - S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/ - DATA (GB( 2, 7,IC),IC=1,3) / - S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/ - DATA (GA( 2, 8,IC),IC=1,3) / - S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/ - DATA (GB( 2, 8,IC),IC=1,3) / - S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/ +C DATA (GA( 2, 7,IC),IC=1,3) / +C S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/ +C DATA (GB( 2, 7,IC),IC=1,3) / +C S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/ +C DATA (GA( 2, 8,IC),IC=1,3) / +C S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/ +C DATA (GB( 2, 8,IC),IC=1,3) / +C S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 3, 7,IC),IC=1,3) / - S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/ - DATA (GB( 3, 7,IC),IC=1,3) / - S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/ - DATA (GA( 3, 8,IC),IC=1,3) / - S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/ - DATA (GB( 3, 8,IC),IC=1,3) / - S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/ +C DATA (GA( 3, 7,IC),IC=1,3) / +C S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/ +C DATA (GB( 3, 7,IC),IC=1,3) / +C S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/ +C DATA (GA( 3, 8,IC),IC=1,3) / +C S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/ +C DATA (GB( 3, 8,IC),IC=1,3) / +C S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 4, 7,IC),IC=1,3) / - S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/ - DATA (GB( 4, 7,IC),IC=1,3) / - S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/ - DATA (GA( 4, 8,IC),IC=1,3) / - S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/ - DATA (GB( 4, 8,IC),IC=1,3) / - S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/ +C DATA (GA( 4, 7,IC),IC=1,3) / +C S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/ +C DATA (GB( 4, 7,IC),IC=1,3) / +C S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/ +C DATA (GA( 4, 8,IC),IC=1,3) / +C S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/ +C DATA (GB( 4, 8,IC),IC=1,3) / +C S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 5, 7,IC),IC=1,3) / - S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/ - DATA (GB( 5, 7,IC),IC=1,3) / - S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/ - DATA (GA( 5, 8,IC),IC=1,3) / - S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/ - DATA (GB( 5, 8,IC),IC=1,3) / - S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/ +C DATA (GA( 5, 7,IC),IC=1,3) / +C S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/ +C DATA (GB( 5, 7,IC),IC=1,3) / +C S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/ +C DATA (GA( 5, 8,IC),IC=1,3) / +C S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/ +C DATA (GB( 5, 8,IC),IC=1,3) / +C S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 6, 7,IC),IC=1,3) / - S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/ - DATA (GB( 6, 7,IC),IC=1,3) / - S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/ - DATA (GA( 6, 8,IC),IC=1,3) / - S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/ - DATA (GB( 6, 8,IC),IC=1,3) / - S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/ +C DATA (GA( 6, 7,IC),IC=1,3) / +C S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/ +C DATA (GB( 6, 7,IC),IC=1,3) / +C S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/ +C DATA (GA( 6, 8,IC),IC=1,3) / +C S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/ +C DATA (GB( 6, 8,IC),IC=1,3) / +C S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 7, 7,IC),IC=1,3) / - S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/ - DATA (GB( 7, 7,IC),IC=1,3) / - S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/ - DATA (GA( 7, 8,IC),IC=1,3) / - S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/ - DATA (GB( 7, 8,IC),IC=1,3) / - S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/ +C DATA (GA( 7, 7,IC),IC=1,3) / +C S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/ +C DATA (GB( 7, 7,IC),IC=1,3) / +C S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/ +C DATA (GA( 7, 8,IC),IC=1,3) / +C S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/ +C DATA (GB( 7, 8,IC),IC=1,3) / +C S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 8, 7,IC),IC=1,3) / - S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/ - DATA (GB( 8, 7,IC),IC=1,3) / - S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/ - DATA (GA( 8, 8,IC),IC=1,3) / - S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/ - DATA (GB( 8, 8,IC),IC=1,3) / - S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/ +C DATA (GA( 8, 7,IC),IC=1,3) / +C S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/ +C DATA (GB( 8, 7,IC),IC=1,3) / +C S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/ +C DATA (GA( 8, 8,IC),IC=1,3) / +C S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/ +C DATA (GB( 8, 8,IC),IC=1,3) / +C S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 9, 7,IC),IC=1,3) / - S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/ - DATA (GB( 9, 7,IC),IC=1,3) / - S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/ - DATA (GA( 9, 8,IC),IC=1,3) / - S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/ - DATA (GB( 9, 8,IC),IC=1,3) / - S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/ +C DATA (GA( 9, 7,IC),IC=1,3) / +C S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/ +C DATA (GB( 9, 7,IC),IC=1,3) / +C S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/ +C DATA (GA( 9, 8,IC),IC=1,3) / +C S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/ +C DATA (GB( 9, 8,IC),IC=1,3) / +C S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(10, 7,IC),IC=1,3) / - S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/ - DATA (GB(10, 7,IC),IC=1,3) / - S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/ - DATA (GA(10, 8,IC),IC=1,3) / - S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/ - DATA (GB(10, 8,IC),IC=1,3) / - S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/ +C DATA (GA(10, 7,IC),IC=1,3) / +C S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/ +C DATA (GB(10, 7,IC),IC=1,3) / +C S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/ +C DATA (GA(10, 8,IC),IC=1,3) / +C S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/ +C DATA (GB(10, 8,IC),IC=1,3) / +C S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/ C C----- INTERVAL = 3 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(11, 7,IC),IC=1,3) / - S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/ - DATA (GB(11, 7,IC),IC=1,3) / - S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/ - DATA (GA(11, 8,IC),IC=1,3) / - S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/ - DATA (GB(11, 8,IC),IC=1,3) / - S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/ +C DATA (GA(11, 7,IC),IC=1,3) / +C S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/ +C DATA (GB(11, 7,IC),IC=1,3) / +C S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/ +C DATA (GA(11, 8,IC),IC=1,3) / +C S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/ +C DATA (GB(11, 8,IC),IC=1,3) / +C S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/ C C @@ -4053,132 +3868,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 1, 9,IC),IC=1,3) / - S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/ - DATA (GB( 1, 9,IC),IC=1,3) / - S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/ - DATA (GA( 1,10,IC),IC=1,3) / - S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/ - DATA (GB( 1,10,IC),IC=1,3) / - S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/ +C DATA (GA( 1, 9,IC),IC=1,3) / +C S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/ +C DATA (GB( 1, 9,IC),IC=1,3) / +C S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/ +C DATA (GA( 1,10,IC),IC=1,3) / +C S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/ +C DATA (GB( 1,10,IC),IC=1,3) / +C S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 2, 9,IC),IC=1,3) / - S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/ - DATA (GB( 2, 9,IC),IC=1,3) / - S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/ - DATA (GA( 2,10,IC),IC=1,3) / - S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/ - DATA (GB( 2,10,IC),IC=1,3) / - S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/ +C DATA (GA( 2, 9,IC),IC=1,3) / +C S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/ +C DATA (GB( 2, 9,IC),IC=1,3) / +C S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/ +C DATA (GA( 2,10,IC),IC=1,3) / +C S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/ +C DATA (GB( 2,10,IC),IC=1,3) / +C S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 3, 9,IC),IC=1,3) / - S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/ - DATA (GB( 3, 9,IC),IC=1,3) / - S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/ - DATA (GA( 3,10,IC),IC=1,3) / - S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/ - DATA (GB( 3,10,IC),IC=1,3) / - S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/ +C DATA (GA( 3, 9,IC),IC=1,3) / +C S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/ +C DATA (GB( 3, 9,IC),IC=1,3) / +C S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/ +C DATA (GA( 3,10,IC),IC=1,3) / +C S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/ +C DATA (GB( 3,10,IC),IC=1,3) / +C S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 4, 9,IC),IC=1,3) / - S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/ - DATA (GB( 4, 9,IC),IC=1,3) / - S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/ - DATA (GA( 4,10,IC),IC=1,3) / - S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/ - DATA (GB( 4,10,IC),IC=1,3) / - S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/ +C DATA (GA( 4, 9,IC),IC=1,3) / +C S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/ +C DATA (GB( 4, 9,IC),IC=1,3) / +C S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/ +C DATA (GA( 4,10,IC),IC=1,3) / +C S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/ +C DATA (GB( 4,10,IC),IC=1,3) / +C S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 5, 9,IC),IC=1,3) / - S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/ - DATA (GB( 5, 9,IC),IC=1,3) / - S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/ - DATA (GA( 5,10,IC),IC=1,3) / - S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/ - DATA (GB( 5,10,IC),IC=1,3) / - S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/ +C DATA (GA( 5, 9,IC),IC=1,3) / +C S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/ +C DATA (GB( 5, 9,IC),IC=1,3) / +C S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/ +C DATA (GA( 5,10,IC),IC=1,3) / +C S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/ +C DATA (GB( 5,10,IC),IC=1,3) / +C S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 6, 9,IC),IC=1,3) / - S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/ - DATA (GB( 6, 9,IC),IC=1,3) / - S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/ - DATA (GA( 6,10,IC),IC=1,3) / - S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/ - DATA (GB( 6,10,IC),IC=1,3) / - S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/ +C DATA (GA( 6, 9,IC),IC=1,3) / +C S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/ +C DATA (GB( 6, 9,IC),IC=1,3) / +C S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/ +C DATA (GA( 6,10,IC),IC=1,3) / +C S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/ +C DATA (GB( 6,10,IC),IC=1,3) / +C S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 7, 9,IC),IC=1,3) / - S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/ - DATA (GB( 7, 9,IC),IC=1,3) / - S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/ - DATA (GA( 7,10,IC),IC=1,3) / - S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/ - DATA (GB( 7,10,IC),IC=1,3) / - S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/ +C DATA (GA( 7, 9,IC),IC=1,3) / +C S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/ +C DATA (GB( 7, 9,IC),IC=1,3) / +C S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/ +C DATA (GA( 7,10,IC),IC=1,3) / +C S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/ +C DATA (GB( 7,10,IC),IC=1,3) / +C S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 8, 9,IC),IC=1,3) / - S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/ - DATA (GB( 8, 9,IC),IC=1,3) / - S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/ - DATA (GA( 8,10,IC),IC=1,3) / - S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/ - DATA (GB( 8,10,IC),IC=1,3) / - S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/ +C DATA (GA( 8, 9,IC),IC=1,3) / +C S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/ +C DATA (GB( 8, 9,IC),IC=1,3) / +C S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/ +C DATA (GA( 8,10,IC),IC=1,3) / +C S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/ +C DATA (GB( 8,10,IC),IC=1,3) / +C S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA( 9, 9,IC),IC=1,3) / - S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/ - DATA (GB( 9, 9,IC),IC=1,3) / - S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/ - DATA (GA( 9,10,IC),IC=1,3) / - S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/ - DATA (GB( 9,10,IC),IC=1,3) / - S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/ +C DATA (GA( 9, 9,IC),IC=1,3) / +C S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/ +C DATA (GB( 9, 9,IC),IC=1,3) / +C S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/ +C DATA (GA( 9,10,IC),IC=1,3) / +C S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/ +C DATA (GB( 9,10,IC),IC=1,3) / +C S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(10, 9,IC),IC=1,3) / - S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/ - DATA (GB(10, 9,IC),IC=1,3) / - S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/ - DATA (GA(10,10,IC),IC=1,3) / - S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/ - DATA (GB(10,10,IC),IC=1,3) / - S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/ +C DATA (GA(10, 9,IC),IC=1,3) / +C S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/ +C DATA (GB(10, 9,IC),IC=1,3) / +C S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/ +C DATA (GA(10,10,IC),IC=1,3) / +C S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/ +C DATA (GB(10,10,IC),IC=1,3) / +C S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 - DATA (GA(11, 9,IC),IC=1,3) / - S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/ - DATA (GB(11, 9,IC),IC=1,3) / - S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/ - DATA (GA(11,10,IC),IC=1,3) / - S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/ - DATA (GB(11,10,IC),IC=1,3) / - S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/ +C DATA (GA(11, 9,IC),IC=1,3) / +C S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/ +C DATA (GB(11, 9,IC),IC=1,3) / +C S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/ +C DATA (GA(11,10,IC),IC=1,3) / +C S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/ +C DATA (GB(11,10,IC),IC=1,3) / +C S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/ C C @@ -4193,132 +4008,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 1, 5,IC),IC=1,3) / - S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/ - DATA (GB( 1, 5,IC),IC=1,3) / - S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/ - DATA (GA( 1, 6,IC),IC=1,3) / - S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/ - DATA (GB( 1, 6,IC),IC=1,3) / - S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/ +C DATA (GA( 1, 5,IC),IC=1,3) / +C S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/ +C DATA (GB( 1, 5,IC),IC=1,3) / +C S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/ +C DATA (GA( 1, 6,IC),IC=1,3) / +C S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/ +C DATA (GB( 1, 6,IC),IC=1,3) / +C S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 2, 5,IC),IC=1,3) / - S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/ - DATA (GB( 2, 5,IC),IC=1,3) / - S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/ - DATA (GA( 2, 6,IC),IC=1,3) / - S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/ - DATA (GB( 2, 6,IC),IC=1,3) / - S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/ +C DATA (GA( 2, 5,IC),IC=1,3) / +C S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/ +C DATA (GB( 2, 5,IC),IC=1,3) / +C S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/ +C DATA (GA( 2, 6,IC),IC=1,3) / +C S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/ +C DATA (GB( 2, 6,IC),IC=1,3) / +C S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 3, 5,IC),IC=1,3) / - S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/ - DATA (GB( 3, 5,IC),IC=1,3) / - S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/ - DATA (GA( 3, 6,IC),IC=1,3) / - S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/ - DATA (GB( 3, 6,IC),IC=1,3) / - S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/ +C DATA (GA( 3, 5,IC),IC=1,3) / +C S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/ +C DATA (GB( 3, 5,IC),IC=1,3) / +C S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/ +C DATA (GA( 3, 6,IC),IC=1,3) / +C S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/ +C DATA (GB( 3, 6,IC),IC=1,3) / +C S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 4, 5,IC),IC=1,3) / - S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/ - DATA (GB( 4, 5,IC),IC=1,3) / - S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/ - DATA (GA( 4, 6,IC),IC=1,3) / - S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/ - DATA (GB( 4, 6,IC),IC=1,3) / - S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/ +C DATA (GA( 4, 5,IC),IC=1,3) / +C S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/ +C DATA (GB( 4, 5,IC),IC=1,3) / +C S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/ +C DATA (GA( 4, 6,IC),IC=1,3) / +C S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/ +C DATA (GB( 4, 6,IC),IC=1,3) / +C S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 5, 5,IC),IC=1,3) / - S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/ - DATA (GB( 5, 5,IC),IC=1,3) / - S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/ - DATA (GA( 5, 6,IC),IC=1,3) / - S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/ - DATA (GB( 5, 6,IC),IC=1,3) / - S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/ +C DATA (GA( 5, 5,IC),IC=1,3) / +C S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/ +C DATA (GB( 5, 5,IC),IC=1,3) / +C S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/ +C DATA (GA( 5, 6,IC),IC=1,3) / +C S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/ +C DATA (GB( 5, 6,IC),IC=1,3) / +C S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 6, 5,IC),IC=1,3) / - S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/ - DATA (GB( 6, 5,IC),IC=1,3) / - S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/ - DATA (GA( 6, 6,IC),IC=1,3) / - S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/ - DATA (GB( 6, 6,IC),IC=1,3) / - S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/ +C DATA (GA( 6, 5,IC),IC=1,3) / +C S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/ +C DATA (GB( 6, 5,IC),IC=1,3) / +C S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/ +C DATA (GA( 6, 6,IC),IC=1,3) / +C S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/ +C DATA (GB( 6, 6,IC),IC=1,3) / +C S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 7, 5,IC),IC=1,3) / - S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/ - DATA (GB( 7, 5,IC),IC=1,3) / - S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/ - DATA (GA( 7, 6,IC),IC=1,3) / - S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/ - DATA (GB( 7, 6,IC),IC=1,3) / - S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/ +C DATA (GA( 7, 5,IC),IC=1,3) / +C S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/ +C DATA (GB( 7, 5,IC),IC=1,3) / +C S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/ +C DATA (GA( 7, 6,IC),IC=1,3) / +C S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/ +C DATA (GB( 7, 6,IC),IC=1,3) / +C S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 8, 5,IC),IC=1,3) / - S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/ - DATA (GB( 8, 5,IC),IC=1,3) / - S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/ - DATA (GA( 8, 6,IC),IC=1,3) / - S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/ - DATA (GB( 8, 6,IC),IC=1,3) / - S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/ +C DATA (GA( 8, 5,IC),IC=1,3) / +C S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/ +C DATA (GB( 8, 5,IC),IC=1,3) / +C S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/ +C DATA (GA( 8, 6,IC),IC=1,3) / +C S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/ +C DATA (GB( 8, 6,IC),IC=1,3) / +C S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 9, 5,IC),IC=1,3) / - S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/ - DATA (GB( 9, 5,IC),IC=1,3) / - S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/ - DATA (GA( 9, 6,IC),IC=1,3) / - S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/ - DATA (GB( 9, 6,IC),IC=1,3) / - S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/ +C DATA (GA( 9, 5,IC),IC=1,3) / +C S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/ +C DATA (GB( 9, 5,IC),IC=1,3) / +C S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/ +C DATA (GA( 9, 6,IC),IC=1,3) / +C S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/ +C DATA (GB( 9, 6,IC),IC=1,3) / +C S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA(10, 5,IC),IC=1,3) / - S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/ - DATA (GB(10, 5,IC),IC=1,3) / - S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/ - DATA (GA(10, 6,IC),IC=1,3) / - S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/ - DATA (GB(10, 6,IC),IC=1,3) / - S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/ +C DATA (GA(10, 5,IC),IC=1,3) / +C S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/ +C DATA (GB(10, 5,IC),IC=1,3) / +C S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/ +C DATA (GA(10, 6,IC),IC=1,3) / +C S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/ +C DATA (GB(10, 6,IC),IC=1,3) / +C S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/ C C----- INTERVAL = 5 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA(11, 5,IC),IC=1,3) / - S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/ - DATA (GB(11, 5,IC),IC=1,3) / - S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/ - DATA (GA(11, 6,IC),IC=1,3) / - S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/ - DATA (GB(11, 6,IC),IC=1,3) / - S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/ +C DATA (GA(11, 5,IC),IC=1,3) / +C S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/ +C DATA (GB(11, 5,IC),IC=1,3) / +C S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/ +C DATA (GA(11, 6,IC),IC=1,3) / +C S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/ +C DATA (GB(11, 6,IC),IC=1,3) / +C S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/ C C @@ -4332,132 +4147,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 1,11,IC),IC=1,3) / - S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/ - DATA (GB( 1,11,IC),IC=1,3) / - S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/ - DATA (GA( 1,12,IC),IC=1,3) / - S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/ - DATA (GB( 1,12,IC),IC=1,3) / - S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/ +C DATA (GA( 1,11,IC),IC=1,3) / +C S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/ +C DATA (GB( 1,11,IC),IC=1,3) / +C S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/ +C DATA (GA( 1,12,IC),IC=1,3) / +C S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/ +C DATA (GB( 1,12,IC),IC=1,3) / +C S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 2,11,IC),IC=1,3) / - S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/ - DATA (GB( 2,11,IC),IC=1,3) / - S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/ - DATA (GA( 2,12,IC),IC=1,3) / - S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/ - DATA (GB( 2,12,IC),IC=1,3) / - S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/ +C DATA (GA( 2,11,IC),IC=1,3) / +C S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/ +C DATA (GB( 2,11,IC),IC=1,3) / +C S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/ +C DATA (GA( 2,12,IC),IC=1,3) / +C S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/ +C DATA (GB( 2,12,IC),IC=1,3) / +C S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 3,11,IC),IC=1,3) / - S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/ - DATA (GB( 3,11,IC),IC=1,3) / - S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/ - DATA (GA( 3,12,IC),IC=1,3) / - S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/ - DATA (GB( 3,12,IC),IC=1,3) / - S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/ +C DATA (GA( 3,11,IC),IC=1,3) / +C S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/ +C DATA (GB( 3,11,IC),IC=1,3) / +C S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/ +C DATA (GA( 3,12,IC),IC=1,3) / +C S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/ +C DATA (GB( 3,12,IC),IC=1,3) / +C S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 4,11,IC),IC=1,3) / - S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/ - DATA (GB( 4,11,IC),IC=1,3) / - S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/ - DATA (GA( 4,12,IC),IC=1,3) / - S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/ - DATA (GB( 4,12,IC),IC=1,3) / - S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/ +C DATA (GA( 4,11,IC),IC=1,3) / +C S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/ +C DATA (GB( 4,11,IC),IC=1,3) / +C S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/ +C DATA (GA( 4,12,IC),IC=1,3) / +C S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/ +C DATA (GB( 4,12,IC),IC=1,3) / +C S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 5,11,IC),IC=1,3) / - S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/ - DATA (GB( 5,11,IC),IC=1,3) / - S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/ - DATA (GA( 5,12,IC),IC=1,3) / - S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/ - DATA (GB( 5,12,IC),IC=1,3) / - S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/ +C DATA (GA( 5,11,IC),IC=1,3) / +C S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/ +C DATA (GB( 5,11,IC),IC=1,3) / +C S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/ +C DATA (GA( 5,12,IC),IC=1,3) / +C S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/ +C DATA (GB( 5,12,IC),IC=1,3) / +C S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 6,11,IC),IC=1,3) / - S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/ - DATA (GB( 6,11,IC),IC=1,3) / - S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/ - DATA (GA( 6,12,IC),IC=1,3) / - S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/ - DATA (GB( 6,12,IC),IC=1,3) / - S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/ +C DATA (GA( 6,11,IC),IC=1,3) / +C S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/ +C DATA (GB( 6,11,IC),IC=1,3) / +C S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/ +C DATA (GA( 6,12,IC),IC=1,3) / +C S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/ +C DATA (GB( 6,12,IC),IC=1,3) / +C S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 7,11,IC),IC=1,3) / - S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/ - DATA (GB( 7,11,IC),IC=1,3) / - S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/ - DATA (GA( 7,12,IC),IC=1,3) / - S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/ - DATA (GB( 7,12,IC),IC=1,3) / - S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/ +C DATA (GA( 7,11,IC),IC=1,3) / +C S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/ +C DATA (GB( 7,11,IC),IC=1,3) / +C S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/ +C DATA (GA( 7,12,IC),IC=1,3) / +C S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/ +C DATA (GB( 7,12,IC),IC=1,3) / +C S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 8,11,IC),IC=1,3) / - S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/ - DATA (GB( 8,11,IC),IC=1,3) / - S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/ - DATA (GA( 8,12,IC),IC=1,3) / - S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/ - DATA (GB( 8,12,IC),IC=1,3) / - S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/ +C DATA (GA( 8,11,IC),IC=1,3) / +C S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/ +C DATA (GB( 8,11,IC),IC=1,3) / +C S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/ +C DATA (GA( 8,12,IC),IC=1,3) / +C S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/ +C DATA (GB( 8,12,IC),IC=1,3) / +C S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA( 9,11,IC),IC=1,3) / - S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/ - DATA (GB( 9,11,IC),IC=1,3) / - S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/ - DATA (GA( 9,12,IC),IC=1,3) / - S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/ - DATA (GB( 9,12,IC),IC=1,3) / - S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/ +C DATA (GA( 9,11,IC),IC=1,3) / +C S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/ +C DATA (GB( 9,11,IC),IC=1,3) / +C S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/ +C DATA (GA( 9,12,IC),IC=1,3) / +C S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/ +C DATA (GB( 9,12,IC),IC=1,3) / +C S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA(10,11,IC),IC=1,3) / - S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/ - DATA (GB(10,11,IC),IC=1,3) / - S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/ - DATA (GA(10,12,IC),IC=1,3) / - S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/ - DATA (GB(10,12,IC),IC=1,3) / - S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/ +C DATA (GA(10,11,IC),IC=1,3) / +C S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/ +C DATA (GB(10,11,IC),IC=1,3) / +C S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/ +C DATA (GA(10,12,IC),IC=1,3) / +C S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/ +C DATA (GB(10,12,IC),IC=1,3) / +C S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/ C C----- INTERVAL = 6 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 - DATA (GA(11,11,IC),IC=1,3) / - S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/ - DATA (GB(11,11,IC),IC=1,3) / - S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/ - DATA (GA(11,12,IC),IC=1,3) / - S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/ - DATA (GB(11,12,IC),IC=1,3) / - S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/ +C DATA (GA(11,11,IC),IC=1,3) / +C S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/ +C DATA (GB(11,11,IC),IC=1,3) / +C S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/ +C DATA (GA(11,12,IC),IC=1,3) / +C S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/ +C DATA (GB(11,12,IC),IC=1,3) / +C S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/ C C @@ -4477,132 +4292,132 @@ C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 1,13,IC),IC=1,3) / - S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/ - DATA (GB( 1,13,IC),IC=1,3) / - S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/ - DATA (GA( 1,14,IC),IC=1,3) / - S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/ - DATA (GB( 1,14,IC),IC=1,3) / - S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/ +C DATA (GA( 1,13,IC),IC=1,3) / +C S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/ +C DATA (GB( 1,13,IC),IC=1,3) / +C S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/ +C DATA (GA( 1,14,IC),IC=1,3) / +C S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/ +C DATA (GB( 1,14,IC),IC=1,3) / +C S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 2,13,IC),IC=1,3) / - S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/ - DATA (GB( 2,13,IC),IC=1,3) / - S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/ - DATA (GA( 2,14,IC),IC=1,3) / - S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/ - DATA (GB( 2,14,IC),IC=1,3) / - S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/ +C DATA (GA( 2,13,IC),IC=1,3) / +C S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/ +C DATA (GB( 2,13,IC),IC=1,3) / +C S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/ +C DATA (GA( 2,14,IC),IC=1,3) / +C S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/ +C DATA (GB( 2,14,IC),IC=1,3) / +C S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 3,13,IC),IC=1,3) / - S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/ - DATA (GB( 3,13,IC),IC=1,3) / - S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/ - DATA (GA( 3,14,IC),IC=1,3) / - S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/ - DATA (GB( 3,14,IC),IC=1,3) / - S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/ +C DATA (GA( 3,13,IC),IC=1,3) / +C S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/ +C DATA (GB( 3,13,IC),IC=1,3) / +C S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/ +C DATA (GA( 3,14,IC),IC=1,3) / +C S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/ +C DATA (GB( 3,14,IC),IC=1,3) / +C S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 4,13,IC),IC=1,3) / - S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/ - DATA (GB( 4,13,IC),IC=1,3) / - S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/ - DATA (GA( 4,14,IC),IC=1,3) / - S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/ - DATA (GB( 4,14,IC),IC=1,3) / - S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/ +C DATA (GA( 4,13,IC),IC=1,3) / +C S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/ +C DATA (GB( 4,13,IC),IC=1,3) / +C S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/ +C DATA (GA( 4,14,IC),IC=1,3) / +C S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/ +C DATA (GB( 4,14,IC),IC=1,3) / +C S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 5,13,IC),IC=1,3) / - S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/ - DATA (GB( 5,13,IC),IC=1,3) / - S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/ - DATA (GA( 5,14,IC),IC=1,3) / - S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/ - DATA (GB( 5,14,IC),IC=1,3) / - S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/ +C DATA (GA( 5,13,IC),IC=1,3) / +C S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/ +C DATA (GB( 5,13,IC),IC=1,3) / +C S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/ +C DATA (GA( 5,14,IC),IC=1,3) / +C S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/ +C DATA (GB( 5,14,IC),IC=1,3) / +C S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 6,13,IC),IC=1,3) / - S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/ - DATA (GB( 6,13,IC),IC=1,3) / - S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/ - DATA (GA( 6,14,IC),IC=1,3) / - S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/ - DATA (GB( 6,14,IC),IC=1,3) / - S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/ +C DATA (GA( 6,13,IC),IC=1,3) / +C S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/ +C DATA (GB( 6,13,IC),IC=1,3) / +C S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/ +C DATA (GA( 6,14,IC),IC=1,3) / +C S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/ +C DATA (GB( 6,14,IC),IC=1,3) / +C S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 7,13,IC),IC=1,3) / - S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/ - DATA (GB( 7,13,IC),IC=1,3) / - S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/ - DATA (GA( 7,14,IC),IC=1,3) / - S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/ - DATA (GB( 7,14,IC),IC=1,3) / - S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/ +C DATA (GA( 7,13,IC),IC=1,3) / +C S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/ +C DATA (GB( 7,13,IC),IC=1,3) / +C S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/ +C DATA (GA( 7,14,IC),IC=1,3) / +C S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/ +C DATA (GB( 7,14,IC),IC=1,3) / +C S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 8,13,IC),IC=1,3) / - S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/ - DATA (GB( 8,13,IC),IC=1,3) / - S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/ - DATA (GA( 8,14,IC),IC=1,3) / - S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/ - DATA (GB( 8,14,IC),IC=1,3) / - S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/ +C DATA (GA( 8,13,IC),IC=1,3) / +C S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/ +C DATA (GB( 8,13,IC),IC=1,3) / +C S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/ +C DATA (GA( 8,14,IC),IC=1,3) / +C S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/ +C DATA (GB( 8,14,IC),IC=1,3) / +C S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA( 9,13,IC),IC=1,3) / - S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/ - DATA (GB( 9,13,IC),IC=1,3) / - S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/ - DATA (GA( 9,14,IC),IC=1,3) / - S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/ - DATA (GB( 9,14,IC),IC=1,3) / - S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/ +C DATA (GA( 9,13,IC),IC=1,3) / +C S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/ +C DATA (GB( 9,13,IC),IC=1,3) / +C S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/ +C DATA (GA( 9,14,IC),IC=1,3) / +C S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/ +C DATA (GB( 9,14,IC),IC=1,3) / +C S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA(10,13,IC),IC=1,3) / - S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/ - DATA (GB(10,13,IC),IC=1,3) / - S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/ - DATA (GA(10,14,IC),IC=1,3) / - S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/ - DATA (GB(10,14,IC),IC=1,3) / - S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/ +C DATA (GA(10,13,IC),IC=1,3) / +C S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/ +C DATA (GB(10,13,IC),IC=1,3) / +C S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/ +C DATA (GA(10,14,IC),IC=1,3) / +C S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/ +C DATA (GB(10,14,IC),IC=1,3) / +C S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/ C C----- INTERVAL = 2 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 - DATA (GA(11,13,IC),IC=1,3) / - S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/ - DATA (GB(11,13,IC),IC=1,3) / - S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/ - DATA (GA(11,14,IC),IC=1,3) / - S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/ - DATA (GB(11,14,IC),IC=1,3) / - S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/ +C DATA (GA(11,13,IC),IC=1,3) / +C S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/ +C DATA (GB(11,13,IC),IC=1,3) / +C S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/ +C DATA (GA(11,14,IC),IC=1,3) / +C S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/ +C DATA (GB(11,14,IC),IC=1,3) / +C S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/ C C @@ -4624,147 +4439,148 @@ C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 1,15,IC),IC=1,3) / - S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/ - DATA (GB( 1,15,IC),IC=1,3) / - S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/ - DATA (GA( 1,16,IC),IC=1,3) / - S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/ - DATA (GB( 1,16,IC),IC=1,3) / - S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/ +C DATA (GA( 1,15,IC),IC=1,3) / +C S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/ +C DATA (GB( 1,15,IC),IC=1,3) / +C S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/ +C DATA (GA( 1,16,IC),IC=1,3) / +C S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/ +C DATA (GB( 1,16,IC),IC=1,3) / +C S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 200.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 2,15,IC),IC=1,3) / - S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/ - DATA (GB( 2,15,IC),IC=1,3) / - S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/ - DATA (GA( 2,16,IC),IC=1,3) / - S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/ - DATA (GB( 2,16,IC),IC=1,3) / - S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/ +C DATA (GA( 2,15,IC),IC=1,3) / +C S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/ +C DATA (GB( 2,15,IC),IC=1,3) / +C S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/ +C DATA (GA( 2,16,IC),IC=1,3) / +C S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/ +C DATA (GB( 2,16,IC),IC=1,3) / +C S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 212.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 3,15,IC),IC=1,3) / - S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/ - DATA (GB( 3,15,IC),IC=1,3) / - S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/ - DATA (GA( 3,16,IC),IC=1,3) / - S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/ - DATA (GB( 3,16,IC),IC=1,3) / - S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/ +C DATA (GA( 3,15,IC),IC=1,3) / +C S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/ +C DATA (GB( 3,15,IC),IC=1,3) / +C S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/ +C DATA (GA( 3,16,IC),IC=1,3) / +C S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/ +C DATA (GB( 3,16,IC),IC=1,3) / +C S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 225.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 4,15,IC),IC=1,3) / - S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/ - DATA (GB( 4,15,IC),IC=1,3) / - S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/ - DATA (GA( 4,16,IC),IC=1,3) / - S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/ - DATA (GB( 4,16,IC),IC=1,3) / - S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/ +C DATA (GA( 4,15,IC),IC=1,3) / +C S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/ +C DATA (GB( 4,15,IC),IC=1,3) / +C S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/ +C DATA (GA( 4,16,IC),IC=1,3) / +C S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/ +C DATA (GB( 4,16,IC),IC=1,3) / +C S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 237.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 5,15,IC),IC=1,3) / - S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/ - DATA (GB( 5,15,IC),IC=1,3) / - S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/ - DATA (GA( 5,16,IC),IC=1,3) / - S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/ - DATA (GB( 5,16,IC),IC=1,3) / - S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/ +C DATA (GA( 5,15,IC),IC=1,3) / +C S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/ +C DATA (GB( 5,15,IC),IC=1,3) / +C S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/ +C DATA (GA( 5,16,IC),IC=1,3) / +C S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/ +C DATA (GB( 5,16,IC),IC=1,3) / +C S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 250.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 6,15,IC),IC=1,3) / - S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/ - DATA (GB( 6,15,IC),IC=1,3) / - S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/ - DATA (GA( 6,16,IC),IC=1,3) / - S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/ - DATA (GB( 6,16,IC),IC=1,3) / - S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/ +C DATA (GA( 6,15,IC),IC=1,3) / +C S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/ +C DATA (GB( 6,15,IC),IC=1,3) / +C S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/ +C DATA (GA( 6,16,IC),IC=1,3) / +C S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/ +C DATA (GB( 6,16,IC),IC=1,3) / +C S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 262.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 7,15,IC),IC=1,3) / - S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/ - DATA (GB( 7,15,IC),IC=1,3) / - S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/ - DATA (GA( 7,16,IC),IC=1,3) / - S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/ - DATA (GB( 7,16,IC),IC=1,3) / - S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/ +C DATA (GA( 7,15,IC),IC=1,3) / +C S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/ +C DATA (GB( 7,15,IC),IC=1,3) / +C S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/ +C DATA (GA( 7,16,IC),IC=1,3) / +C S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/ +C DATA (GB( 7,16,IC),IC=1,3) / +C S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 275.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 8,15,IC),IC=1,3) / - S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/ - DATA (GB( 8,15,IC),IC=1,3) / - S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/ - DATA (GA( 8,16,IC),IC=1,3) / - S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/ - DATA (GB( 8,16,IC),IC=1,3) / - S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/ +C DATA (GA( 8,15,IC),IC=1,3) / +C S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/ +C DATA (GB( 8,15,IC),IC=1,3) / +C S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/ +C DATA (GA( 8,16,IC),IC=1,3) / +C S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/ +C DATA (GB( 8,16,IC),IC=1,3) / +C S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 287.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA( 9,15,IC),IC=1,3) / - S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/ - DATA (GB( 9,15,IC),IC=1,3) / - S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/ - DATA (GA( 9,16,IC),IC=1,3) / - S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/ - DATA (GB( 9,16,IC),IC=1,3) / - S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/ +C DATA (GA( 9,15,IC),IC=1,3) / +C S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/ +C DATA (GB( 9,15,IC),IC=1,3) / +C S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/ +C DATA (GA( 9,16,IC),IC=1,3) / +C S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/ +C DATA (GB( 9,16,IC),IC=1,3) / +C S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 300.0 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA(10,15,IC),IC=1,3) / - S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/ - DATA (GB(10,15,IC),IC=1,3) / - S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/ - DATA (GA(10,16,IC),IC=1,3) / - S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/ - DATA (GB(10,16,IC),IC=1,3) / - S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/ +C DATA (GA(10,15,IC),IC=1,3) / +C S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/ +C DATA (GB(10,15,IC),IC=1,3) / +C S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/ +C DATA (GA(10,16,IC),IC=1,3) / +C S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/ +C DATA (GB(10,16,IC),IC=1,3) / +C S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/ C C----- INTERVAL = 4 ----- T = 312.5 C C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 - DATA (GA(11,15,IC),IC=1,3) / - S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/ - DATA (GB(11,15,IC),IC=1,3) / - S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/ - DATA (GA(11,16,IC),IC=1,3) / - S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/ - DATA (GB(11,16,IC),IC=1,3) / - S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/ +C DATA (GA(11,15,IC),IC=1,3) / +C S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/ +C DATA (GB(11,15,IC),IC=1,3) / +C S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/ +C DATA (GA(11,16,IC),IC=1,3) / +C S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/ +C DATA (GB(11,16,IC),IC=1,3) / +C S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/ +C +C ------------------------------------------------------------------ +C DATA (( XP( J,K),J=1,6), K=1,6) / +C S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03, +C S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03, +C S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03, +C S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02, +C S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03, +C S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02, +C S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03, +C S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02, +C S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02, +C S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01, +C S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03, +C S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 / -C ------------------------------------------------------------------ - DATA (( XP( J,K),J=1,6), K=1,6) / - S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03, - S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03, - S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03, - S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02, - S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03, - S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02, - S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03, - S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02, - S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02, - S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01, - S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03, - S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 / C C @@ -4774,4 +4590,5 @@ 100 CONTINUE C +!cdir collapse DO 102 JK = 1 , KFLEV+1 DO 101 JL = 1, KDLON @@ -5625,4 +5442,5 @@ S , PADJD,PADJU,PCNTRB,PDBDT) USE dimphy + USE radiation_AR4_param, only : WG1 IMPLICIT none cym#include "dimensions.h" @@ -5685,10 +5503,5 @@ REAL(KIND=8) zwtr c -C* Data Block: -c - REAL(KIND=8) WG1(2) - SAVE WG1 -c$OMP THREADPRIVATE(WG1) - DATA (WG1(jk),jk=1,2) /1.0, 1.0/ + C----------------------------------------------------------------------- C @@ -5903,4 +5716,5 @@ C C +!cdir collapse DO 130 JA = 1 , 8 DO 120 JL = 1, KDLON @@ -6080,4 +5894,8 @@ C C + +!CDIR ON_ADB(PUU1) +!CDIR ON_ADB(PUU2) +!CDIR COLLAPSE DO 130 JA = 1 , 8 DO 120 JL = 1, KDLON Index: /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4_param.f90 =================================================================== --- /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4_param.f90 (revision 1250) +++ /LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4_param.f90 (revision 1250) @@ -0,0 +1,240 @@ +MODULE radiation_AR4_param + + REAL*8, parameter :: ZPDH2O = 0.8 + REAL*8, parameter :: ZPDUMG = 0.75 + REAL*8, parameter :: ZPRH2O = 30000.0 + REAL*8, parameter :: ZPRUMG = 30000.0 + REAL*8, parameter :: RTDH2O = 0.40 + REAL*8, parameter :: RTDUMG = 0.375 + REAL*8, parameter :: RTH2O = 240.0 + REAL*8, parameter :: RTUMG = 240.0 + + REAL*8, dimension(2), parameter :: WG1 = (/1.0, 1.0/) + REAL*8, dimension(11), parameter :: TINTP = (/ 187.5, 200., 212.5, 225., 237.5, 250., 262.5, 275., 287.5, 300., 312.5 /) + + real*8, dimension(11,16,3), parameter :: GA = reshape ( (/ & + 0.63499072E-02, 0.65566348E-02, 0.67849730E-02, 0.70481947E-02, 0.73585943E-02, 0.77242818E-02, 0.81472693E-02, 0.86227527E-02,& + 0.91396814E-02, 0.96825438E-02, 0.10233955E-01, 0.77266491E-02, 0.81323287E-02, 0.86507620E-02, 0.92776391E-02, 0.99806312E-02,& + 0.10709803E-01, 0.11414739E-01, 0.12058772E-01, 0.12623992E-01, 0.13108146E-01, 0.13518390E-01, 0.11644593E+01, 0.11747203E+01,& + 0.11837872E+01, 0.11918561E+01, 0.11990757E+01, 0.12055643E+01, 0.12114186E+01, 0.12167192E+01, 0.12215344E+01, 0.12259226E+01,& + 0.12299344E+01, 0.12006968E+01, 0.12108196E+01, 0.12196717E+01, 0.12274493E+01, 0.12343189E+01, 0.12404147E+01, 0.12458431E+01,& + 0.12506907E+01, 0.12550299E+01, 0.12589256E+01, 0.12624402E+01, 0.15750172E+00, 0.16174076E+00, 0.16548628E+00, 0.16881124E+00,& + 0.17177839E+00, 0.17443933E+00, 0.17683622E+00, 0.17900375E+00, 0.18097099E+00, 0.18276283E+00, 0.18440117E+00, 0.17770551E+00,& + 0.18176757E+00, 0.18527967E+00, 0.18833348E+00, 0.19100108E+00, 0.19334122E+00, 0.19540288E+00, 0.19722732E+00, 0.19884918E+00,& + 0.20029696E+00, 0.20159300E+00, 0.10192131E+02, 0.97258602E+01, 0.92992890E+01, 0.89154021E+01, 0.85730084E+01, 0.82685838E+01,& + 0.79978921E+01, 0.77568055E+01, 0.75416266E+01, 0.73491694E+01, 0.71767400E+01, 0.92439050E+01, 0.87567422E+01, 0.83270144E+01,& + 0.79528337E+01, 0.76286839E+01, 0.73477879E+01, 0.71035818E+01, 0.68903312E+01, 0.67032875E+01, 0.65386461E+01, 0.63934377E+01,& + 0.24870635E+02, 0.24725591E+02, 0.24600320E+02, 0.24487300E+02, 0.24384935E+02, 0.24292341E+02, 0.24208572E+02, 0.24132642E+02,& + 0.24063614E+02, 0.24000649E+02, 0.23943021E+02, 0.24586283E+02, 0.24441465E+02, 0.24311657E+02, 0.24196167E+02, 0.24093406E+02,& + 0.24001597E+02, 0.23919098E+02, 0.23844511E+02, 0.23776708E+02, 0.23714816E+02, 0.23658197E+02, 0.11990218E+02, 0.10904073E+02,& + 0.89126838E+01, 0.85622405E+01, 0.94892164E+01, 0.13580937E+02,-0.32050918E+03,-0.37133165E+01, 0.18890836E+00, 0.14209226E+01,& + 0.19817679E+01, 0.79709806E+01, 0.75400737E+01, 0.81804377E+01, 0.10564339E+02, 0.46896789E+02,-0.30926524E+01, 0.85742941E+00,& + 0.19164038E+01, 0.23513199E+01, 0.25566644E+01, 0.26555181E+01, 0.87668459E-01, 0.83754276E-01, 0.80460283E-01, 0.77659686E-01,& + 0.75257056E-01, 0.73179175E-01, 0.71369063E-01, 0.69781812E-01, 0.68381606E-01, 0.67139539E-01, 0.66032012E-01, 0.74878820E-01,& + 0.71650966E-01, 0.68979615E-01, 0.66745345E-01, 0.64857571E-01, 0.63248495E-01, 0.61866970E-01, 0.60673632E-01, 0.59637277E-01,& + 0.58732178E-01, 0.57936092E-01, 0.13230067E+02, 0.13213564E+02, 0.13209140E+02, 0.13213894E+02, 0.13225963E+02, 0.13243806E+02,& + 0.13266104E+02, 0.13291782E+02, 0.13319961E+02, 0.13349927E+02, 0.13381108E+02, 0.13183816E+02, 0.13189991E+02, 0.13209485E+02,& + 0.13238789E+02, 0.13275017E+02, 0.13316096E+02, 0.13360555E+02, 0.13407324E+02, 0.13455544E+02, 0.13504450E+02, 0.13553282E+02,& +-0.99506586E-03,-0.10184169E-02,-0.10404730E-02,-0.10621792E-02,-0.10847662E-02,-0.11094726E-02,-0.11372949E-02,-0.11687683E-02,& +-0.12038314E-02,-0.12418367E-02,-0.12817135E-02,-0.11661515E-02,-0.11886130E-02,-0.12139929E-02,-0.12445811E-02,-0.12807672E-02,& +-0.13208251E-02,-0.13619034E-02,-0.14014165E-02,-0.14378639E-02,-0.14708488E-02,-0.15006791E-02, 0.41243390E+00, 0.43407282E+00,& + 0.45331413E+00, 0.47048604E+00, 0.48586286E+00, 0.49968044E+00, 0.51214132E+00, 0.52341830E+00, 0.53365803E+00, 0.54298448E+00,& + 0.55150227E+00, 0.48318936E+00, 0.50501827E+00, 0.52409502E+00, 0.54085277E+00, 0.55565422E+00, 0.56878618E+00, 0.58047395E+00,& + 0.59089894E+00, 0.60021475E+00, 0.60856112E+00, 0.61607594E+00,-0.22159303E-01,-0.22748917E-01,-0.23269898E-01,-0.23732392E-01,& +-0.24145123E-01,-0.24515269E-01,-0.24848690E-01,-0.25150210E-01,-0.25423873E-01,-0.25673139E-01,-0.25901055E-01,-0.24972399E-01,& +-0.25537247E-01,-0.26025624E-01,-0.26450280E-01,-0.26821236E-01,-0.27146657E-01,-0.27433354E-01,-0.27687065E-01,-0.27912608E-01,& +-0.28113944E-01,-0.28294180E-01, 0.80737799E+01, 0.79171158E+01, 0.77609605E+01, 0.76087371E+01, 0.74627112E+01, 0.73239981E+01,& + 0.71929934E+01, 0.70697065E+01, 0.69539626E+01, 0.68455144E+01, 0.67441020E+01, 0.77425778E+01, 0.75443460E+01, 0.73526151E+01,& + 0.71711188E+01, 0.70015571E+01, 0.68442532E+01, 0.66987996E+01, 0.65644820E+01, 0.64405267E+01, 0.63262376E+01, 0.62210701E+01,& + 0.10542131E+02, 0.10515895E+02, 0.10492949E+02, 0.10472049E+02, 0.10452961E+02, 0.10435562E+02, 0.10419710E+02, 0.10405247E+02,& + 0.10392022E+02, 0.10379892E+02, 0.10368736E+02, 0.10490353E+02, 0.10463512E+02, 0.10439183E+02, 0.10417324E+02, 0.10397704E+02,& + 0.10380038E+02, 0.10364052E+02, 0.10349509E+02, 0.10336215E+02, 0.10324018E+02, 0.10312808E+02,-0.12823142E+01,-0.10571588E+01,& +-0.74864953E+00,-0.58705980E+00,-0.49305772E+00,-0.51461431E+00, 0.12373350E+02, 0.44809588E+00, 0.46548918E+00, 0.59121475E+00,& + 0.74676119E+00,-0.74805226E+00,-0.56252739E+00,-0.46188072E+00,-0.40712065E+00,-0.15295996E+01, 0.43555255E+00, 0.50380874E+00,& + 0.68537352E+00, 0.89437630E+00, 0.11127003E+01, 0.13329782E+01, 0.13845511E+01, 0.13187042E+01, 0.12644396E+01, 0.12191543E+01,& + 0.11809511E+01, 0.11484154E+01, 0.11204723E+01, 0.10962918E+01, 0.10752229E+01, 0.10567474E+01, 0.10404465E+01, 0.11718758E+01,& + 0.11216131E+01, 0.10809473E+01, 0.10476396E+01, 0.10200373E+01, 0.99692726E+00, 0.97740923E+00, 0.96080188E+00, 0.94657562E+00,& + 0.93430511E+00, 0.92363528E+00, 0.22042132E+02, 0.22107298E+02, 0.22180915E+02, 0.22259478E+02, 0.22341039E+02, 0.22424247E+02,& + 0.22508089E+02, 0.22591771E+02, 0.22674661E+02, 0.22756246E+02, 0.22836093E+02, 0.22169501E+02, 0.22270075E+02, 0.22379193E+02,& + 0.22492992E+02, 0.22608508E+02, 0.22723843E+02, 0.22837837E+02, 0.22949751E+02, 0.23059032E+02, 0.23165146E+02, 0.23267456E+02,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00,& + 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00, 0.00000000E+00 /) & + , (/ 11,16,3 /) ) + + real*8, dimension(11,16,3), parameter :: GB = reshape ( (/ & + 0.63499072E-02, 0.65566348E-02, 0.67849730E-02, 0.70481947E-02, 0.73585943E-02, 0.77242818E-02, 0.81472693E-02, 0.86227527E-02,& + 0.91396814E-02, 0.96825438E-02, 0.10233955E-01, 0.77266491E-02, 0.81323287E-02, 0.86507620E-02, 0.92776391E-02, 0.99806312E-02,& + 0.10709803E-01, 0.11414739E-01, 0.12058772E-01, 0.12623992E-01, 0.13108146E-01, 0.13518390E-01, 0.11644593E+01, 0.11747203E+01,& + 0.11837872E+01, 0.11918561E+01, 0.11990757E+01, 0.12055643E+01, 0.12114186E+01, 0.12167192E+01, 0.12215344E+01, 0.12259226E+01,& + 0.12299344E+01, 0.12006968E+01, 0.12108196E+01, 0.12196717E+01, 0.12274493E+01, 0.12343189E+01, 0.12404147E+01, 0.12458431E+01,& + 0.12506907E+01, 0.12550299E+01, 0.12589256E+01, 0.12624402E+01, 0.15750172E+00, 0.16174076E+00, 0.16548628E+00, 0.16881124E+00,& + 0.17177839E+00, 0.17443933E+00, 0.17683622E+00, 0.17900375E+00, 0.18097099E+00, 0.18276283E+00, 0.18440117E+00, 0.17770551E+00,& + 0.18176757E+00, 0.18527967E+00, 0.18833348E+00, 0.19100108E+00, 0.19334122E+00, 0.19540288E+00, 0.19722732E+00, 0.19884918E+00,& + 0.20029696E+00, 0.20159300E+00, 0.10192131E+02, 0.97258602E+01, 0.92992890E+01, 0.89154021E+01, 0.85730084E+01, 0.82685838E+01,& + 0.79978921E+01, 0.77568055E+01, 0.75416266E+01, 0.73491694E+01, 0.71767400E+01, 0.92439050E+01, 0.87567422E+01, 0.83270144E+01,& + 0.79528337E+01, 0.76286839E+01, 0.73477879E+01, 0.71035818E+01, 0.68903312E+01, 0.67032875E+01, 0.65386461E+01, 0.63934377E+01,& + 0.24870635E+02, 0.24725591E+02, 0.24600320E+02, 0.24487300E+02, 0.24384935E+02, 0.24292341E+02, 0.24208572E+02, 0.24132642E+02,& + 0.24063614E+02, 0.24000649E+02, 0.23943021E+02, 0.24586283E+02, 0.24441465E+02, 0.24311657E+02, 0.24196167E+02, 0.24093406E+02,& + 0.24001597E+02, 0.23919098E+02, 0.23844511E+02, 0.23776708E+02, 0.23714816E+02, 0.23658197E+02, 0.11990218E+02, 0.10904073E+02,& + 0.89126838E+01, 0.85622405E+01, 0.94892164E+01, 0.13580937E+02,-0.32050918E+03,-0.37133165E+01, 0.18890836E+00, 0.14209226E+01,& + 0.19817679E+01, 0.79709806E+01, 0.75400737E+01, 0.81804377E+01, 0.10564339E+02, 0.46896789E+02,-0.30926524E+01, 0.85742941E+00,& + 0.19164038E+01, 0.23513199E+01, 0.25566644E+01, 0.26555181E+01, 0.87668459E-01, 0.83754276E-01, 0.80460283E-01, 0.77659686E-01,& + 0.75257056E-01, 0.73179175E-01, 0.71369063E-01, 0.69781812E-01, 0.68381606E-01, 0.67139539E-01, 0.66032012E-01, 0.74878820E-01,& + 0.71650966E-01, 0.68979615E-01, 0.66745345E-01, 0.64857571E-01, 0.63248495E-01, 0.61866970E-01, 0.60673632E-01, 0.59637277E-01,& + 0.58732178E-01, 0.57936092E-01, 0.13230067E+02, 0.13213564E+02, 0.13209140E+02, 0.13213894E+02, 0.13225963E+02, 0.13243806E+02,& + 0.13266104E+02, 0.13291782E+02, 0.13319961E+02, 0.13349927E+02, 0.13381108E+02, 0.13183816E+02, 0.13189991E+02, 0.13209485E+02,& + 0.13238789E+02, 0.13275017E+02, 0.13316096E+02, 0.13360555E+02, 0.13407324E+02, 0.13455544E+02, 0.13504450E+02, 0.13553282E+02,& + 0.97222852E-01, 0.98862238E-01, 0.10061504E+00, 0.10256222E+00, 0.10475952E+00, 0.10720986E+00, 0.10985370E+00, 0.11257633E+00,& + 0.11522980E+00, 0.11766343E+00, 0.11975320E+00, 0.10681591E+00, 0.10921298E+00, 0.11198225E+00, 0.11487826E+00, 0.11751113E+00,& + 0.11951535E+00, 0.12069945E+00, 0.12108524E+00, 0.12084229E+00, 0.12019005E+00, 0.11932684E+00, 0.10346097E+01, 0.10433655E+01,& + 0.10511933E+01, 0.10582150E+01, 0.10645317E+01, 0.10702313E+01, 0.10753907E+01, 0.10800762E+01, 0.10843446E+01, 0.10882439E+01,& + 0.10918144E+01, 0.10626130E+01, 0.10716026E+01, 0.10795108E+01, 0.10865006E+01, 0.10927103E+01, 0.10982489E+01, 0.11032019E+01,& + 0.11076379E+01, 0.11116160E+01, 0.11151910E+01, 0.11184188E+01, 0.38103212E+00, 0.38913800E+00, 0.39613651E+00, 0.40222421E+00,& + 0.40756010E+00, 0.41226954E+00, 0.41645142E+00, 0.42018474E+00, 0.42353379E+00, 0.42655211E+00, 0.42928533E+00, 0.41646579E+00,& + 0.42345095E+00, 0.42937476E+00, 0.43444062E+00, 0.43880316E+00, 0.44258354E+00, 0.44587882E+00, 0.44876776E+00, 0.45131451E+00,& + 0.45357095E+00, 0.45557797E+00, 0.82623280E+01, 0.81072291E+01, 0.79523834E+01, 0.78012527E+01, 0.76561458E+01, 0.75182174E+01,& + 0.73878952E+01, 0.72652133E+01, 0.71500151E+01, 0.70420667E+01, 0.69411177E+01, 0.79342219E+01, 0.77373458E+01, 0.75467334E+01,& + 0.73661786E+01, 0.71974319E+01, 0.70408543E+01, 0.68960649E+01, 0.67623672E+01, 0.66389989E+01, 0.65252707E+01, 0.64206412E+01,& + 0.10656640E+02, 0.10630910E+02, 0.10608399E+02, 0.10587891E+02, 0.10569156E+02, 0.10552075E+02, 0.10536510E+02, 0.10522307E+02,& + 0.10509317E+02, 0.10497402E+02, 0.10486443E+02, 0.10605856E+02, 0.10579514E+02, 0.10555632E+02, 0.10534169E+02, 0.10514900E+02,& + 0.10497547E+02, 0.10481842E+02, 0.10467553E+02, 0.10454488E+02, 0.10442501E+02, 0.10431483E+02, 0.26681588E+02, 0.24728346E+02,& + 0.20551342E+02, 0.19955244E+02, 0.22227100E+02, 0.31770288E+02,-0.74061287E+03,-0.81329826E+01, 0.90279822E+00, 0.37532746E+01,& + 0.50437916E+01, 0.18377807E+02, 0.17643148E+02, 0.19296161E+02, 0.24951120E+02, 0.10957372E+03,-0.67432659E+01, 0.24550746E+01,& + 0.49089917E+01, 0.59008712E+01, 0.63532616E+01, 0.65558627E+01, 0.23203798E+01, 0.22288925E+01, 0.21515593E+01, 0.20855896E+01,& + 0.20288489E+01, 0.19796791E+01, 0.19367778E+01, 0.18991112E+01, 0.18658501E+01, 0.18363226E+01, 0.18099779E+01, 0.20206726E+01,& + 0.19441824E+01, 0.18807257E+01, 0.18275618E+01, 0.17825910E+01, 0.17442308E+01, 0.17112809E+01, 0.16828137E+01, 0.16580908E+01,& + 0.16365014E+01, 0.16175164E+01, 0.22051750E+02, 0.22116850E+02, 0.22190410E+02, 0.22268925E+02, 0.22350445E+02, 0.22433617E+02,& + 0.22517429E+02, 0.22601086E+02, 0.22683956E+02, 0.22765522E+02, 0.22845354E+02, 0.22178972E+02, 0.22279484E+02, 0.22388551E+02,& + 0.22502309E+02, 0.22617792E+02, 0.22733099E+02, 0.22847071E+02, 0.22958967E+02, 0.23068234E+02, 0.23174336E+02, 0.23276638E+02,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01,& + 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01, 0.10000000E+01 /) & + , (/ 11,16,3 /) ) + + real*8, dimension(6,6), parameter :: XP = reshape ( (/ & + 0.46430621E+02, 0.12928299E+03, 0.20732648E+03, 0.31398411E+03, 0.18373177E+03, -0.11412303E+03, & + 0.73604774E+02, 0.27887914E+03, 0.27076947E+03, -0.57322111E+02, -0.64742459E+02, 0.87238280E+02, & + 0.37050866E+02, 0.20498759E+03, 0.37558029E+03, 0.17401171E+03, -0.13350302E+03, -0.37651795E+02, & + 0.14930141E+02, 0.89161160E+02, 0.17793062E+03, 0.93433860E+02, -0.70646020E+02, -0.26373150E+02, & + 0.40386780E+02, 0.10855270E+03, 0.50755010E+02, -0.31496190E+02, 0.12791300E+00, 0.18017770E+01, & + 0.90811926E+01, 0.75073923E+02, 0.24654438E+03, 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /) , (/ 6,6 /) ) + + REAL*8, dimension(2), parameter :: RSUN = (/ 0.441676 , 0.558324 /) + REAL*8, dimension(2,6), parameter :: RRAY = reshape ( & + (/ .428937E-01, .697200E-02,& + .890743E+00, .173297E-01,& + -.288555E+01,-.850903E-01,& + .522744E+01, .248261E+00,& + -.469173E+01,-.302031E+00,& + .161645E+01, .129662E+00 /) , (/2,6/) ) + + REAL*8, dimension(2,5), parameter :: TAUA = reshape ( & + (/ 0.730719, 0.730719, 0.912819, 0.912819, 0.725059, & + 0.725059, 0.745405, 0.745405, 0.682188, 0.682188 /),(/2,5/) ) + REAL*8, dimension(2,5), parameter :: RPIZA = reshape ( & + (/ 0.872212, 0.872212, 0.982545, 0.982545, 0.623143, & + 0.623143, 0.944887, 0.944887, 0.997975, 0.997975 /),(/2,5/) ) + REAL*8, dimension(2,5), parameter :: RCGA = reshape ( & + (/ 0.647596, 0.647596, 0.739002, 0.739002, 0.580845, & + 0.580845, 0.662657, 0.662657, 0.624246, 0.624246 /),(/2,5/) ) + + REAL*8, dimension(2,3,7), parameter :: APAD = reshape ( & + (/ 0.912418292E+05, 0.376655383E-08, 0.000000000E-00,& + 0.739646016E-08, 0.925887084E-04, 0.410177786E+03,& + 0.723613782E+05, 0.978576773E-04, 0.000000000E-00,& + 0.131849595E-03, 0.129353723E-01, 0.672595424E+02,& + 0.596037057E+04, 0.387714006E+00, 0.000000000E-00,& + 0.437772681E+00, 0.800821928E+00, 0.000000000E-00,& + 0.000000000E-00, 0.118461660E+03, 0.000000000E-00,& + 0.151345118E+03, 0.242715973E+02, 0.000000000E-00,& + 0.000000000E-00, 0.119079797E+04, 0.000000000E-00,& + 0.233628890E+04, 0.878331486E+02, 0.000000000E-00,& + 0.000000000E-00, 0.293353397E+03, 0.000000000E-00,& + 0.797219934E+03, 0.191559725E+02, 0.000000000E-00,& + 0.000000000E-00, 0.000000000E+00, 0.000000000E-00,& + 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 /) , (/2,3,7/) ) + REAL*8, dimension(2,3,7), parameter :: BPAD = reshape ( & + (/ 0.912418292E+05, 0.376655383E-08, 0.000000000E-00,& + 0.739646016E-08, 0.925887084E-04, 0.410177786E+03,& + 0.724555318E+05, 0.979023421E-04, 0.000000000E-00,& + 0.131861712E-03, 0.131812683E-01, 0.731185438E+02,& + 0.602593328E+04, 0.388611139E+00, 0.000000000E-00,& + 0.437949001E+00, 0.812706117E+00, 0.100000000E+01,& + 0.100000000E+01, 0.120291383E+03, 0.000000000E-00,& + 0.151692730E+03, 0.249863591E+02, 0.000000000E+00,& + 0.000000000E-00, 0.130531005E+04, 0.000000000E-00,& + 0.237071130E+04, 0.931071925E+02, 0.000000000E+00,& + 0.000000000E-00, 0.415049409E+03, 0.000000000E-00,& + 0.867914360E+03, 0.252233437E+02, 0.000000000E+00,& + 0.000000000E-00, 0.100000000E+01, 0.000000000E-00,& + 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 /) , (/2,3,7/) ) + REAL*8, dimension(2,3), parameter :: D = reshape ( & + (/ 0.0, 0.0, 0.0, 0.0, 0.0, 0.8 /) , (/2,3/) ) + + REAL*8, parameter :: TREF = 250.0 + REAL*8, dimension(2), parameter :: RT1 = (/ -0.577350269, +0.577350269 /) + REAL*8, dimension(5,5), parameter :: RAER= reshape ( & + (/ .038520, .037196, .040532, .054934, .038520 & + , .12613 , .18313 , .10357 , .064106, .126130 & + , .012579, .013649, .018652, .025181, .012579 & + , .011890, .016142, .021105, .028908, .011890 & + , .013792, .026810, .052203, .066338, .013792 /) , (/5,5/) ) + + REAL*8, dimension(8,3), parameter :: AT= reshape ( & + (/ 0.298199E-02,0.143676E-01,0.197861E-01,0.289560E-01,& + 0.103800E-01,0.868859E-02,0.250073E-03,0.307423E-01,& + -.394023E-03,0.366501E-02,0.315541E-02,-.208807E-02,& + 0.436296E-02,-.972752E-03,0.455875E-03,0.110879E-02,& + 0.319566E-04,-.160822E-02,-.174547E-02,-.121943E-02,& + -.161431E-02,0.000000E-00,0.109242E-03,-.322172E-03 /) , (/8,3/) ) + + REAL*8, dimension(8,3), parameter :: BT= reshape ( & + (/ -0.106432E-04,-0.553979E-04,-0.877012E-04,-0.165960E-03,& + -0.276744E-04,-0.278412E-04, 0.199846E-05,-0.108482E-03,& + 0.660324E-06,-0.101701E-04, 0.513302E-04, 0.157704E-03,& + -0.327381E-04,-0.713940E-06,-0.216313E-05, 0.258096E-05,& + 0.174356E-06, 0.920868E-05, 0.523138E-06,-0.146427E-04,& + 0.127646E-04 ,0.117469E-05, 0.175991E-06,-0.814575E-06 /) , (/8,3/) ) + + REAL*8, dimension(4), parameter :: OCT = (/ -.326E-03, -.102E-05, .137E-02, -.535E-05 /) + + end module radiation_AR4_param Index: /LMDZ4/branches/LMDZ4-dev/makelmdz_fcm =================================================================== --- /LMDZ4/branches/LMDZ4-dev/makelmdz_fcm (revision 1249) +++ /LMDZ4/branches/LMDZ4-dev/makelmdz_fcm (revision 1250) @@ -162,7 +162,8 @@ # lecture des chemins propres à l'architecture de la machine # ############################################################### - -rm -f ./arch.path -ln -s ./arch/arch-${arch}.path ./arch.path +rm -f .void_file +echo > .void_file +rm -f arch.path +ln -s arch/arch-${arch}.path ./arch.path source arch.path @@ -360,4 +361,5 @@ rm -f bin/${code}${SUFF_NAME}.e rm -f arch.fcm +rm -f arch.opt echo "%ARCH $arch" >> $config_fcm @@ -380,4 +382,12 @@ ln -s arch/arch-${arch}.fcm arch.fcm +if test -f arch/arch-${arch}.opt && [ $compil_mod = "prod" ] + then + ln -s arch/arch-${arch}.opt arch.opt +else + ln -s .void_file arch.opt +fi + + rm -f $LIBOGCM/${arch}${SUFF_NAME}/.config/fcm.bld.lock ./build_gcm @@ -387,3 +397,2 @@ ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config config ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config/tmp tmp_src -#\rm -f $libf/grid/dimensions.h