Changeset 1250 for LMDZ4/branches
- Timestamp:
- Oct 29, 2009, 2:55:23 PM (15 years ago)
- Location:
- LMDZ4/branches/LMDZ4-dev
- Files:
-
- 7 added
- 1 deleted
- 16 edited
-
arch/arch-SX8_BRODIE.opt (added)
-
arch/arch-SX8_MERCURE.opt (added)
-
arch/arch-SX9_MERCURE.fcm (added)
-
arch/arch-SX9_MERCURE.opt (added)
-
arch/arch-SX9_MERCURE.path (added)
-
bld.cfg (modified) (3 diffs)
-
libf/dyn3dpar/calfis_p.F (modified) (30 diffs)
-
libf/dyn3dpar/cray.F (modified) (2 diffs)
-
libf/dyn3dpar/groupeun_p.F (modified) (4 diffs)
-
libf/dyn3dpar/leapfrog_p.F (modified) (1 diff)
-
libf/dyn3dpar/pres2lev.F (deleted)
-
libf/dyn3dpar/pres2lev.F90 (added)
-
libf/filtrez/mod_fft_wrapper.F90 (modified) (2 diffs)
-
libf/filtrez/mod_filtre_fft.F90 (modified) (2 diffs)
-
libf/phylmd/aeropt_2bands.F90 (modified) (8 diffs)
-
libf/phylmd/aeropt_5wv.F90 (modified) (6 diffs)
-
libf/phylmd/carbon_cycle_mod.F90 (modified) (3 diffs)
-
libf/phylmd/clouds_gno.F (modified) (11 diffs)
-
libf/phylmd/cvltr.F90 (modified) (3 diffs)
-
libf/phylmd/fisrtilp.F (modified) (8 diffs)
-
libf/phylmd/hbtm.F (modified) (1 diff)
-
libf/phylmd/radiation_AR4.F (modified) (33 diffs)
-
libf/phylmd/radiation_AR4_param.f90 (added)
-
makelmdz_fcm (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
LMDZ4/branches/LMDZ4-dev/bld.cfg
r1228 r1250 27 27 src::filtrez %SRC_PATH/filtrez 28 28 src::bibio %SRC_PATH/bibio 29 src::ext_src %EXT_SRC30 29 31 30 bld::lib::dyn %DYN … … 39 38 bld::target lib%{DYN}.a lib%{PHYS}.a libgrid.a libfiltrez.a libbibio.a 40 39 bld::target %EXEC%SUFF_NAME.e 41 bld::exe_dep %{DYN} %{PHYS} grid filtrez bibio ext_src40 bld::exe_dep %{DYN} %{PHYS} grid filtrez bibio 42 41 43 42 … … 60 59 bld::tool::fppkeys %CPP_KEY %FPP_DEF 61 60 61 62 #bld::tool::fflags::phys::readaerosol %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt -pi auto 63 #bld::tool::fflags::phys::aeropt_2bands %BASE_FFLAGS %PROD_FFLAGS %INCDIR 64 #bld::tool::fflags::phys::radiation_AR4 %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -Wf,-O,extendreorder 65 #bld::tool::fflags::phys::radiation_AR4_param %BASE_FFLAGS %PROD_FFLAGS1 %INCDIR -C hopt -f3 66 #bld::tool::fflags::phys::fisrtilp %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt 67 #bld::tool::fflags::phys::cv30_routines %BASE_FFLAGS %PROD_FFLAGS %INCDIR -Wf,-O,extendreorder 68 #bld::tool::fflags::phys::cvltr %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt 69 #bld::tool::fflags::phys::clouds_gno %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt 70 #bld::tool::fflags::dyn::vlsplt_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt 71 #bld::tool::fflags::dyn::groupeun_p %BASE_FFLAGS %PROD_FFLAGS %INCDIR -C hopt 72 73 74 inc arch.opt 62 75 63 76 # Pre-process code before analysing dependencies -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/calfis_p.F
r1231 r1250 158 158 REAL,ALLOCATABLE,SAVE :: pcvgt(:,:), pcvgq(:,:,:) 159 159 c 160 c REAL,ALLOCATABLE,SAVE :: pvervel(:,:)161 c162 160 REAL,ALLOCATABLE,SAVE :: zdufi(:,:),zdvfi(:,:) 163 161 REAL,ALLOCATABLE,SAVE :: zdtfi(:,:),zdqfi(:,:,:) … … 175 173 REAL,ALLOCATABLE,SAVE :: ztfi_omp(:,:) 176 174 REAL,ALLOCATABLE,SAVE :: zqfi_omp(:,:,:) 177 c REAL,ALLOCATABLE,SAVE :: pvervel_omp(:,:)178 175 REAL,ALLOCATABLE,SAVE :: zdufi_omp(:,:) 179 176 REAL,ALLOCATABLE,SAVE :: zdvfi_omp(:,:) … … 252 249 ALLOCATE(pcvgu(klon,llm), pcvgv(klon,llm)) 253 250 ALLOCATE(pcvgt(klon,llm), pcvgq(klon,llm,2)) 254 c ALLOCATE(pvervel(klon,llm))255 251 ALLOCATE(zdufi(klon,llm),zdvfi(klon,llm)) 256 252 ALLOCATE(zdtfi(klon,llm),zdqfi(klon,llm,nqtot)) … … 278 274 279 275 c$OMP MASTER 276 !CDIR ON_ADB(index_i) 277 !CDIR ON_ADB(index_j) 280 278 do ig0=1,klon 281 279 i=index_i(ig0) … … 300 298 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 301 299 DO l = 1, llmp1 300 !CDIR ON_ADB(index_i) 301 !CDIR ON_ADB(index_j) 302 302 do ig0=1,klon 303 303 i=index_i(ig0) … … 314 314 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 315 315 DO l=1,llm 316 316 !CDIR ON_ADB(index_i) 317 !CDIR ON_ADB(index_j) 317 318 do ig0=1,klon 318 319 i=index_i(ig0) … … 321 322 zplay(ig0,l) = preff * pksurcp ** unskap 322 323 ztfi(ig0,l) = pteta(i,j,l) * pksurcp 323 c pcvgt(ig0,l) = pdteta(i,j,l) * pksurcp / pmasse(i,j,l)324 324 enddo 325 325 … … 335 335 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 336 336 DO l=1,llm 337 !CDIR ON_ADB(index_i) 338 !CDIR ON_ADB(index_j) 337 339 do ig0=1,klon 338 340 i=index_i(ig0) … … 344 346 ENDDO 345 347 346 c convergence dynamique pour les traceurs "EAU"347 ! Earth-specific treatment of first 2 tracers (water)348 if (planet_type=="earth") then349 DO iq=1,2350 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)351 DO l=1,llm352 do ig0=1,klon353 i=index_i(ig0)354 j=index_j(ig0)355 c pcvgq(ig0,l,iq) = pdq(i,j,l,iq) / pmasse(i,j,l)356 enddo357 ENDDO358 c$OMP END DO NOWAIT359 ENDDO360 endif ! of if (planet_type=="earth")361 362 348 363 349 c Geopotentiel calcule par rapport a la surface locale: … … 378 364 c$OMP END DO NOWAIT 379 365 380 c .... Calcul de la vitesse verticale ( en Pa*m*s ou Kg/s ) ....381 c JG : ancien calcule de omega utilise dans physiq.F. Maintenant le flux382 c de masse est calclue dans advtrac_p.F383 c384 cc$OMP DO SCHEDULE(STATIC,OMP_CHUNK)385 c DO l=1,llm386 c do ig0=1,klon387 c i=index_i(ig0)388 c j=index_j(ig0)389 c pvervel(ig0,l) = pw(i,j,l)*g* unsaire(i,j)390 c enddo391 c if (is_north_pole) pvervel(1,l)=pw(1,1,l)*g /apoln392 c if (is_south_pole) pvervel(klon,l)=pw(1,jjp1,l)*g/apols393 c ENDDO394 cc$OMP END DO NOWAIT395 366 396 367 c … … 406 377 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 407 378 DO l=1,llm 379 !CDIR ON_ADB(index_i) 380 !CDIR ON_ADB(index_j) 381 !CDIR SPARSE 408 382 do ig0=kstart,kend 409 383 i=index_i(ig0) … … 412 386 zufi(ig0,l)= 0.5 *( pucov(iim,j,l)/cu(iim,j) 413 387 $ + pucov(1,j,l)/cu(1,j) ) 414 c pcvgu(ig0,l)= 0.5*( pducov(iim,j,l)/cu(iim,j)415 c $ + pducov(1,j,l)/cu(1,j) )416 388 else 417 389 zufi(ig0,l)= 0.5*( pucov(i-1,j,l)/cu(i-1,j) 418 390 $ + pucov(i,j,l)/cu(i,j) ) 419 c pcvgu(ig0,l)= 0.5*( pducov(i-1,j,l)/cu(i-1,j)420 c $ + pducov(i,j,l)/cu(i,j) )421 391 endif 422 392 enddo 423 393 ENDDO 424 394 c$OMP END DO NOWAIT 395 425 396 c 46.champ v: 426 397 c ----------- 398 427 399 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 428 400 DO l=1,llm 401 !CDIR ON_ADB(index_i) 402 !CDIR ON_ADB(index_j) 429 403 DO ig0=kstart,kend 430 404 i=index_i(ig0) … … 433 407 $ + pvcov(i,j,l)/cv(i,j) ) 434 408 435 c pcvgv(ig0+i,l)= 0.5 * ( pdvcov(i,j-1,l)/cv(i,j-1)436 c $ + pdvcov(i,j,l)/cv(i,j) )437 409 ENDDO 438 410 ENDDO … … 449 421 450 422 z1(1) =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,1,l)/cv(1,1) 451 c z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,1,l)/cv(1,1)452 423 DO i=2,iim 453 424 z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,1,l)/cv(i,1) 454 c z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,1,l)/cv(i,1)455 425 ENDDO 456 426 457 427 DO i=1,iim 458 428 zcos(i) = COS(rlonv(i))*z1(i) 459 c zcosbis(i)= COS(rlonv(i))*z1bis(i)460 429 zsin(i) = SIN(rlonv(i))*z1(i) 461 c zsinbis(i)= SIN(rlonv(i))*z1bis(i)462 430 ENDDO 463 431 464 432 zufi(1,l) = SSUM(iim,zcos,1)/pi 465 c pcvgu(1,l) = SSUM(iim,zcosbis,1)/pi466 433 zvfi(1,l) = SSUM(iim,zsin,1)/pi 467 c pcvgv(1,l) = SSUM(iim,zsinbis,1)/pi468 434 469 435 ENDDO … … 482 448 483 449 z1(1) =(rlonu(1)-rlonu(iim)+2.*pi)*pvcov(1,jjm,l)/cv(1,jjm) 484 c z1bis(1)=(rlonu(1)-rlonu(iim)+2.*pi)*pdvcov(1,jjm,l)/cv(1,jjm)485 450 DO i=2,iim 486 z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm) 487 c z1bis(i)=(rlonu(i)-rlonu(i-1))*pdvcov(i,jjm,l)/cv(i,jjm) 451 z1(i) =(rlonu(i)-rlonu(i-1))*pvcov(i,jjm,l)/cv(i,jjm) 488 452 ENDDO 489 453 490 454 DO i=1,iim 491 455 zcos(i) = COS(rlonv(i))*z1(i) 492 c zcosbis(i) = COS(rlonv(i))*z1bis(i)493 456 zsin(i) = SIN(rlonv(i))*z1(i) 494 c zsinbis(i) = SIN(rlonv(i))*z1bis(i)495 457 ENDDO 496 458 497 459 zufi(klon,l) = SSUM(iim,zcos,1)/pi 498 c pcvgu(klon,l) = SSUM(iim,zcosbis,1)/pi499 460 zvfi(klon,l) = SSUM(iim,zsin,1)/pi 500 c pcvgv(klon,l) = SSUM(iim,zsinbis,1)/pi501 502 461 ENDDO 503 462 c$OMP END DO NOWAIT … … 506 465 507 466 IF (is_sequential) THEN 508 if (planet_type=="earth") then509 #ifdef CPP_EARTH510 467 c 511 468 cIM calcul PV a teta=350, 380, 405K … … 514 471 $ ntetaSTD,rtetaSTD,PVteta) 515 472 c 516 #endif517 endif518 473 ENDIF 519 474 … … 525 480 c --------------------- 526 481 527 cc$OMP PARALLEL DEFAULT(NONE)528 cc$OMP+ PRIVATE(i,l,offset,iq)529 cc$OMP+ SHARED(klon_omp_nb,nqtot,klon_omp_begin,530 cc$OMP+ debut,lafin,rdayvrai,heure,dtphys,zplev,zplay,531 cc$OMP+ zphi,zphis,presnivs,clesphy0,zufi,zvfi,ztfi,532 cc$OMP+ zqfi,pvervel,zdufi,zdvfi,zdtfi,zdqfi,zdpsrf)533 534 c PRIVATE(zplev_omp,zplay_omp,zphi_omp,zphis_omp,535 c c$OMP+ presnivs_omp,zufi_omp,zvfi_omp,ztfi_omp,536 c c$OMP+ zqfi_omp,pvervel_omp,zdufi_omp,zdvfi_omp,537 c c$OMP+ zdtfi_omp,zdqfi_omp,zdpsrf_omp)538 482 539 483 c$OMP BARRIER … … 550 494 allocate(ztfi_omp(klon,llm)) 551 495 allocate(zqfi_omp(klon,llm,nqtot)) 552 c allocate(pvervel_omp(klon,llm))553 496 allocate(zdufi_omp(klon,llm)) 554 497 allocate(zdvfi_omp(klon,llm)) … … 617 560 enddo 618 561 619 c do l=1,llm620 c do i=1,klon621 c pvervel_omp(i,l)=pvervel(offset+i,l)622 c enddo623 c enddo624 625 562 do l=1,llm 626 563 do i=1,klon … … 660 597 661 598 c$OMP BARRIER 662 cym call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm)663 599 664 600 if (planet_type=="earth") then … … 681 617 . ztfi_omp, 682 618 . zqfi_omp, 683 c . pvervel_omp,684 619 c#ifdef INCA 685 620 . flxwfi_omp, … … 695 630 #endif 696 631 endif !of if (planet_type=="earth") 697 698 cym call WriteField_phy_p('zdtfi_omp',zdtfi_omp(:,:),llm)699 700 632 c$OMP BARRIER 701 633 … … 753 685 enddo 754 686 enddo 755 756 c do l=1,llm757 c do i=1,klon758 c pvervel(offset+i,l)=pvervel_omp(i,l)759 c enddo760 c enddo761 687 762 688 do l=1,llm … … 791 717 792 718 793 cc$OMP END PARALLEL794 719 klon=klon_mpi 795 720 500 CONTINUE … … 797 722 798 723 c$OMP MASTER 799 cym call WriteField_phy('zdtfi',zdtfi(:,:),llm)800 724 call stop_timer(timer_physic) 801 725 c$OMP END MASTER … … 913 837 DO l=1,llm 914 838 915 !!cdir NODEP 839 !CDIR ON_ADB(index_i) 840 !CDIR ON_ADB(index_j) 841 !cdir NODEP 916 842 do ig0=kstart,kend 917 843 i=index_i(ig0) … … 975 901 976 902 C 977 903 !cdir NODEP 978 904 DO iq=1,nqtot 979 905 iiq=niadv(iq) 980 906 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 981 907 DO l=1,llm 982 983 !!cdir NODEP 908 !CDIR ON_ADB(index_i) 909 !CDIR ON_ADB(index_j) 910 !cdir NODEP 984 911 DO ig0=kstart,kend 985 912 i=index_i(ig0) … … 1009 936 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 1010 937 DO l=1,llm 1011 !!cdir NODEP 938 !CDIR ON_ADB(index_i) 939 !CDIR ON_ADB(index_j) 940 !cdir NODEP 1012 941 do ig0=kstart,kend 1013 942 i=index_i(ig0) … … 1052 981 c$OMP DO SCHEDULE(STATIC,OMP_CHUNK) 1053 982 DO l=1,llm 1054 !!cdir NODEP 983 !CDIR ON_ADB(index_i) 984 !CDIR ON_ADB(index_j) 985 !cdir NODEP 1055 986 do ig0=kstart,kend 1056 987 i=index_i(ig0) -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/cray.F
r774 r1250 13 13 real sx((n-1)*incx+1),sy((n-1)*incy+1) 14 14 c 15 if (incx.eq.1.and.incy.eq.1) then 16 do 10 i=1,n 17 sy(i)=sx(i) 18 10 continue 19 else 15 20 iy=1 16 21 ix=1 17 do 1 0i=1,n22 do 11 i=1,n 18 23 sy(iy)=sx(ix) 19 24 ix=ix+incx 20 25 iy=iy+incy 21 10 continue 26 11 continue 27 endif 22 28 c 23 29 return … … 32 38 c 33 39 ssum=0. 40 if (incx.eq.1) then 41 do 10 i=1,n 42 ssum=ssum+sx(i) 43 10 continue 44 else 34 45 ix=1 35 do 1 0i=1,n46 do 11 i=1,n 36 47 ssum=ssum+sx(ix) 37 48 ix=ix+incx 38 10 continue 49 11 continue 50 endif 39 51 c 40 52 return -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/groupeun_p.F
r1087 r1250 1 1 SUBROUTINE groupeun_p(jjmax,llmax,jjb,jje,q) 2 2 USE parallel 3 USE Write_Field_p 3 4 IMPLICIT NONE 4 5 … … 17 18 REAL airecs,qs 18 19 19 INTEGER i,j,l,ig, j1,j2,i0,jd20 INTEGER i,j,l,ig,ig2,j1,j2,i0,jd 20 21 21 22 c--------------------------------------------------------------------c … … 37 38 LOGICAL, SAVE :: first = .TRUE. 38 39 !$OMP THREADPRIVATE(first) 40 INTEGER,SAVE :: i_index(iim,ngroup) 41 INTEGER :: offset 42 REAL :: qsum(iim/ngroup) 39 43 40 44 IF (first) THEN … … 54 58 j_start = MAX(jjb, j1-jd) 55 59 j_finish = MIN(jje, j2-jd) 56 DO j=j_start, j_finish 57 DO i0=1,iim,2**(ngroup-ig+1) 58 qn=0. 59 DO i=i0,i0+2**(ngroup-ig+1)-1 60 qn=qn+q(i,j,l) 61 ENDDO 62 DO i=i0,i0+2**(ngroup-ig+1)-1 63 q(i,j,l)=qn*airen_tab(i,j,jd) 64 ENDDO 60 DO ig2=1,ngroup-ig+1 61 offset=2**(ig2-1) 62 DO j=j_start, j_finish 63 !CDIR NODEP 64 !CDIR ON_ADB(q) 65 DO i0=1,iim,2**ig2 66 q(i0,j,l)=q(i0,j,l)+q(i0+offset,j,l) 67 ENDDO 68 ENDDO 69 ENDDO 70 71 DO j=j_start, j_finish 72 !CDIR NODEP 73 !CDIR ON_ADB(q) 74 DO i=1,iim 75 q(i,j,l)=q(i-MOD(i-1,2**(ngroup-ig+1)),j,l) 76 ENDDO 77 ENDDO 78 79 DO j=j_start, j_finish 80 !CDIR ON_ADB(airen_tab) 81 !CDIR ON_ADB(q) 82 DO i=1,iim 83 q(i,j,l)=q(i,j,l)*airen_tab(i,j,jd) 65 84 ENDDO 66 85 q(iip1,j,l)=q(1,j,l) 67 86 ENDDO 68 87 69 88 !c Concerne le pole sud 70 89 j_start = MAX(1+jjp1-jje-jd, j1-jd) 71 90 j_finish = MIN(1+jjp1-jjb-jd, j2-jd) 72 DO j=j_start, j_finish 73 DO i0=1,iim,2**(ngroup-ig+1) 74 qs=0. 75 DO i=i0,i0+2**(ngroup-ig+1)-1 76 qs=qs+q(i,jjp1-j+1-jd,l) 77 ENDDO 78 DO i=i0,i0+2**(ngroup-ig+1)-1 79 q(i,jjp1-j+1-jd,l)=qs*aires_tab(i,jjp1-j+1,jd) 80 ENDDO 91 DO ig2=1,ngroup-ig+1 92 offset=2**(ig2-1) 93 DO j=j_start, j_finish 94 !CDIR NODEP 95 !CDIR ON_ADB(q) 96 DO i0=1,iim,2**ig2 97 q(i0,jjp1-j+1-jd,l)= q(i0,jjp1-j+1-jd,l) 98 & +q(i0+offset,jjp1-j+1-jd,l) 99 ENDDO 100 ENDDO 101 ENDDO 102 103 104 DO j=j_start, j_finish 105 !CDIR NODEP 106 !CDIR ON_ADB(q) 107 DO i=1,iim 108 q(i,jjp1-j+1-jd,l)=q(i-MOD(i-1,2**(ngroup-ig+1)), 109 & jjp1-j+1-jd,l) 110 ENDDO 111 ENDDO 112 113 DO j=j_start, j_finish 114 !CDIR ON_ADB(aires_tab) 115 !CDIR ON_ADB(q) 116 DO i=1,iim 117 q(i,jjp1-j+1-jd,l)=q(i,jjp1-j+1-jd,l)* 118 & aires_tab(i,jjp1-j+1,jd) 81 119 ENDDO 82 120 q(iip1,jjp1-j+1-jd,l)=q(1,jjp1-j+1-jd,l) 83 121 ENDDO 122 84 123 85 124 j1=j2+1 -
LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/leapfrog_p.F
r1247 r1250 327 327 c$OMP MASTER 328 328 ItCount=ItCount+1 329 if (MOD(ItCount,1 )==1) then329 if (MOD(ItCount,12)==0) then 330 330 debug=.true. 331 331 else -
LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_fft_wrapper.F90
r986 r1250 19 19 INTEGER,INTENT(IN) :: nb_vect 20 20 REAL,INTENT(IN) :: vect(vsize+inc,nb_vect) 21 COMPLEX ,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect)21 COMPLEX*16,INTENT(INOUT) :: TF_vect(vsize/2+1,nb_vect) 22 22 23 23 STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" … … 29 29 INTEGER,INTENT(IN) :: nb_vect 30 30 REAL,INTENT(INOUT) :: vect(vsize+inc,nb_vect) 31 COMPLEX ,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect)31 COMPLEX*16,INTENT(IN ) :: TF_vect(vsize/2+1,nb_vect) 32 32 33 33 STOP "wrapper fft : une FFT doit etre specifiee a l'aide d'une clee CPP, sinon utiliser le filtre classique" -
LMDZ4/branches/LMDZ4-dev/libf/filtrez/mod_filtre_fft.F90
r1249 r1250 199 199 200 200 REAL :: vect(iim+inc,jj_end-jj_begin+1,nbniv) 201 COMPLEX 201 COMPLEX*16 :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) 202 202 INTEGER :: nb_vect 203 203 INTEGER :: i,j,l … … 260 260 REAL,INTENT(INOUT) :: vect_inout(iim+1,nlat,nbniv) 261 261 262 263 COMPLEX 262 REAL :: vect(iim+inc,jj_end-jj_begin+1,nbniv) 263 COMPLEX*16 :: TF_vect(iim/2+1,jj_end-jj_begin+1,nbniv) 264 264 INTEGER :: nb_vect 265 265 INTEGER :: i,j,l -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_2bands.F90
r1249 r1250 51 51 REAL, DIMENSION(klon,klev,naero_tot,nbands) :: cg_ae 52 52 LOGICAL :: soluble 53 INTEGER :: i, k, ierr, inu, m, mrfspecies53 INTEGER :: i, k,n, ierr, inu, m, mrfspecies 54 54 INTEGER :: spsol, spinsol, spss 55 INTEGER :: RH_num 55 INTEGER :: RH_num(klon,klev) 56 56 INTEGER, PARAMETER :: nb_level=19 ! number of vertical levels in DATA 57 57 … … 122 122 REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./) 123 123 REAL, PARAMETER :: RH_MAX=95. 124 REAL:: DELTA, rh, H 125 REAL:: tau_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of epaisseur optique aerosol 126 REAL:: piz_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of Single scattering albedo 127 REAL:: cg_ae2b_int(KLON,KLEV,nbands) ! Intermediate computation of Assymetry parameter 124 REAL:: DELTA(klon,klev), rh(klon,klev), H 125 REAL:: tau_ae2b_int ! Intermediate computation of epaisseur optique aerosol 126 REAL:: piz_ae2b_int ! Intermediate computation of Single scattering albedo 127 REAL:: cg_ae2b_int ! Intermediate computation of Assymetry parameter 128 REAL :: Fact_RH(nbre_RH) 128 129 REAL :: zrho 129 130 REAL :: fac … … 146 147 REAL:: piz_aeri_2bands(nbands,naero_insoluble) !-- unit 147 148 149 INTEGER :: id 150 LOGICAL :: used_aer(naero_tot) 151 REAL :: tmp_var, tmp_var_pi 148 152 149 153 DATA presnivs_19/& … … 586 590 587 591 DO k=1, klev 588 DO i=1, klon 589 ! IF (t_seri(i,k).EQ.0.) THEN 590 ! WRITE(lunout,*) 't_seri(i,k)=0 for i=',i,'k=',k 591 ! CALL abort_gcm('aeropt_2bands','t_seri=0',1) 592 ! END IF 593 ! IF (pplay(i,k).EQ.0.) THEN 594 ! WRITE(lunout,*) 'pplay(i,k)=0 for i=',i,'k=',k 595 ! CALL abort_gcm('aeropt_2bands','pplay=0',1) 596 ! END IF 597 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 598 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 592 DO i=1, klon 593 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 594 !CDIR UNROLL=naero_spc 595 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 599 596 !RAF zrho 600 mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 601 ENDDO 597 !CDIR UNROLL=naero_spc 598 mass_temp_pi(i,k,:) = m_allaer_pi(i,k,:) / zrho / 1.e+9 599 zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] 600 ENDDO 602 601 ENDDO 603 602 604 603 IF (flag_aerosol .EQ. 1) THEN 605 nb_aer = 1604 nb_aer = 2 606 605 ALLOCATE (aerosol_name(nb_aer)) 607 606 aerosol_name(1) = id_ASSO4M 608 607 aerosol_name(2) = id_CSSO4M 609 608 ELSEIF (flag_aerosol .EQ. 2) THEN 610 609 nb_aer = 2 … … 650 649 ! compute optical_thickness_at_gridpoint_per_species 651 650 652 tau_ae(:,:,:,:)=0. 653 !RAF 654 tau_ae_pi(:,:,:,:)=0. 655 piz_ae(:,:,:,:)=0. 656 cg_ae(:,:,:,:)=0. 657 tau_allaer(:,:,:,:)=0. 658 piz_allaer(:,:,:,:)=0. 659 cg_allaer(:,:,:,:)=0. 660 661 ! 662 ! Calculations that need to be done since we are not in the subroutines INCA 663 ! 664 ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] 665 zdp1(:,:)=pdel(:,:)/(gravit*delt) 666 667 651 652 653 !!CDIR ON_ADB(RH_tab) 654 !CDIR ON_ADB(fact_RH) 655 !CDIR SHORTLOOP 656 DO n=1,nbre_RH-1 657 fact_RH(n)=1./(RH_tab(n+1)-RH_tab(n)) 658 ENDDO 659 660 DO k=1, KLEV 661 !!CDIR ON_ADB(RH_tab) 662 !CDIR ON_ADB(fact_RH) 663 DO i=1, KLON 664 rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) 665 RH_num(i,k) = INT( rh(i,k)/10. + 1.) 666 IF (rh(i,k).GT.85.) RH_num(i,k)=10 667 IF (rh(i,k).GT.90.) RH_num(i,k)=11 668 669 DELTA(i,k)=(rh(i,k)-RH_tab(RH_num(i,k)))*fact_RH(RH_num(i,k)) 670 ENDDO 671 ENDDO 672 673 used_aer(:)=.FALSE. 674 668 675 DO m=1,nb_aer ! tau is only computed for each mass 669 670 fac=1.0 676 fac=1.0 671 677 IF (aerosol_name(m).EQ.id_ASBCM) THEN 672 678 soluble=.TRUE. … … 715 721 ENDIF 716 722 717 718 tau_ae2b_int(:,:,:)=0. 719 piz_ae2b_int(:,:,:)=0. 720 cg_ae2b_int(:,:,:)=0. 721 722 DO inu=1,nbands 723 id=aerosol_name(m) 724 used_aer(id)=.TRUE. 725 726 727 IF (soluble) THEN 728 729 IF (spss.NE.0) THEN 730 731 IF (spss.EQ.1) THEN !accumulation mode 732 DO k=1, KLEV 733 !CDIR ON_ADB(A1_ASSSM_b1) 734 !CDIR ON_ADB(A2_ASSSM_b1) 735 !CDIR ON_ADB(A3_ASSSM_b1) 736 !CDIR ON_ADB(B1_ASSSM_b1) 737 !CDIR ON_ADB(B2_ASSSM_b1) 738 !CDIR ON_ADB(C1_ASSSM_b1) 739 !CDIR ON_ADB(C2_ASSSM_b2) 740 !CDIR ON_ADB(A1_ASSSM_b2) 741 !CDIR ON_ADB(A2_ASSSM_b2) 742 !CDIR ON_ADB(A3_ASSSM_b2) 743 !CDIR ON_ADB(B1_ASSSM_b2) 744 !CDIR ON_ADB(B2_ASSSM_b2) 745 !CDIR ON_ADB(C1_ASSSM_b2) 746 !CDIR ON_ADB(C2_ASSSM_b2) 747 DO i=1, KLON 748 H=rh(i,k)/100 749 tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 750 tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 751 752 ! band 1 753 tau_ae2b_int=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05) 754 piz_ae2b_int=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H 755 cg_ae2b_int=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H 756 757 tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int 758 tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int 759 piz_ae(i,k,id,1) = piz_ae2b_int 760 cg_ae(i,k,id,1)= cg_ae2b_int 761 762 !band 2 763 tau_ae2b_int=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05) 764 piz_ae2b_int=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H 765 cg_ae2b_int=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H 766 767 tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int 768 tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int 769 piz_ae(i,k,id,2) = piz_ae2b_int 770 cg_ae(i,k,id,2)= cg_ae2b_int 771 772 ENDDO 773 ENDDO 774 ENDIF 775 776 IF (spss.EQ.2) THEN !coarse mode 777 DO k=1, KLEV 778 !CDIR ON_ADB(A1_CSSSM_b1) 779 !CDIR ON_ADB(A2_CSSSM_b1) 780 !CDIR ON_ADB(A3_CSSSM_b1) 781 !CDIR ON_ADB(B1_CSSSM_b1) 782 !CDIR ON_ADB(B2_CSSSM_b1) 783 !CDIR ON_ADB(C1_CSSSM_b1) 784 !CDIR ON_ADB(C2_CSSSM_b2) 785 !CDIR ON_ADB(A1_CSSSM_b2) 786 !CDIR ON_ADB(A2_CSSSM_b2) 787 !CDIR ON_ADB(A3_CSSSM_b2) 788 !CDIR ON_ADB(B1_CSSSM_b2) 789 !CDIR ON_ADB(B2_CSSSM_b2) 790 !CDIR ON_ADB(C1_CSSSM_b2) 791 !CDIR ON_ADB(C2_CSSSM_b2) 792 DO i=1, KLON 793 H=rh(i,k)/100 794 tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 795 tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 796 ! band 1 797 tau_ae2b_int=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05) 798 piz_ae2b_int=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H 799 cg_ae2b_int=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H 800 801 tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int 802 tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int 803 piz_ae(i,k,id,1) = piz_ae2b_int 804 cg_ae(i,k,id,1)= cg_ae2b_int 805 806 ! band 2 807 tau_ae2b_int=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05) 808 piz_ae2b_int=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H 809 cg_ae2b_int=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H 810 811 tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int 812 tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int 813 piz_ae(i,k,id,2) = piz_ae2b_int 814 cg_ae(i,k,id,2)= cg_ae2b_int 815 816 ENDDO 817 ENDDO 818 ENDIF 819 820 IF (spss.EQ.3) THEN !super coarse mode 821 DO k=1, KLEV 822 !CDIR ON_ADB(A1_SSSSM_b1) 823 !CDIR ON_ADB(A2_SSSSM_b1) 824 !CDIR ON_ADB(A3_SSSSM_b1) 825 !CDIR ON_ADB(B1_SSSSM_b1) 826 !CDIR ON_ADB(B2_SSSSM_b1) 827 !CDIR ON_ADB(C1_SSSSM_b1) 828 !CDIR ON_ADB(C2_SSSSM_b2) 829 !CDIR ON_ADB(A1_SSSSM_b2) 830 !CDIR ON_ADB(A2_SSSSM_b2) 831 !CDIR ON_ADB(A3_SSSSM_b2) 832 !CDIR ON_ADB(B1_SSSSM_b2) 833 !CDIR ON_ADB(B2_SSSSM_b2) 834 !CDIR ON_ADB(C1_SSSSM_b2) 835 !CDIR ON_ADB(C2_SSSSM_b2) 836 DO i=1, KLON 837 H=rh(i,k)/100 838 tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 839 tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 840 841 ! band 1 842 tau_ae2b_int=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05) 843 piz_ae2b_int=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H 844 cg_ae2b_int=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H 845 846 tau_ae(i,k,id,1) = tmp_var*tau_ae2b_int 847 tau_ae_pi(i,k,id,1) = tmp_var_pi* tau_ae2b_int 848 piz_ae(i,k,id,1) = piz_ae2b_int 849 cg_ae(i,k,id,1)= cg_ae2b_int 850 851 ! band 2 852 tau_ae2b_int=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05) 853 piz_ae2b_int=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H 854 cg_ae2b_int=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H 855 856 tau_ae(i,k,id,2) = tmp_var*tau_ae2b_int 857 tau_ae_pi(i,k,id,2) = tmp_var_pi* tau_ae2b_int 858 piz_ae(i,k,id,2) = piz_ae2b_int 859 cg_ae(i,k,id,2)= cg_ae2b_int 860 861 ENDDO 862 ENDDO 863 ENDIF 864 865 ELSE 866 867 !CDIR ON_ADB(alpha_aers_2bands) 868 !CDIR ON_ADB(piz_aers_2bands) 869 !CDIR ON_ADB(cg_aers_2bands) 870 DO k=1, KLEV 871 DO i=1, KLON 872 tmp_var=mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 873 tmp_var_pi=mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt*fac 874 !CDIR UNROLL=nbands 875 DO inu=1,nbands 876 877 tau_ae2b_int= alpha_aers_2bands(RH_num(i,k),inu,spsol)+ & 878 DELTA(i,k)* (alpha_aers_2bands(RH_num(i,k)+1,inu,spsol) - & 879 alpha_aers_2bands(RH_num(i,k),inu,spsol)) 880 881 piz_ae2b_int = piz_aers_2bands(RH_num(i,k),inu,spsol) + & 882 DELTA(i,k)* (piz_aers_2bands(RH_num(i,k)+1,inu,spsol) - & 883 piz_aers_2bands(RH_num(i,k),inu,spsol)) 884 885 cg_ae2b_int = cg_aers_2bands(RH_num(i,k),inu,spsol) + & 886 DELTA(i,k)* (cg_aers_2bands(RH_num(i,k)+1,inu,spsol) - & 887 cg_aers_2bands(RH_num(i,k),inu,spsol)) 888 889 tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int 890 tau_ae_pi(i,k,id,inu) = tmp_var_pi* tau_ae2b_int 891 piz_ae(i,k,id,inu) = piz_ae2b_int 892 cg_ae(i,k,id,inu)= cg_ae2b_int 893 894 ENDDO 895 ENDDO 896 ENDDO 897 898 ENDIF 899 900 ELSE ! For all aerosol insoluble components 901 902 !CDIR ON_ADB(alpha_aers_2bands) 903 !CDIR ON_ADB(piz_aers_2bands) 904 !CDIR ON_ADB(cg_aers_2bands) 723 905 DO k=1, KLEV 724 DO i=1, KLON 725 726 rh=MIN(RHcl(i,k)*100.,RH_MAX) 727 RH_num = INT( rh/10. + 1.) 728 729 IF (rh.GT.85.) RH_num=10 730 IF (rh.GT.90.) RH_num=11 731 DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num)) 732 733 ! DO inu=1,nbands 734 IF (soluble) THEN 735 736 ! First optical parameters are computed for seasalt 737 IF (spss.NE.0) THEN 738 H=rh/100 739 IF (spss.EQ.1) THEN !accumulation mode 740 ! band 1 741 tau_ae2b_int(i,k,1)=A1_ASSSM_b1(k)+A2_ASSSM_b1(k)*H+A3_ASSSM_b1(k)/(H-1.05) 742 piz_ae2b_int(i,k,1)=1-B1_ASSSM_b1(k)-B2_ASSSM_b1(k)*H 743 cg_ae2b_int(i,k,1)=C1_ASSSM_b1(k)+C2_ASSSM_b1(k)*H 744 !band 2 745 tau_ae2b_int(i,k,2)=A1_ASSSM_b2(k)+A2_ASSSM_b2(k)*H+A3_ASSSM_b2(k)/(H-1.05) 746 piz_ae2b_int(i,k,2)=1-B1_ASSSM_b2(k)-B2_ASSSM_b2(k)*H 747 cg_ae2b_int(i,k,2)=C1_ASSSM_b2(k)+C2_ASSSM_b2(k)*H 748 ENDIF 749 IF (spss.EQ.2) THEN !coarse mode 750 ! band 1 751 tau_ae2b_int(i,k,1)=A1_CSSSM_b1(k)+A2_CSSSM_b1(k)*H+A3_CSSSM_b1(k)/(H-1.05) 752 piz_ae2b_int(i,k,1)=1-B1_CSSSM_b1(k)-B2_CSSSM_b1(k)*H 753 cg_ae2b_int(i,k,1)=C1_CSSSM_b1(k)+C2_CSSSM_b1(k)*H 754 ! band 2 755 tau_ae2b_int(i,k,2)=A1_CSSSM_b2(k)+A2_CSSSM_b2(k)*H+A3_CSSSM_b2(k)/(H-1.05) 756 piz_ae2b_int(i,k,2)=1-B1_CSSSM_b2(k)-B2_CSSSM_b2(k)*H 757 cg_ae2b_int(i,k,2)=C1_CSSSM_b2(k)+C2_CSSSM_b2(k)*H 758 ENDIF 759 IF (spss.EQ.3) THEN !super coarse mode 760 ! band 1 761 tau_ae2b_int(i,k,1)=A1_SSSSM_b1(k)+A2_SSSSM_b1(k)*H+A3_SSSSM_b1(k)/(H-1.05) 762 piz_ae2b_int(i,k,1)=1-B1_SSSSM_b1(k)-B2_SSSSM_b1(k)*H 763 cg_ae2b_int(i,k,1)=C1_SSSSM_b1(k)+C2_SSSSM_b1(k)*H 764 ! band 2 765 tau_ae2b_int(i,k,2)=A1_SSSSM_b2(k)+A2_SSSSM_b2(k)*H+A3_SSSSM_b2(k)/(H-1.05) 766 piz_ae2b_int(i,k,2)=1-B1_SSSSM_b2(k)-B2_SSSSM_b2(k)*H 767 cg_ae2b_int(i,k,2)=C1_SSSSM_b2(k)+C2_SSSSM_b2(k)*H 768 ENDIF 769 ELSE 770 tau_ae2b_int(i,k,inu)= & 771 alpha_aers_2bands(RH_num,inu,spsol)+ & 772 DELTA* (alpha_aers_2bands(RH_num+1,inu,spsol) - & 773 alpha_aers_2bands(RH_num,inu,spsol)) 774 775 piz_ae2b_int(i,k,inu)= & 776 piz_aers_2bands(RH_num,inu,spsol) + & 777 DELTA* (piz_aers_2bands(RH_num+1,inu,spsol) - & 778 piz_aers_2bands(RH_num,inu,spsol)) 779 780 cg_ae2b_int(i,k,inu)= & 781 cg_aers_2bands(RH_num,inu,spsol) + & 782 DELTA* (cg_aers_2bands(RH_num+1,inu,spsol) - & 783 cg_aers_2bands(RH_num,inu,spsol)) 784 ENDIF 785 786 tau_ae(i,k,aerosol_name(m),inu) = & 787 mass_temp(i,k,spsol)*1000.*zdp1(i,k)*delt* & 788 tau_ae2b_int(i,k,inu)*fac 789 !RAF fac 790 tau_ae_pi(i,k,aerosol_name(m),inu) = & 791 mass_temp_pi(i,k,spsol)*1000.*zdp1(i,k)*delt* & 792 tau_ae2b_int(i,k,inu)*fac 793 794 ELSE ! For all aerosol insoluble components 795 tau_ae2b_int(i,k,inu) = alpha_aeri_2bands(inu,spinsol) 796 piz_ae2b_int(i,k,inu) = piz_aeri_2bands(inu,spinsol) 797 cg_ae2b_int(i,k,inu) = cg_aeri_2bands(inu,spinsol) 798 799 tau_ae(i,k,aerosol_name(m),inu) = & 800 mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* & 801 delt*tau_ae2b_int(i,k,inu)*fac 802 !RAF fac 803 tau_ae_pi(i,k,aerosol_name(m),inu) = & 804 mass_temp_pi(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)* & 805 delt*tau_ae2b_int(i,k,inu)*fac 806 ENDIF 807 808 piz_ae(i,k,aerosol_name(m),inu) = piz_ae2b_int(i,k,inu) 809 810 cg_ae(i,k,aerosol_name(m),inu)= cg_ae2b_int(i,k,inu) 811 812 ENDDO ! nbands : boucle sur les bandes spectrale 813 ENDDO ! klon : Boucle sur les points geographiques (grille horizontale) 814 ENDDO ! klev : Boucle sur les niveaux verticaux 815 ENDDO ! nb_aer : Boucle sur les masses de traceurs 816 906 DO i=1, KLON 907 tmp_var=mass_temp(i,k,naero_soluble+ spinsol)*1000.*zdp1(i,k)*delt*fac 908 tmp_var_pi=mass_temp_pi(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)*delt*fac 909 !CDIR UNROLL=nbands 910 DO inu=1,nbands 911 tau_ae2b_int = alpha_aeri_2bands(inu,spinsol) 912 piz_ae2b_int = piz_aeri_2bands(inu,spinsol) 913 cg_ae2b_int = cg_aeri_2bands(inu,spinsol) 914 915 tau_ae(i,k,id,inu) = tmp_var*tau_ae2b_int 916 tau_ae_pi(i,k,id,inu) = tmp_var_pi*tau_ae2b_int 917 piz_ae(i,k,id,inu) = piz_ae2b_int 918 cg_ae(i,k,id,inu)= cg_ae2b_int 919 ENDDO 920 ENDDO 921 ENDDO 922 923 ENDIF ! soluble 924 925 ENDDO ! nb_aer 926 927 DO m=1,nb_aer 928 IF (.NOT. used_aer(m)) THEN 929 tau_ae(:,:,:,:)=0. 930 tau_ae_pi(:,:,:,:)=0. 931 piz_ae(:,:,:,:)=0. 932 cg_ae(:,:,:,:)=0. 933 ENDIF 934 ENDDO 817 935 818 936 DO inu=1, nbands 819 DO mrfspecies=1,naero_grp 937 DO mrfspecies=1,naero_grp 938 IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER 820 939 DO k=1, KLEV 821 DO i=1, KLON 822 IF (mrfspecies .EQ. 2) THEN ! = total aerosol AER 823 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ & 824 tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ & 825 tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ & 826 tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+tau_ae(i,k,id_SSSSM,inu)+ & 827 tau_ae(i,k,id_CIDUSTM,inu) 828 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 940 DO i=1, KLON 941 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu)+ & 942 tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu)+ & 943 tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu)+ & 944 tau_ae(i,k,id_ASSSM,inu)+tau_ae(i,k,id_CSSSM,inu)+ & 945 tau_ae(i,k,id_SSSSM,inu)+ tau_ae(i,k,id_CIDUSTM,inu) 946 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 829 947 830 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & 831 tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ & 832 tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ & 833 tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ & 834 tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ & 835 tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ & 836 tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ & 837 tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ & 838 tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ & 839 tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu))/tau_allaer(i,k,mrfspecies,inu) 840 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 841 842 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)+ & 948 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)+ & 949 tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)+ & 950 tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)+ & 951 tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)+ & 952 tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu)+ & 953 tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)+ & 954 tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu)+ & 955 tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)+ & 956 tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)+ & 957 tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)) & 958 /tau_allaer(i,k,mrfspecies,inu) 959 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 960 961 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)+ & 843 962 tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu)*cg_ae(i,k,id_CSSO4M,inu)+ & 844 963 tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu)*cg_ae(i,k,id_ASBCM,inu)+ & … … 851 970 tau_ae(i,k,id_CIDUSTM,inu)*piz_ae(i,k,id_CIDUSTM,inu)*cg_ae(i,k,id_CIDUSTM,inu))/ & 852 971 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 853 854 ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT 972 ENDDO 973 ENDDO 974 975 ELSEIF (mrfspecies .EQ. 3) THEN ! = natural aerosol NAT 976 977 DO k=1, KLEV 978 DO i=1, KLON 855 979 !RAF 856 980 tau_allaer(i,k,mrfspecies,inu)=tau_ae_pi(i,k,id_ASSO4M,inu)+ & … … 892 1016 cg_ae(i,k,id_CIDUSTM,inu))/ & 893 1017 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 894 895 !END RAF 896 ELSEIF (mrfspecies .EQ. 4) THEN ! = BC 897 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu) 898 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 899 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) & 1018 ENDDO 1019 ENDDO 1020 1021 ELSEIF (mrfspecies .EQ. 4) THEN ! = BC 1022 DO k=1, KLEV 1023 DO i=1, KLON 1024 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASBCM,inu)+tau_ae(i,k,id_AIBCM,inu) 1025 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 1026 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASBCM,inu)*piz_ae(i,k,id_ASBCM,inu) & 900 1027 +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu))/ & 901 1028 tau_allaer(i,k,mrfspecies,inu) 902 903 1029 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 1030 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)& 904 1031 +tau_ae(i,k,id_AIBCM,inu)*piz_ae(i,k,id_AIBCM,inu)*cg_ae(i,k,id_AIBCM,inu))/ & 905 1032 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 906 907 ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4 908 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu) 909 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 910 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) & 1033 ENDDO 1034 ENDDO 1035 1036 ELSEIF (mrfspecies .EQ. 5) THEN ! = SO4 1037 1038 DO k=1, KLEV 1039 DO i=1, KLON 1040 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASSO4M,inu)+tau_ae(i,k,id_CSSO4M,inu) 1041 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 1042 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_CSSO4M,inu)*piz_ae(i,k,id_CSSO4M,inu) & 911 1043 +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu))/ & 912 1044 tau_allaer(i,k,mrfspecies,inu) 913 914 1045 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 1046 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)& 915 1047 +tau_ae(i,k,id_ASSO4M,inu)*piz_ae(i,k,id_ASSO4M,inu)*cg_ae(i,k,id_ASSO4M,inu))/ & 916 1048 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 917 918 ELSEIF (mrfspecies .EQ. 6) THEN ! = POM 919 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu) 920 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 921 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) & 1049 ENDDO 1050 ENDDO 1051 1052 ELSEIF (mrfspecies .EQ. 6) THEN ! = POM 1053 1054 DO k=1, KLEV 1055 DO i=1, KLON 1056 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_ASPOMM,inu)+tau_ae(i,k,id_AIPOMM,inu) 1057 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 1058 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASPOMM,inu)*piz_ae(i,k,id_ASPOMM,inu) & 922 1059 +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu))/ & 923 1060 tau_allaer(i,k,mrfspecies,inu) 924 925 1061 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 1062 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)& 926 1063 +tau_ae(i,k,id_AIPOMM,inu)*piz_ae(i,k,id_AIPOMM,inu)*cg_ae(i,k,id_AIPOMM,inu))/ & 927 1064 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 928 929 ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST 930 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu) 931 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 932 piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu) 933 cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu) 934 935 ELSEIF (mrfspecies .EQ. 8) THEN ! = SS 936 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) 937 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 938 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) & 939 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) & 940 +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ & 941 tau_allaer(i,k,mrfspecies,inu) 942 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 943 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)& 944 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) & 945 +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ & 946 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 947 948 ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3 949 tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary 950 piz_allaer(i,k,mrfspecies,inu)=0. 951 cg_allaer(i,k,mrfspecies,inu)=0. 952 ENDIF 953 ENDDO 1065 ENDDO 954 1066 ENDDO 955 ENDDO 1067 1068 ELSEIF (mrfspecies .EQ. 7) THEN ! = DUST 1069 1070 DO k=1, KLEV 1071 DO i=1, KLON 1072 tau_allaer(i,k,mrfspecies,inu)=tau_ae(i,k,id_CIDUSTM,inu) 1073 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 1074 piz_allaer(i,k,mrfspecies,inu)=piz_ae(i,k,id_CIDUSTM,inu) 1075 cg_allaer(i,k,mrfspecies,inu)=cg_ae(i,k,id_CIDUSTM,inu) 1076 ENDDO 1077 ENDDO 1078 1079 ELSEIF (mrfspecies .EQ. 8) THEN ! = SS 1080 1081 DO k=1, KLEV 1082 DO i=1, KLON 1083 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) 1084 tau_allaer(i,k,mrfspecies,inu)=MAX(tau_allaer(i,k,mrfspecies,inu),1e-20) 1085 piz_allaer(i,k,mrfspecies,inu)=(tau_ae(i,k,id_ASSSM,inu)*piz_ae(i,k,id_ASSSM,inu) & 1086 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu) & 1087 +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu))/ & 1088 tau_allaer(i,k,mrfspecies,inu) 1089 piz_allaer(i,k,mrfspecies,inu)=MAX(piz_allaer(i,k,mrfspecies,inu),1e-20) 1090 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)& 1091 +tau_ae(i,k,id_CSSSM,inu)*piz_ae(i,k,id_CSSSM,inu)*cg_ae(i,k,id_CSSSM,inu) & 1092 +tau_ae(i,k,id_SSSSM,inu)*piz_ae(i,k,id_SSSSM,inu)*cg_ae(i,k,id_SSSSM,inu))/ & 1093 (tau_allaer(i,k,mrfspecies,inu)*piz_allaer(i,k,mrfspecies,inu)) 1094 ENDDO 1095 ENDDO 1096 1097 ELSEIF (mrfspecies .EQ. 9) THEN ! = NO3 1098 1099 DO k=1, KLEV 1100 DO i=1, KLON 1101 tau_allaer(i,k,mrfspecies,inu)=0. ! preliminary 1102 piz_allaer(i,k,mrfspecies,inu)=0. 1103 cg_allaer(i,k,mrfspecies,inu)=0. 1104 ENDDO 1105 ENDDO 1106 1107 ELSE 1108 1109 DO k=1, KLEV 1110 DO i=1, KLON 1111 tau_allaer(i,k,mrfspecies,inu)=0. 1112 piz_allaer(i,k,mrfspecies,inu)=0. 1113 cg_allaer(i,k,mrfspecies,inu)=0. 1114 ENDDO 1115 ENDDO 1116 1117 ENDIF 1118 1119 ENDDO 956 1120 ENDDO 957 1121 -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/aeropt_5wv.F90
r1249 r1250 126 126 127 127 REAL,PARAMETER :: RH_tab(nbre_RH)=(/0.,10.,20.,30.,40.,50.,60.,70.,80.,85.,90.,95./) 128 REAL :: DELTA , rh, H129 REAL :: tau_ae5wv_int (KLON,KLEV,las)! Intermediate computation of epaisseur optique aerosol130 REAL :: piz_ae5wv_int (KLON,KLEV,las)! Intermediate single scattering albedo aerosol131 REAL :: cg_ae5wv_int (KLON,KLEV,las)! Intermediate asymmetry parameter aerosol128 REAL :: DELTA(klon,klev), rh(klon,klev), H 129 REAL :: tau_ae5wv_int ! Intermediate computation of epaisseur optique aerosol 130 REAL :: piz_ae5wv_int ! Intermediate single scattering albedo aerosol 131 REAL :: cg_ae5wv_int ! Intermediate asymmetry parameter aerosol 132 132 REAL, PARAMETER :: RH_MAX=95. 133 133 REAL :: taue670(KLON) ! epaisseur optique aerosol absorption 550 nm … … 158 158 ! Proprietes optiques 159 159 ! 160 REAL :: radry = 287.054 ! dry air mass constant 161 160 REAL :: radry = 287.054 161 REAL :: tau_tmp ! dry air mass constant 162 REAL :: fact_RH(nbre_RH) 163 LOGICAL :: used_tau(naero_spc) 164 162 165 DATA presnivs_19/& 163 166 100426.5, 98327.6, 95346.5, 90966.8, 84776.9, & … … 477 480 ! Initialisations 478 481 ai(:) = 0. 479 tau_ae5wv_int(:,:,:) = 0.480 piz_ae5wv_int(:,:,:) = 0.481 cg_ae5wv_int(:,:,:) = 0.482 482 tausum(:,:,:) = 0. 483 tau(:,:,:,:) = 0.484 483 485 484 … … 489 488 ! IF (pplay(i,k).EQ.0) stop 'stop aeropt_5wv p ' 490 489 zrho=pplay(i,k)/t_seri(i,k)/RD ! kg/m3 490 !CDIR UNROLL=naero_spc 491 491 mass_temp(i,k,:) = m_allaer(i,k,:) / zrho / 1.e+9 492 zdp1(i,k)=pdel(i,k)/(gravit*delt) ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] 493 492 494 ENDDO 493 495 ENDDO … … 545 547 ! Calculations that need to be done since we are not in the subroutines INCA 546 548 ! 547 ! air mass auxiliary variable --> zdp1 [kg/(m^2 *s)] 548 zdp1=pdel(:,:)/(gravit*delt) 549 549 550 !!CDIR ON_ADB(RH_tab) 551 !CDIR ON_ADB(fact_RH) 552 !CDIR NOVECTOR 553 DO RH_Num=1,nbre_RH-1 554 fact_RH(RH_num)=1./(RH_tab(RH_num+1)-RH_tab(RH_num)) 555 ENDDO 556 557 DO k=1, KLEV 558 !!CDIR ON_ADB(RH_tab) 559 !CDIR ON_ADB(fact_RH) 560 DO i=1, KLON 561 rh(i,k)=MIN(RHcl(i,k)*100.,RH_MAX) 562 RH_num = INT( rh(i,k)/10. + 1.) 563 IF (rh(i,k).GT.85.) RH_num=10 564 IF (rh(i,k).GT.90.) RH_num=11 565 DELTA(i,k)=(rh(i,k)-RH_tab(RH_num))*fact_RH(RH_num) 566 ENDDO 567 ENDDO 568 569 !CDIR SHORTLOOP 570 used_tau(:)=.FALSE. 571 550 572 DO m=1,nb_aer ! tau is only computed for each mass 551 573 fac=1.0 … … 595 617 CYCLE 596 618 ENDIF 597 619 620 used_tau(spsol)=.TRUE. 598 621 DO la=1,las 599 tau3d(:,:)=0. 600 piz3d(:,:)=0. 601 cg3d(:,:)=0. 602 abs3d(:,:)=0. 603 604 DO k=1, KLEV 605 DO i=1, KLON 606 607 rh=MIN(RHcl(i,k)*100.,RH_MAX) 608 RH_num = INT( rh/10. + 1.) 609 610 IF (rh.GT.85.) RH_num=10 611 IF (rh.GT.90.) RH_num=11 612 DELTA=(rh-RH_tab(RH_num))/(RH_tab(RH_num+1)-RH_tab(RH_num)) 613 614 IF (soluble) THEN 615 616 IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm 617 H=rh/100 618 IF (spss.EQ.1) THEN !accumulation mode 619 tau_ae5wv_int(i,k,la)=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05) 620 piz_ae5wv_int(i,k,la)=1-B1_ASSSM(k)-B2_ASSSM(k)*H 621 cg_ae5wv_int(i,k,la)=C1_ASSSM(k)+C2_ASSSM(k)*H 622 ENDIF 623 IF (spss.EQ.2) THEN !coarse mode 624 tau_ae5wv_int(i,k,la)=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05) 625 piz_ae5wv_int(i,k,la)=1-B1_CSSSM(k)-B2_CSSSM(k)*H 626 cg_ae5wv_int(i,k,la)=C1_CSSSM(k)+C2_CSSSM(k)*H 627 ENDIF 628 IF (spss.EQ.3) THEN !super coarse mode 629 tau_ae5wv_int(i,k,la)=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05) 630 piz_ae5wv_int(i,k,la)=1-B1_SSSSM(k)-B2_SSSSM(k)*H 631 cg_ae5wv_int(i,k,la)=C1_SSSSM(k)+C2_SSSSM(k)*H 632 ENDIF 633 ELSE 634 635 tau_ae5wv_int(i,k,la) = & 636 alpha_aers_5wv(RH_num,la,spsol)+DELTA* & 637 (alpha_aers_5wv(RH_num+1,la,spsol) - & 638 alpha_aers_5wv(RH_num,la,spsol)) 639 640 piz_ae5wv_int(i,k,la) = & 641 piz_aers_5wv(RH_num,la,spsol)+DELTA* & 642 (piz_aers_5wv(RH_num+1,la,spsol) - & 643 piz_aers_5wv(RH_num,la,spsol)) 644 645 cg_ae5wv_int(i,k,la) = & 646 cg_aers_5wv(RH_num,la,spsol)+DELTA* & 647 (cg_aers_5wv(RH_num+1,la,spsol) - & 648 cg_aers_5wv(RH_num,la,spsol)) 649 ENDIF 650 651 tau3d(i,k) = & 652 mass_temp(i,k,spsol)*1000.*zdp1(i,k)*tau_ae5wv_int(i,k,la)*delt*fac 653 654 ELSE ! For insoluble aerosol 655 tau_ae5wv_int(i,k,la) = alpha_aeri_5wv(la,spinsol) 656 piz_ae5wv_int(i,k,la) = piz_aeri_5wv(la,spinsol) 657 cg_ae5wv_int(i,k,la) = cg_aeri_5wv(la,spinsol) 658 659 tau3d(i,k) = & 660 mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* & 661 tau_ae5wv_int(i,k,la)*delt*fac 622 623 IF (soluble) THEN 624 625 IF((la.EQ.2).AND.(spss.NE.0)) THEN !la=2 corresponds to 550 nm 626 IF (spss.EQ.1) THEN !accumulation mode 627 DO k=1, KLEV 628 !CDIR ON_ADB(A1_ASSSM) 629 !CDIR ON_ADB(A2_ASSSM) 630 !CDIR ON_ADB(A3_ASSSM) 631 DO i=1, KLON 632 H=rh(i,k)/100 633 tau_ae5wv_int=A1_ASSSM(k)+A2_ASSSM(k)*H+A3_ASSSM(k)/(H-1.05) 634 tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & 635 *tau_ae5wv_int*delt*fac 636 tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) 637 ENDDO 638 ENDDO 662 639 ENDIF 663 664 ENDDO ! Boucle sur les points géographiques (grille horizontale) 665 ENDDO ! Boucle sur les niveaux verticaux 666 667 IF (soluble) THEN 668 669 tau(:,:,la,spsol)=tau3d(:,:) 670 640 641 IF (spss.EQ.2) THEN !coarse mode 642 DO k=1, KLEV 643 !CDIR ON_ADB(A1_CSSSM) 644 !CDIR ON_ADB(A2_CSSSM) 645 !CDIR ON_ADB(A3_CSSSM) 646 DO i=1, KLON 647 H=rh(i,k)/100 648 tau_ae5wv_int=A1_CSSSM(k)+A2_CSSSM(k)*H+A3_CSSSM(k)/(H-1.05) 649 tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & 650 *tau_ae5wv_int*delt*fac 651 tausum(i,la,spsol) = tausum(i,la,spsol)+tau(i,k,la,spsol) 652 ENDDO 653 ENDDO 654 ENDIF 655 656 IF (spss.EQ.3) THEN !super coarse mode 657 DO k=1, KLEV 658 !CDIR ON_ADB(A1_SSSSM) 659 !CDIR ON_ADB(A2_SSSSM) 660 !CDIR ON_ADB(A3_SSSSM) 661 DO i=1, KLON 662 H=rh(i,k)/100 663 tau_ae5wv_int=A1_SSSSM(k)+A2_SSSSM(k)*H+A3_SSSSM(k)/(H-1.05) 664 tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & 665 *tau_ae5wv_int*delt*fac 666 tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) 667 ENDDO 668 ENDDO 669 ENDIF 670 671 ELSE 671 672 DO k=1, KLEV 672 DO i=1,KLON 673 tausum(i,la,spsol)=tausum(i,la,spsol)+tau3d(i,k) 673 !CDIR ON_ADB(alpha_aers_5wv) 674 DO i=1, KLON 675 tau_ae5wv_int = alpha_aers_5wv(RH_num,la,spsol)+DELTA(i,k)* & 676 (alpha_aers_5wv(RH_num+1,la,spsol) - & 677 alpha_aers_5wv(RH_num,la,spsol)) 678 679 tau(i,k,la,spsol) = mass_temp(i,k,spsol)*1000.*zdp1(i,k) & 680 *tau_ae5wv_int*delt*fac 681 tausum(i,la,spsol)=tausum(i,la,spsol)+tau(i,k,la,spsol) 674 682 ENDDO 675 683 ENDDO 676 ELSE ! For insoluble aerosol 677 tau(:,:,la,naero_soluble+spinsol)=tau3d(:,:) 678 679 DO k=1, KLEV 680 DO i=1,KLON 681 tausum(i,la,naero_soluble+spinsol)= & 682 tausum(i,la,naero_soluble+spinsol)+tau3d(i,k) 683 ENDDO ! Boucle sur les points géographiques (grille horizontale) 684 ENDDO ! Boucle sur les niveaux verticaux 684 ENDIF 685 686 ELSE ! For insoluble aerosol 687 DO k=1, KLEV 688 !CDIR ON_ADB(alpha_aeri_5wv) 689 DO i=1, KLON 690 tau_ae5wv_int = alpha_aeri_5wv(la,spinsol) 691 tau(i,k,la,naero_soluble+spinsol) = mass_temp(i,k,naero_soluble+spinsol)*1000.*zdp1(i,k)* & 692 tau_ae5wv_int*delt*fac 693 tausum(i,la,naero_soluble+spinsol)= tausum(i,la,naero_soluble+spinsol) & 694 +tau(i,k,la,naero_soluble+spinsol) 695 ENDDO 696 ENDDO 685 697 ENDIF 686 687 698 ENDDO ! boucle sur les longueurs d'onde 688 699 ENDDO ! Boucle sur les masses de traceurs 700 701 DO m=1,naero_spc 702 IF (.NOT.used_tau(m)) tau(:,:,:,m)=0. 703 ENDDO 689 704 ! 690 705 ! -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/carbon_cycle_mod.F90
r1249 r1250 9 9 10 10 ! Variables read from parmeter file physiq.def 11 LOGICAL, SAVE,PUBLIC :: carbon_cycle_tr ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys11 LOGICAL, PUBLIC :: carbon_cycle_tr ! 3D transport of CO2 in the atmosphere, parameter read in conf_phys 12 12 !$OMP THREADPRIVATE(carbon_cycle_tr) 13 LOGICAL, SAVE,PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES)13 LOGICAL, PUBLIC :: carbon_cycle_cpl ! Coupling of CO2 fluxes between LMDZ/ORCHIDEE and LMDZ/OCEAN(PISCES) 14 14 !$OMP THREADPRIVATE(carbon_cycle_cpl) 15 15 LOGICAL :: carbon_cycle_emis_comp=.FALSE. ! Calculation of emission compatible … … 17 17 ! Scalare values when no transport, from physiq.def 18 18 REAL :: fos_fuel_s ! carbon_cycle_fos_fuel dans physiq.def 19 !$OMP THREADPRIVATE(fos_fuel_s) 19 20 REAL :: emis_land_s ! not yet implemented 21 !$OMP THREADPRIVATE(emis_land_s) 20 22 21 23 INTEGER :: ntr_co2 ! Number of tracers concerning the carbon cycle … … 25 27 INTEGER :: id_fco2_land_use ! - " - 26 28 INTEGER :: id_fco2_fos_fuel ! - " - 27 29 !$OMP THREADPRIVATE(ntr_co2, id_fco2_tot, id_fco2_ocn, id_fco2_land, id_fco2_land_use, id_fco2_fos_fuel) 30 28 31 REAL, DIMENSION(:), ALLOCATABLE :: fos_fuel ! CO2 fossil fuel emission from file [gC/m2/d] 32 !$OMP THREADPRIVATE(fos_fuel) 29 33 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_ocn_day ! flux CO2 from ocean for 1 day (cumulated) [gC/m2/d] 34 !$OMP THREADPRIVATE(fco2_ocn_day) 30 35 REAL, DIMENSION(:), ALLOCATABLE :: fco2_land_day ! flux CO2 from land for 1 day (cumulated) [gC/m2/d] 36 !$OMP THREADPRIVATE(fco2_land_day) 31 37 REAL, DIMENSION(:), ALLOCATABLE :: fco2_lu_day ! Emission from land use change for 1 day (cumulated) [gC/m2/d] 38 !$OMP THREADPRIVATE(fco2_lu_day) 32 39 33 40 ! Following 2 fields will be initialized in surf_land_orchidee at each time step 34 41 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_land_inst ! flux CO2 from land at one time step 42 !$OMP THREADPRIVATE(fco2_land_inst) 35 43 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: fco2_lu_inst ! Emission from land use change at one time step 44 !$OMP THREADPRIVATE(fco2_lu_inst) 36 45 37 46 ! Calculated co2 field to be send to the ocean via the coupler and to ORCHIDEE 38 47 REAL, DIMENSION(:), ALLOCATABLE, PUBLIC :: co2_send 39 48 !$OMP THREADPRIVATE(co2_send) 40 49 41 50 CONTAINS -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/clouds_gno.F
r1162 r1250 47 47 48 48 INTEGER i,K, n, m 49 REAL mu(klon), qsat (klon), delta(klon), beta(klon)50 real zu2 (klon),zv2(klon)51 REAL xx(klon), aux(klon), coeff (klon), block(klon)52 REAL dist (klon), fprime(klon), det(klon)53 REAL pi, u (klon), v(klon), erfcu(klon), erfcv(klon)54 REAL xx1 (klon), xx2(klon)55 real erf, kkk49 REAL mu(klon), qsat, delta(klon), beta(klon) 50 real zu2,zv2 51 REAL xx(klon), aux(klon), coeff, block 52 REAL dist, fprime, det 53 REAL pi, u, v, erfcu, erfcv 54 REAL xx1, xx2 55 real erf,hsqrtlog_2,v2 56 56 real sqrtpi,sqrt2,zx1,zx2,exdel 57 57 c lconv = true si le calcul a converge (entre autre si qsub < min_q) 58 58 LOGICAL lconv(klon) 59 59 60 cym 61 cldf(:,:)=0.0 60 !cdir arraycomb 61 cldf (1:klon,1:ND)=0.0 ! cym 62 ratqsc(1:klon,1:ND)=0.0 63 ptconv(1:klon,1:ND)=.false. 64 !cdir end arraycomb 62 65 63 66 pi = ACOS(-1.) 64 67 sqrtpi=sqrt(pi) 65 68 sqrt2=sqrt(2.) 66 67 ptconv=.false. 68 ratqsc=0. 69 69 hsqrtlog_2=0.5*SQRT(log(2.)) 70 70 71 71 DO 500 K = 1, ND … … 74 74 mu(i) = R(i,K) 75 75 mu(i) = MAX(mu(i),min_mu) 76 qsat (i)= RS(i,K)77 qsat (i) = MAX(qsat(i),min_mu)78 delta(i) = log(mu(i)/qsat (i))79 76 qsat = RS(i,K) 77 qsat = MAX(qsat,min_mu) 78 delta(i) = log(mu(i)/qsat) 79 c enddo ! vector 80 80 81 81 C … … 106 106 c suffisamment d'eau nuageuse. 107 107 108 108 c do i=1,klon ! vector 109 109 110 110 IF ( QSUB(i,K) .lt. min_Q ) THEN … … 124 124 c -- roots of equation v > vmax: 125 125 126 det (i) = delta(i) + vmax(i)**2.127 if (det (i).LE.0.0) vmax(i) = vmax0 + 1.0128 det (i) = delta(i) + vmax(i)**2.129 130 if (det (i).LE.0.) then126 det = delta(i) + vmax(i)*vmax(i) 127 if (det.LE.0.0) vmax(i) = vmax0 + 1.0 128 det = delta(i) + vmax(i)*vmax(i) 129 130 if (det.LE.0.) then 131 131 xx(i) = -0.0001 132 132 else 133 133 zx1=-sqrt2*vmax(i) 134 zx2=SQRT(1.0+delta(i)/(vmax(i)* *2.))135 xx1 (i)=zx1*(1.0-zx2)136 xx2 (i)=zx1*(1.0+zx2)137 xx(i) = 1.01 * xx1 (i)138 if ( xx1 (i) .GE. 0.0 ) xx(i) = 0.5*xx2(i)134 zx2=SQRT(1.0+delta(i)/(vmax(i)*vmax(i))) 135 xx1=zx1*(1.0-zx2) 136 xx2=zx1*(1.0+zx2) 137 xx(i) = 1.01 * xx1 138 if ( xx1 .GE. 0.0 ) xx(i) = 0.5*xx2 139 139 endif 140 if (delta(i).LT.0.) xx(i) = - 0.5*SQRT(log(2.))140 if (delta(i).LT.0.) xx(i) = -hsqrtlog_2 141 141 142 142 ENDIF … … 153 153 if (.not.lconv(i)) then 154 154 155 u(i) = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2) 156 v(i) = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2) 157 158 IF ( v(i) .GT. vmax(i) ) THEN 159 160 IF ( ABS(u(i)) .GT. vmax(i) 155 u = delta(i)/(xx(i)*sqrt2) + xx(i)/(2.*sqrt2) 156 v = delta(i)/(xx(i)*sqrt2) - xx(i)/(2.*sqrt2) 157 v2 = v*v 158 159 IF ( v .GT. vmax(i) ) THEN 160 161 IF ( ABS(u) .GT. vmax(i) 161 162 : .AND. delta(i) .LT. 0. ) THEN 162 163 … … 171 172 endif 172 173 xx(i) = -SQRT(aux(i)) 173 block (i) = EXP(-v(i)*v(i)) / v(i)/ sqrtpi174 dist (i)= 0.0175 fprime (i)= 1.0174 block = EXP(-v*v) / v / sqrtpi 175 dist = 0.0 176 fprime = 1.0 176 177 177 178 ELSE … … 179 180 c -- erfv -> 1.0, use an asymptotic expression of erfv for v large: 180 181 181 erfcu (i) = 1.0-ERF(u(i))182 erfcu = 1.0-ERF(u) 182 183 c !!! ATTENTION : rajout d'un seuil pour l'exponentiel 183 aux(i) = sqrtpi*erfcu (i)*EXP(min(v(i)*v(i),100.))184 coeff (i) = 1.0 - 1./2./(v(i)**2.) + 3./4./(v(i)**4.)185 block (i) = coeff(i) * EXP(-v(i)*v(i)) / v(i)/ sqrtpi186 dist (i) = v(i) * aux(i) / coeff(i)- beta(i)187 fprime (i) = 2.0 / xx(i) * (v(i)**2.)188 : * ( coeff(i)*EXP(-delta(i)) - u(i) * aux(i))189 : / coeff (i) / coeff(i)184 aux(i) = sqrtpi*erfcu*EXP(min(v2,100.)) 185 coeff = 1.0 - 0.5/(v2) + 0.75/(v2*v2) 186 block = coeff * EXP(-v2) / v / sqrtpi 187 dist = v * aux(i) / coeff - beta(i) 188 fprime = 2.0 / xx(i) * (v2) 189 : * ( EXP(-delta(i)) - u * aux(i) / coeff ) 190 : / coeff 190 191 191 192 ENDIF ! ABS(u) … … 195 196 c -- general case: 196 197 197 erfcu (i) = 1.0-ERF(u(i))198 erfcv (i) = 1.0-ERF(v(i))199 block (i) = erfcv(i)200 dist (i) = erfcu(i) / erfcv(i)- beta(i)201 zu2 (i)=u(i)*u(i)202 zv2 (i)=v(i)*v(i)203 if(zu2 (i).gt.20..or. zv2(i).gt.20.) then198 erfcu = 1.0-ERF(u) 199 erfcv = 1.0-ERF(v) 200 block = erfcv 201 dist = erfcu / erfcv - beta(i) 202 zu2=u*u 203 zv2=v2 204 if(zu2.gt.20..or. zv2.gt.20.) then 204 205 c print*,'ATTENTION !!! xx(',i,') =', xx(i) 205 206 c print*,'ATTENTION !!! klon,ND,R,RS,QSUB,PTCONV,RATQSC,CLDF', … … 207 208 c .CLDF(i,k) 208 209 c print*,'ATTENTION !!! zu2 zv2 =',zu2(i),zv2(i) 209 zu2 (i)=20.210 zv2 (i)=20.211 fprime (i)= 0.210 zu2=20. 211 zv2=20. 212 fprime = 0. 212 213 else 213 fprime (i) = 2. /sqrtpi /xx(i) /erfcv(i)**2.214 : * ( erfcv (i)*v(i)*EXP(-zu2(i))215 : - erfcu (i)*u(i)*EXP(-zv2(i)) )214 fprime = 2. /sqrtpi /xx(i) /(erfcv*erfcv) 215 : * ( erfcv*v*EXP(-zu2) 216 : - erfcu*u*EXP(-zv2) ) 216 217 endif 217 218 ENDIF ! x … … 223 224 ! stop 224 225 ! endif 225 if (abs(fprime (i)).lt.1.e-11) then226 if (abs(fprime).lt.1.e-11) then 226 227 ! print*,'avant test fprime<.e-11 ' 227 228 ! s ,i,k,lconv(i),u(i),v(i),beta(i),fprime(i) 228 229 ! print*,'klon,ND,R,RS,QSUB', 229 230 ! s klon,ND,R(i,k),rs(i,k),qsub(i,k) 230 fprime (i)=sign(1.e-11,fprime(i))231 fprime=sign(1.e-11,fprime) 231 232 endif 232 233 233 234 234 if ( ABS(dist (i)/beta(i)) .LT. epsilon ) then235 if ( ABS(dist/beta(i)) .LT. epsilon ) then 235 236 c print*,'v-u **2',(v(i)-u(i))**2 236 237 c print*,'exp v-u **2',exp((v(i)-u(i))**2) … … 238 239 lconv(i)=.true. 239 240 c borne pour l'exponentielle 240 ratqsc(i,k)=min(2.*(v (i)-u(i))**2,20.)241 ratqsc(i,k)=min(2.*(v-u)*(v-u),20.) 241 242 ratqsc(i,k)=sqrt(exp(ratqsc(i,k))-1.) 242 CLDF(i,K) = 0.5 * block (i)243 CLDF(i,K) = 0.5 * block 243 244 else 244 xx(i) = xx(i) - dist (i)/fprime(i)245 xx(i) = xx(i) - dist/fprime 245 246 endif 246 247 c print*,'apres test ',i,k,lconv(i) -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/cvltr.F90
r1191 r1250 28 28 29 29 ! Variables locales 30 30 ! REAL,DIMENSION(klon,klev) :: zed 31 31 REAL,DIMENSION(klon,klev,klev) :: zmd 32 32 REAL,DIMENSION(klon,klev,klev) :: za … … 34 34 REAL,DIMENSION(klon,klev) :: zmfp,zmfu 35 35 INTEGER :: i,k,j 36 REAL :: pdtimeRG 36 37 37 38 ! ========================================= 38 39 ! calcul des tendances liees au downdraft 39 40 ! ========================================= 40 zed(:,:)=0. 41 zmfd(:,:)=0. 42 zmfa(:,:)=0. 43 zmfu(:,:)=0. 44 zmfp(:,:)=0. 45 zmd(:,:,:)=0. 46 za(:,:,:)=0. 41 !cdir collapse 42 DO j=1,klev 43 DO i=1,klon 44 ! zed(i,j)=0. 45 zmfd(i,j)=0. 46 zmfa(i,j)=0. 47 zmfu(i,j)=0. 48 zmfp(i,j)=0. 49 END DO 50 END DO 51 !cdir collapse 52 DO k=1,klev 53 DO j=1,klev 54 DO i=1,klon 55 zmd(i,j,k)=0. 56 za (i,j,k)=0. 57 END DO 58 END DO 59 END DO 47 60 ! entrainement 48 49 50 51 52 61 ! DO k=1,klev-1 62 ! DO i=1,klon 63 ! zed(i,k)=max(0.,mp(i,k)-mp(i,k+1)) 64 ! END DO 65 ! END DO 53 66 54 67 ! calcul de la matrice d echange 55 68 ! matrice de distribution de la masse entrainee en k 56 69 57 DO k=1,klev 70 DO k=1,klev-1 58 71 DO i=1,klon 59 zmd(i,k,k)= zed(i,k)72 zmd(i,k,k)=max(0.,mp(i,k)-mp(i,k+1)) 60 73 END DO 61 74 END DO … … 126 139 DO k=1, klev 127 140 DO i=1, klon 141 dx(i,k)=paprs(i,k)-paprs(i,k+1) 142 ENDDO 143 ENDDO 144 pdtimeRG=pdtime*RG 145 !cdir collapse 146 DO k=1, klev 147 DO i=1, klon 128 148 dx(i,k)=(zmfd(i,k)+zmfu(i,k) & 129 +zmfa(i,k)+zmfp(i,k))*pdtime & 130 *RG/(paprs(i,k)-paprs(i,k+1)) 149 +zmfa(i,k)+zmfp(i,k))*pdtimeRG/dx(i,k) 131 150 ! print*,'dx',k,dx(i,k) 132 151 ENDDO -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/fisrtilp.F
r1091 r1250 54 54 REAL frac_impa(klon,klev) 55 55 REAL frac_nucl(klon,klev) 56 real zct (klon),zcl(klon)56 real zct ,zcl 57 57 cAA 58 58 c … … 87 87 REAL ztglace, zt(klon) 88 88 INTEGER nexpo ! exponentiel pour glace/eau 89 REAL zdz(klon),zrho(klon),ztot (klon), zrhol(klon)90 REAL zchau (klon),zfroi(klon),zfice(klon),zneb(klon)89 REAL zdz(klon),zrho(klon),ztot , zrhol(klon) 90 REAL zchau ,zfroi ,zfice(klon),zneb(klon) 91 91 c 92 92 LOGICAL appel1er … … 150 150 cAA Initialisation a 1 des coefs des fractions lessivees 151 151 c 152 !cdir collapse 152 153 DO k = 1, klev 153 154 DO i = 1, klon … … 161 162 c 162 163 cMAf Initialisation a 0 de zoliq 163 164 165 164 c DO i = 1, klon 165 c zoliq(i)=0. 166 c ENDDO 166 167 c Determiner les nuages froids par leur temperature 167 168 c nexpo regle la raideur de la transition eau liquide / eau glace. … … 173 174 c Initialiser les sorties: 174 175 c 176 !cdir collapse 175 177 DO k = 1, klev+1 176 178 DO i = 1, klon … … 180 182 ENDDO 181 183 184 !cdir collapse 182 185 DO k = 1, klev 183 186 DO i = 1, klon … … 194 197 rain(i) = 0.0 195 198 snow(i) = 0.0 196 ENDDO 199 zoliq(i)=0. 200 c ENDDO 197 201 c 198 202 c Initialiser le flux de precipitation a zero 199 203 c 200 204 c DO i = 1, klon 201 205 zrfl(i) = 0.0 202 206 zneb(i) = seuil_neb … … 441 445 zrhol(i) = zrho(i) * zoliq(i) / zneb(i) 442 446 443 if (ptconv(i,k)) then 444 zcl(i)=cld_lc_con 445 zct(i)=1./cld_tau_con 446 else 447 zcl(i)=cld_lc_lsc 448 zct(i)=1./cld_tau_lsc 449 endif 450 c quantit�d'eau ��minier. 451 zchau(i) = zct(i)*dtime/FLOAT(ninter) * zoliq(i) 452 . *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl(i))**2)) *(1.-zfice(i)) 453 c meme chose pour la glace. 454 if (ptconv(i,k)) then 455 zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) 447 IF (zneb(i).EQ.seuil_neb) THEN 448 ztot = 0.0 449 ELSE 450 c quantite d'eau a eliminer: zchau 451 c meme chose pour la glace: zfroi 452 if (ptconv(i,k)) then 453 zcl =cld_lc_con 454 zct =1./cld_tau_con 455 zfroi = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) 456 456 . *fallvc(zrhol(i)) * zfice(i) 457 else 458 zfroi(i) = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) 457 else 458 zcl =cld_lc_lsc 459 zct =1./cld_tau_lsc 460 zfroi = dtime/FLOAT(ninter)/zdz(i)*zoliq(i) 459 461 . *fallvs(zrhol(i)) * zfice(i) 460 endif 461 ztot(i) = zchau(i) + zfroi(i) 462 IF (zneb(i).EQ.seuil_neb) ztot(i) = 0.0 463 ztot(i) = MIN(MAX(ztot(i),0.0),zoliq(i)) 464 zoliq(i) = MAX(zoliq(i)-ztot(i), 0.0) 462 endif 463 zchau = zct *dtime/FLOAT(ninter) * zoliq(i) 464 . *(1.0-EXP(-(zoliq(i)/zneb(i)/zcl )**2)) *(1.-zfice(i)) 465 ztot = zchau + zfroi 466 ztot = MAX(ztot ,0.0) 467 ENDIF 468 ztot = MIN(ztot,zoliq(i)) 469 zoliq(i) = MAX(zoliq(i)-ztot , 0.0) 465 470 radliq(i,k) = radliq(i,k) + zoliq(i)/FLOAT(ninter+1) 466 471 ENDIF -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/hbtm.F
r776 r1250 762 762 endif 763 763 c 764 qsatbef(i) = qqsat ! bug dans la version orig ??? 764 765 endif 765 qsatbef(i) = qqsat766 766 camn ???? cette ligne a deja ete faite normalement ? 767 767 endif -
LMDZ4/branches/LMDZ4-dev/libf/phylmd/radiation_AR4.F
r1231 r1250 155 155 c$OMP THREADPRIVATE(initialized) 156 156 cjq-end 157 INTEGER JK 158 REAL tmp_ 157 159 if(.not.initialized) then 158 160 flag_aer=0. … … 162 164 allocate(ZFSUPAI(KDLON,KFLEV+1)) 163 165 allocate(ZFSDNAI(KDLON,KFLEV+1)) 164 ZFSUPAD(:,:)=0. 165 ZFSDNAD(:,:)=0. 166 ZFSUPAI(:,:)=0. 167 ZFSDNAI(:,:)=0. 168 166 DO JK = 1 , KDLON*(KFLEV+1) 167 ZFSUPAD(JK,1) = 0.0 ! ZFSUPAD(:,:)=0. 168 ZFSDNAD(JK,1) = 0.0 ! ZFSDNAD(:,:)=0. 169 ZFSUPAI(JK,1) = 0.0 ! ZFSUPAI(:,:)=0. 170 ZFSDNAI(JK,1) = 0.0 ! ZFSDNAI(:,:)=0. 171 END DO 169 172 endif 170 173 !rv … … 179 182 IF (MOD(itapsw,swpas).EQ.0) THEN 180 183 c 184 tmp_ = 1./( dobson_u * 1e3 * RG) 185 !cdir collapse 181 186 DO JK = 1 , KFLEV 182 DO JL = 1, KDLON183 ZCLDSW0(JL,JK) = 0.0184 ZOZ(JL,JK) = POZON(JL,JK) / dobson_u / 1e3 / RG *PDP(JL,JK)185 ENDDO187 DO JL = 1, KDLON 188 ZCLDSW0(JL,JK) = 0.0 189 ZOZ(JL,JK) = POZON(JL,JK)*tmp_*PDP(JL,JK) 190 ENDDO 186 191 ENDDO 187 192 C … … 342 347 S PRMU,PSEC,PUD) 343 348 USE dimphy 349 USE radiation_AR4_param, only : 350 S ZPDH2O,ZPDUMG,ZPRH2O,ZPRUMG,RTDH2O,RTDUMG,RTH2O,RTUMG 344 351 IMPLICIT none 345 352 cym#include "dimensions.h" … … 389 396 INTEGER jl, jk, jkp1, jkl, jklp1, ja 390 397 C 391 C* Prescribed Data:392 c393 REAL(KIND=8) ZPDH2O,ZPDUMG394 SAVE ZPDH2O,ZPDUMG395 c$OMP THREADPRIVATE(ZPDH2O,ZPDUMG)396 REAL(KIND=8) ZPRH2O,ZPRUMG397 SAVE ZPRH2O,ZPRUMG398 c$OMP THREADPRIVATE(ZPRH2O,ZPRUMG)399 REAL(KIND=8) RTDH2O,RTDUMG400 SAVE RTDH2O,RTDUMG401 c$OMP THREADPRIVATE(RTDH2O,RTDUMG)402 REAL(KIND=8) RTH2O ,RTUMG403 SAVE RTH2O ,RTUMG404 c$OMP THREADPRIVATE(RTH2O ,RTUMG)405 DATA ZPDH2O,ZPDUMG / 0.8 , 0.75 /406 DATA ZPRH2O,ZPRUMG / 30000., 30000. /407 DATA RTDH2O,RTDUMG / 0.40 , 0.375 /408 DATA RTH2O ,RTUMG / 240. , 240. /409 398 C ------------------------------------------------------------------ 410 399 C … … 537 526 S , PFD , PFU) 538 527 USE dimphy 528 USE radiation_AR4_param, only : RSUN, RRAY 539 529 IMPLICIT none 540 530 cym#include "dimensions.h" … … 623 613 C 624 614 INTEGER jl, jk, k, jaj, ikm1, ikl 625 c 626 c Prescribed Data: 627 c 628 REAL(KIND=8) RSUN(2) 629 SAVE RSUN 630 c$OMP THREADPRIVATE(RSUN) 631 REAL(KIND=8) RRAY(2,6) 632 SAVE RRAY 633 c$OMP THREADPRIVATE(RRAY) 634 DATA RSUN(1) / 0.441676 / 635 DATA RSUN(2) / 0.558324 / 636 DATA (RRAY(1,K),K=1,6) / 637 S .428937E-01, .890743E+00,-.288555E+01, 638 S .522744E+01,-.469173E+01, .161645E+01/ 639 DATA (RRAY(2,K),K=1,6) / 640 S .697200E-02, .173297E-01,-.850903E-01, 641 S .248261E+00,-.302031E+00, .129662E+00/ 615 642 616 C ------------------------------------------------------------------ 643 617 C … … 780 754 S , PFDOWN,PFUP ) 781 755 USE dimphy 756 USE radiation_AR4_param, only : RSUN, RRAY 782 757 IMPLICIT none 783 758 cym#include "dimensions.h" … … 900 875 REAL(KIND=8) ZRMUM1, ZWH2O, ZCNEB, ZAA, ZBB, ZRKI, ZRE11 901 876 C 902 C* Prescribed Data: 903 C 904 REAL(KIND=8) RSUN(2) 905 SAVE RSUN 906 c$OMP THREADPRIVATE(RSUN) 907 REAL(KIND=8) RRAY(2,6) 908 SAVE RRAY 909 c$OMP THREADPRIVATE(RRAY) 910 DATA RSUN(1) / 0.441676 / 911 DATA RSUN(2) / 0.558324 / 912 DATA (RRAY(1,K),K=1,6) / 913 S .428937E-01, .890743E+00,-.288555E+01, 914 S .522744E+01,-.469173E+01, .161645E+01/ 915 DATA (RRAY(2,K),K=1,6) / 916 S .697200E-02, .173297E-01,-.850903E-01, 917 S .248261E+00,-.302031E+00, .129662E+00/ 877 918 878 C 919 879 C ------------------------------------------------------------------ … … 1329 1289 S , PRK , PRMU0 , PTAUAZ, PTRA1 , PTRA2 ) 1330 1290 USE dimphy 1291 USE radiation_AR4_param, only : TAUA, RPIZA, RCGA 1331 1292 IMPLICIT none 1332 1293 cym#include "dimensions.h" … … 1398 1359 REAL(KIND=8) ZBMU0, ZBMU1, ZRE11 1399 1360 C 1400 C* Prescribed Data for Aerosols: 1401 C 1402 REAL(KIND=8) TAUA(2,5), RPIZA(2,5), RCGA(2,5) 1403 SAVE TAUA, RPIZA, RCGA 1404 c$OMP THREADPRIVATE(TAUA, RPIZA, RCGA) 1405 DATA ((TAUA(IN,JA),JA=1,5),IN=1,2) / 1406 S .730719, .912819, .725059, .745405, .682188 , 1407 S .730719, .912819, .725059, .745405, .682188 / 1408 DATA ((RPIZA(IN,JA),JA=1,5),IN=1,2) / 1409 S .872212, .982545, .623143, .944887, .997975 , 1410 S .872212, .982545, .623143, .944887, .997975 / 1411 DATA ((RCGA (IN,JA),JA=1,5),IN=1,2) / 1412 S .647596, .739002, .580845, .662657, .624246 , 1413 S .647596, .739002, .580845, .662657, .624246 / 1361 1414 1362 C ------------------------------------------------------------------ 1415 1363 C … … 1419 1367 100 CONTINUE 1420 1368 C 1369 !cdir collapse 1421 1370 DO 103 JK = 1 , KFLEV+1 1422 1371 DO 102 JA = 1 , 6 … … 2179 2128 SUBROUTINE SWTT_LMDAR4 (KNU,KA,PU,PTR) 2180 2129 USE dimphy 2130 USE radiation_AR4_param, only : APAD, BPAD, D 2181 2131 IMPLICIT none 2182 2132 cym#include "dimensions.h" … … 2225 2175 INTEGER jl, i,j 2226 2176 C 2227 C* Prescribed Data: 2228 C 2229 REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3) 2230 SAVE APAD, BPAD, D 2231 c$OMP THREADPRIVATE(APAD, BPAD, D) 2232 DATA ((APAD(1,I,J),I=1,3),J=1,7) / 2233 S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, 2234 S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01, 2235 S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00, 2236 S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02, 2237 S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02, 2238 S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02, 2239 S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 / 2240 DATA ((APAD(2,I,J),I=1,3),J=1,7) / 2241 S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, 2242 S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02, 2243 S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00, 2244 S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00, 2245 S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00, 2246 S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00, 2247 S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 / 2248 C 2249 DATA ((BPAD(1,I,J),I=1,3),J=1,7) / 2250 S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, 2251 S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01, 2252 S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00, 2253 S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02, 2254 S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02, 2255 S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02, 2256 S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 / 2257 DATA ((BPAD(2,I,J),I=1,3),J=1,7) / 2258 S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, 2259 S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02, 2260 S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01, 2261 S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00, 2262 S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00, 2263 S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00, 2264 S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 / 2265 c 2266 DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 / 2267 DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 / 2177 2268 2178 C 2269 2179 C----------------------------------------------------------------------- … … 2297 2207 SUBROUTINE SWTT1_LMDAR4(KNU,KABS,KIND, PU, PTR) 2298 2208 USE dimphy 2209 USE radiation_AR4_param, only : APAD, BPAD, D 2299 2210 IMPLICIT none 2300 2211 cym#include "dimensions.h" … … 2345 2256 INTEGER jl, ja, i, j, ia 2346 2257 C 2347 C* Prescribed Data: 2348 C 2349 REAL(KIND=8) APAD(2,3,7), BPAD(2,3,7), D(2,3) 2350 SAVE APAD, BPAD, D 2351 c$OMP THREADPRIVATE(APAD, BPAD, D) 2352 DATA ((APAD(1,I,J),I=1,3),J=1,7) / 2353 S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, 2354 S 0.723613782E+05, 0.000000000E-00, 0.129353723E-01, 2355 S 0.596037057E+04, 0.000000000E-00, 0.800821928E+00, 2356 S 0.000000000E-00, 0.000000000E-00, 0.242715973E+02, 2357 S 0.000000000E-00, 0.000000000E-00, 0.878331486E+02, 2358 S 0.000000000E-00, 0.000000000E-00, 0.191559725E+02, 2359 S 0.000000000E-00, 0.000000000E-00, 0.000000000E+00 / 2360 DATA ((APAD(2,I,J),I=1,3),J=1,7) / 2361 S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, 2362 S 0.978576773E-04, 0.131849595E-03, 0.672595424E+02, 2363 S 0.387714006E+00, 0.437772681E+00, 0.000000000E-00, 2364 S 0.118461660E+03, 0.151345118E+03, 0.000000000E-00, 2365 S 0.119079797E+04, 0.233628890E+04, 0.000000000E-00, 2366 S 0.293353397E+03, 0.797219934E+03, 0.000000000E-00, 2367 S 0.000000000E+00, 0.000000000E+00, 0.000000000E+00 / 2368 C 2369 DATA ((BPAD(1,I,J),I=1,3),J=1,7) / 2370 S 0.912418292E+05, 0.000000000E-00, 0.925887084E-04, 2371 S 0.724555318E+05, 0.000000000E-00, 0.131812683E-01, 2372 S 0.602593328E+04, 0.000000000E-00, 0.812706117E+00, 2373 S 0.100000000E+01, 0.000000000E-00, 0.249863591E+02, 2374 S 0.000000000E-00, 0.000000000E-00, 0.931071925E+02, 2375 S 0.000000000E-00, 0.000000000E-00, 0.252233437E+02, 2376 S 0.000000000E-00, 0.000000000E-00, 0.100000000E+01 / 2377 DATA ((BPAD(2,I,J),I=1,3),J=1,7) / 2378 S 0.376655383E-08, 0.739646016E-08, 0.410177786E+03, 2379 S 0.979023421E-04, 0.131861712E-03, 0.731185438E+02, 2380 S 0.388611139E+00, 0.437949001E+00, 0.100000000E+01, 2381 S 0.120291383E+03, 0.151692730E+03, 0.000000000E+00, 2382 S 0.130531005E+04, 0.237071130E+04, 0.000000000E+00, 2383 S 0.415049409E+03, 0.867914360E+03, 0.000000000E+00, 2384 S 0.100000000E+01, 0.100000000E+01, 0.000000000E+00 / 2385 c 2386 DATA (D(1,I),I=1,3) / 0.00, 0.00, 0.00 / 2387 DATA (D(2,I),I=1,3) / 0.000000000, 0.000000000, 0.800000000 / 2258 2388 2259 C----------------------------------------------------------------------- 2389 2260 C … … 2614 2485 S PABCU) 2615 2486 USE dimphy 2487 USE radiation_AR4_param, only : TREF, RT1, RAER, AT, BT, OCT 2616 2488 IMPLICIT none 2617 2489 cym#include "dimensions.h" … … 2707 2579 REAL(KIND=8) R10E ! DECIMAL/NATURAL LOG.FACTOR 2708 2580 PARAMETER (R10E=0.4342945) 2709 c 2710 c Used Data Block: 2711 c 2712 REAL(KIND=8) TREF 2713 SAVE TREF 2714 c$OMP THREADPRIVATE(TREF) 2715 REAL(KIND=8) RT1(2) 2716 SAVE RT1 2717 c$OMP THREADPRIVATE(RT1) 2718 REAL(KIND=8) RAER(5,5) 2719 SAVE RAER 2720 c$OMP THREADPRIVATE(RAER) 2721 REAL(KIND=8) AT(8,3), BT(8,3) 2722 SAVE AT, BT 2723 c$OMP THREADPRIVATE(AT, BT) 2724 REAL(KIND=8) OCT(4) 2725 SAVE OCT 2726 c$OMP THREADPRIVATE(OCT) 2727 DATA TREF /250.0/ 2728 DATA (RT1(IG1),IG1=1,2) / -0.577350269, +0.577350269 / 2729 DATA RAER / .038520, .037196, .040532, .054934, .038520 2730 1 , .12613 , .18313 , .10357 , .064106, .126130 2731 2 , .012579, .013649, .018652, .025181, .012579 2732 3 , .011890, .016142, .021105, .028908, .011890 2733 4 , .013792, .026810, .052203, .066338, .013792 / 2734 DATA (AT(1,IR),IR=1,3) / 2735 S 0.298199E-02,-.394023E-03,0.319566E-04 / 2736 DATA (BT(1,IR),IR=1,3) / 2737 S-0.106432E-04,0.660324E-06,0.174356E-06 / 2738 DATA (AT(2,IR),IR=1,3) / 2739 S 0.143676E-01,0.366501E-02,-.160822E-02 / 2740 DATA (BT(2,IR),IR=1,3) / 2741 S-0.553979E-04,-.101701E-04,0.920868E-05 / 2742 DATA (AT(3,IR),IR=1,3) / 2743 S 0.197861E-01,0.315541E-02,-.174547E-02 / 2744 DATA (BT(3,IR),IR=1,3) / 2745 S-0.877012E-04,0.513302E-04,0.523138E-06 / 2746 DATA (AT(4,IR),IR=1,3) / 2747 S 0.289560E-01,-.208807E-02,-.121943E-02 / 2748 DATA (BT(4,IR),IR=1,3) / 2749 S-0.165960E-03,0.157704E-03,-.146427E-04 / 2750 DATA (AT(5,IR),IR=1,3) / 2751 S 0.103800E-01,0.436296E-02,-.161431E-02 / 2752 DATA (BT(5,IR),IR=1,3) / 2753 S -.276744E-04,-.327381E-04,0.127646E-04 / 2754 DATA (AT(6,IR),IR=1,3) / 2755 S 0.868859E-02,-.972752E-03,0.000000E-00 / 2756 DATA (BT(6,IR),IR=1,3) / 2757 S -.278412E-04,-.713940E-06,0.117469E-05 / 2758 DATA (AT(7,IR),IR=1,3) / 2759 S 0.250073E-03,0.455875E-03,0.109242E-03 / 2760 DATA (BT(7,IR),IR=1,3) / 2761 S 0.199846E-05,-.216313E-05,0.175991E-06 / 2762 DATA (AT(8,IR),IR=1,3) / 2763 S 0.307423E-01,0.110879E-02,-.322172E-03 / 2764 DATA (BT(8,IR),IR=1,3) / 2765 S-0.108482E-03,0.258096E-05,-.814575E-06 / 2766 c 2767 DATA OCT /-.326E-03, -.102E-05, .137E-02, -.535E-05/ 2581 2768 2582 C----------------------------------------------------------------------- 2769 2583 c … … 3510 3324 S , PGA,PGB,PGASUR,PGBSUR,PGATOP,PGBTOP) 3511 3325 USE dimphy 3326 USE radiation_AR4_param, only : TINTP, XP, GA, GB 3512 3327 IMPLICIT none 3513 3328 cym#include "dimensions.h" … … 3604 3419 C 3605 3420 C* Used Data Block: 3606 REAL(KIND=8)TINTP(11)3607 3421 c REAL*8 TINTP(11) 3422 c SAVE TINTP 3608 3423 c$OMP THREADPRIVATE(TINTP) 3609 REAL(KIND=8)GA(11,16,3), GB(11,16,3)3610 3424 c REAL*8 GA(11,16,3), GB(11,16,3) 3425 c SAVE GA, GB 3611 3426 c$OMP THREADPRIVATE(GA, GB) 3612 REAL(KIND=8)XP(6,6)3613 3427 c REAL*8 XP(6,6) 3428 c SAVE XP 3614 3429 c$OMP THREADPRIVATE(XP) 3615 3430 c 3616 3617 3431 c DATA TINTP / 187.5, 200., 212.5, 225., 237.5, 250., 3432 c S 262.5, 275., 287.5, 300., 312.5 / 3618 3433 C----------------------------------------------------------------------- 3619 3434 C-- WATER VAPOR -- INT.1 -- 0- 500 CM-1 -- FROM ABS225 ---------------- … … 3628 3443 C 3629 3444 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3630 3631 3632 3633 3634 3635 3636 3637 3445 C DATA (GA( 1, 1,IC),IC=1,3) / 3446 C S 0.63499072E-02,-0.99506586E-03, 0.00000000E+00/ 3447 C DATA (GB( 1, 1,IC),IC=1,3) / 3448 C S 0.63499072E-02, 0.97222852E-01, 0.10000000E+01/ 3449 C DATA (GA( 1, 2,IC),IC=1,3) / 3450 C S 0.77266491E-02,-0.11661515E-02, 0.00000000E+00/ 3451 C DATA (GB( 1, 2,IC),IC=1,3) / 3452 C S 0.77266491E-02, 0.10681591E+00, 0.10000000E+01/ 3638 3453 C 3639 3454 C----- INTERVAL = 1 ----- T = 200.0 3640 3455 C 3641 3456 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3642 3643 3644 3645 3646 3647 3648 3649 3457 C DATA (GA( 2, 1,IC),IC=1,3) / 3458 C S 0.65566348E-02,-0.10184169E-02, 0.00000000E+00/ 3459 C DATA (GB( 2, 1,IC),IC=1,3) / 3460 C S 0.65566348E-02, 0.98862238E-01, 0.10000000E+01/ 3461 C DATA (GA( 2, 2,IC),IC=1,3) / 3462 C S 0.81323287E-02,-0.11886130E-02, 0.00000000E+00/ 3463 C DATA (GB( 2, 2,IC),IC=1,3) / 3464 C S 0.81323287E-02, 0.10921298E+00, 0.10000000E+01/ 3650 3465 C 3651 3466 C----- INTERVAL = 1 ----- T = 212.5 3652 3467 C 3653 3468 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3654 3655 3656 3657 3658 3659 3660 3661 3469 C DATA (GA( 3, 1,IC),IC=1,3) / 3470 C S 0.67849730E-02,-0.10404730E-02, 0.00000000E+00/ 3471 C DATA (GB( 3, 1,IC),IC=1,3) / 3472 C S 0.67849730E-02, 0.10061504E+00, 0.10000000E+01/ 3473 C DATA (GA( 3, 2,IC),IC=1,3) / 3474 C S 0.86507620E-02,-0.12139929E-02, 0.00000000E+00/ 3475 C DATA (GB( 3, 2,IC),IC=1,3) / 3476 C S 0.86507620E-02, 0.11198225E+00, 0.10000000E+01/ 3662 3477 C 3663 3478 C----- INTERVAL = 1 ----- T = 225.0 3664 3479 C 3665 3480 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3666 3667 3668 3669 3670 3671 3672 3673 3481 C DATA (GA( 4, 1,IC),IC=1,3) / 3482 C S 0.70481947E-02,-0.10621792E-02, 0.00000000E+00/ 3483 C DATA (GB( 4, 1,IC),IC=1,3) / 3484 C S 0.70481947E-02, 0.10256222E+00, 0.10000000E+01/ 3485 C DATA (GA( 4, 2,IC),IC=1,3) / 3486 C S 0.92776391E-02,-0.12445811E-02, 0.00000000E+00/ 3487 C DATA (GB( 4, 2,IC),IC=1,3) / 3488 C S 0.92776391E-02, 0.11487826E+00, 0.10000000E+01/ 3674 3489 C 3675 3490 C----- INTERVAL = 1 ----- T = 237.5 3676 3491 C 3677 3492 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3678 3679 3680 3681 3682 3683 3684 3685 3493 C DATA (GA( 5, 1,IC),IC=1,3) / 3494 C S 0.73585943E-02,-0.10847662E-02, 0.00000000E+00/ 3495 C DATA (GB( 5, 1,IC),IC=1,3) / 3496 C S 0.73585943E-02, 0.10475952E+00, 0.10000000E+01/ 3497 C DATA (GA( 5, 2,IC),IC=1,3) / 3498 C S 0.99806312E-02,-0.12807672E-02, 0.00000000E+00/ 3499 C DATA (GB( 5, 2,IC),IC=1,3) / 3500 C S 0.99806312E-02, 0.11751113E+00, 0.10000000E+01/ 3686 3501 C 3687 3502 C----- INTERVAL = 1 ----- T = 250.0 3688 3503 C 3689 3504 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3690 3691 3692 3693 3694 3695 3696 3697 3505 C DATA (GA( 6, 1,IC),IC=1,3) / 3506 C S 0.77242818E-02,-0.11094726E-02, 0.00000000E+00/ 3507 C DATA (GB( 6, 1,IC),IC=1,3) / 3508 C S 0.77242818E-02, 0.10720986E+00, 0.10000000E+01/ 3509 C DATA (GA( 6, 2,IC),IC=1,3) / 3510 C S 0.10709803E-01,-0.13208251E-02, 0.00000000E+00/ 3511 C DATA (GB( 6, 2,IC),IC=1,3) / 3512 C S 0.10709803E-01, 0.11951535E+00, 0.10000000E+01/ 3698 3513 C 3699 3514 C----- INTERVAL = 1 ----- T = 262.5 3700 3515 C 3701 3516 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3702 3703 3704 3705 3706 3707 3708 3709 3517 C DATA (GA( 7, 1,IC),IC=1,3) / 3518 C S 0.81472693E-02,-0.11372949E-02, 0.00000000E+00/ 3519 C DATA (GB( 7, 1,IC),IC=1,3) / 3520 C S 0.81472693E-02, 0.10985370E+00, 0.10000000E+01/ 3521 C DATA (GA( 7, 2,IC),IC=1,3) / 3522 C S 0.11414739E-01,-0.13619034E-02, 0.00000000E+00/ 3523 C DATA (GB( 7, 2,IC),IC=1,3) / 3524 C S 0.11414739E-01, 0.12069945E+00, 0.10000000E+01/ 3710 3525 C 3711 3526 C----- INTERVAL = 1 ----- T = 275.0 3712 3527 C 3713 3528 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3714 3715 3716 3717 3718 3719 3720 3721 3529 C DATA (GA( 8, 1,IC),IC=1,3) / 3530 C S 0.86227527E-02,-0.11687683E-02, 0.00000000E+00/ 3531 C DATA (GB( 8, 1,IC),IC=1,3) / 3532 C S 0.86227527E-02, 0.11257633E+00, 0.10000000E+01/ 3533 C DATA (GA( 8, 2,IC),IC=1,3) / 3534 C S 0.12058772E-01,-0.14014165E-02, 0.00000000E+00/ 3535 C DATA (GB( 8, 2,IC),IC=1,3) / 3536 C S 0.12058772E-01, 0.12108524E+00, 0.10000000E+01/ 3722 3537 C 3723 3538 C----- INTERVAL = 1 ----- T = 287.5 3724 3539 C 3725 3540 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3726 3727 3728 3729 3730 3731 3732 3733 3541 C DATA (GA( 9, 1,IC),IC=1,3) / 3542 C S 0.91396814E-02,-0.12038314E-02, 0.00000000E+00/ 3543 C DATA (GB( 9, 1,IC),IC=1,3) / 3544 C S 0.91396814E-02, 0.11522980E+00, 0.10000000E+01/ 3545 C DATA (GA( 9, 2,IC),IC=1,3) / 3546 C S 0.12623992E-01,-0.14378639E-02, 0.00000000E+00/ 3547 C DATA (GB( 9, 2,IC),IC=1,3) / 3548 C S 0.12623992E-01, 0.12084229E+00, 0.10000000E+01/ 3734 3549 C 3735 3550 C----- INTERVAL = 1 ----- T = 300.0 3736 3551 C 3737 3552 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3738 3739 3740 3741 3742 3743 3744 3745 3553 C DATA (GA(10, 1,IC),IC=1,3) / 3554 C S 0.96825438E-02,-0.12418367E-02, 0.00000000E+00/ 3555 C DATA (GB(10, 1,IC),IC=1,3) / 3556 C S 0.96825438E-02, 0.11766343E+00, 0.10000000E+01/ 3557 C DATA (GA(10, 2,IC),IC=1,3) / 3558 C S 0.13108146E-01,-0.14708488E-02, 0.00000000E+00/ 3559 C DATA (GB(10, 2,IC),IC=1,3) / 3560 C S 0.13108146E-01, 0.12019005E+00, 0.10000000E+01/ 3746 3561 C 3747 3562 C----- INTERVAL = 1 ----- T = 312.5 3748 3563 C 3749 3564 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 3750 3751 3752 3753 3754 3755 3756 3757 3565 C DATA (GA(11, 1,IC),IC=1,3) / 3566 C S 0.10233955E-01,-0.12817135E-02, 0.00000000E+00/ 3567 C DATA (GB(11, 1,IC),IC=1,3) / 3568 C S 0.10233955E-01, 0.11975320E+00, 0.10000000E+01/ 3569 C DATA (GA(11, 2,IC),IC=1,3) / 3570 C S 0.13518390E-01,-0.15006791E-02, 0.00000000E+00/ 3571 C DATA (GB(11, 2,IC),IC=1,3) / 3572 C S 0.13518390E-01, 0.11932684E+00, 0.10000000E+01/ 3758 3573 C 3759 3574 C … … 3770 3585 C 3771 3586 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3772 3773 3774 3775 3776 3777 3778 3779 3587 C DATA (GA( 1, 3,IC),IC=1,3) / 3588 C S 0.11644593E+01, 0.41243390E+00, 0.00000000E+00/ 3589 C DATA (GB( 1, 3,IC),IC=1,3) / 3590 C S 0.11644593E+01, 0.10346097E+01, 0.10000000E+01/ 3591 C DATA (GA( 1, 4,IC),IC=1,3) / 3592 C S 0.12006968E+01, 0.48318936E+00, 0.00000000E+00/ 3593 C DATA (GB( 1, 4,IC),IC=1,3) / 3594 C S 0.12006968E+01, 0.10626130E+01, 0.10000000E+01/ 3780 3595 C 3781 3596 C----- INTERVAL = 2 ----- T = 200.0 3782 3597 C 3783 3598 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3784 3785 3786 3787 3788 3789 3790 3791 3599 C DATA (GA( 2, 3,IC),IC=1,3) / 3600 C S 0.11747203E+01, 0.43407282E+00, 0.00000000E+00/ 3601 C DATA (GB( 2, 3,IC),IC=1,3) / 3602 C S 0.11747203E+01, 0.10433655E+01, 0.10000000E+01/ 3603 C DATA (GA( 2, 4,IC),IC=1,3) / 3604 C S 0.12108196E+01, 0.50501827E+00, 0.00000000E+00/ 3605 C DATA (GB( 2, 4,IC),IC=1,3) / 3606 C S 0.12108196E+01, 0.10716026E+01, 0.10000000E+01/ 3792 3607 C 3793 3608 C----- INTERVAL = 2 ----- T = 212.5 3794 3609 C 3795 3610 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3796 3797 3798 3799 3800 3801 3802 3803 3611 C DATA (GA( 3, 3,IC),IC=1,3) / 3612 C S 0.11837872E+01, 0.45331413E+00, 0.00000000E+00/ 3613 C DATA (GB( 3, 3,IC),IC=1,3) / 3614 C S 0.11837872E+01, 0.10511933E+01, 0.10000000E+01/ 3615 C DATA (GA( 3, 4,IC),IC=1,3) / 3616 C S 0.12196717E+01, 0.52409502E+00, 0.00000000E+00/ 3617 C DATA (GB( 3, 4,IC),IC=1,3) / 3618 C S 0.12196717E+01, 0.10795108E+01, 0.10000000E+01/ 3804 3619 C 3805 3620 C----- INTERVAL = 2 ----- T = 225.0 3806 3621 C 3807 3622 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3808 3809 3810 3811 3812 3813 3814 3815 3623 C DATA (GA( 4, 3,IC),IC=1,3) / 3624 C S 0.11918561E+01, 0.47048604E+00, 0.00000000E+00/ 3625 C DATA (GB( 4, 3,IC),IC=1,3) / 3626 C S 0.11918561E+01, 0.10582150E+01, 0.10000000E+01/ 3627 C DATA (GA( 4, 4,IC),IC=1,3) / 3628 C S 0.12274493E+01, 0.54085277E+00, 0.00000000E+00/ 3629 C DATA (GB( 4, 4,IC),IC=1,3) / 3630 C S 0.12274493E+01, 0.10865006E+01, 0.10000000E+01/ 3816 3631 C 3817 3632 C----- INTERVAL = 2 ----- T = 237.5 3818 3633 C 3819 3634 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3820 3821 3822 3823 3824 3825 3826 3827 3635 C DATA (GA( 5, 3,IC),IC=1,3) / 3636 C S 0.11990757E+01, 0.48586286E+00, 0.00000000E+00/ 3637 C DATA (GB( 5, 3,IC),IC=1,3) / 3638 C S 0.11990757E+01, 0.10645317E+01, 0.10000000E+01/ 3639 C DATA (GA( 5, 4,IC),IC=1,3) / 3640 C S 0.12343189E+01, 0.55565422E+00, 0.00000000E+00/ 3641 C DATA (GB( 5, 4,IC),IC=1,3) / 3642 C S 0.12343189E+01, 0.10927103E+01, 0.10000000E+01/ 3828 3643 C 3829 3644 C----- INTERVAL = 2 ----- T = 250.0 3830 3645 C 3831 3646 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3832 3833 3834 3835 3836 3837 3838 3839 3647 C DATA (GA( 6, 3,IC),IC=1,3) / 3648 C S 0.12055643E+01, 0.49968044E+00, 0.00000000E+00/ 3649 C DATA (GB( 6, 3,IC),IC=1,3) / 3650 C S 0.12055643E+01, 0.10702313E+01, 0.10000000E+01/ 3651 C DATA (GA( 6, 4,IC),IC=1,3) / 3652 C S 0.12404147E+01, 0.56878618E+00, 0.00000000E+00/ 3653 C DATA (GB( 6, 4,IC),IC=1,3) / 3654 C S 0.12404147E+01, 0.10982489E+01, 0.10000000E+01/ 3840 3655 C 3841 3656 C----- INTERVAL = 2 ----- T = 262.5 3842 3657 C 3843 3658 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3844 3845 3846 3847 3848 3849 3850 3851 3659 C DATA (GA( 7, 3,IC),IC=1,3) / 3660 C S 0.12114186E+01, 0.51214132E+00, 0.00000000E+00/ 3661 C DATA (GB( 7, 3,IC),IC=1,3) / 3662 C S 0.12114186E+01, 0.10753907E+01, 0.10000000E+01/ 3663 C DATA (GA( 7, 4,IC),IC=1,3) / 3664 C S 0.12458431E+01, 0.58047395E+00, 0.00000000E+00/ 3665 C DATA (GB( 7, 4,IC),IC=1,3) / 3666 C S 0.12458431E+01, 0.11032019E+01, 0.10000000E+01/ 3852 3667 C 3853 3668 C----- INTERVAL = 2 ----- T = 275.0 3854 3669 C 3855 3670 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3856 3857 3858 3859 3860 3861 3862 3863 3671 C DATA (GA( 8, 3,IC),IC=1,3) / 3672 C S 0.12167192E+01, 0.52341830E+00, 0.00000000E+00/ 3673 C DATA (GB( 8, 3,IC),IC=1,3) / 3674 C S 0.12167192E+01, 0.10800762E+01, 0.10000000E+01/ 3675 C DATA (GA( 8, 4,IC),IC=1,3) / 3676 C S 0.12506907E+01, 0.59089894E+00, 0.00000000E+00/ 3677 C DATA (GB( 8, 4,IC),IC=1,3) / 3678 C S 0.12506907E+01, 0.11076379E+01, 0.10000000E+01/ 3864 3679 C 3865 3680 C----- INTERVAL = 2 ----- T = 287.5 3866 3681 C 3867 3682 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3868 3869 3870 3871 3872 3873 3874 3875 3683 C DATA (GA( 9, 3,IC),IC=1,3) / 3684 C S 0.12215344E+01, 0.53365803E+00, 0.00000000E+00/ 3685 C DATA (GB( 9, 3,IC),IC=1,3) / 3686 C S 0.12215344E+01, 0.10843446E+01, 0.10000000E+01/ 3687 C DATA (GA( 9, 4,IC),IC=1,3) / 3688 C S 0.12550299E+01, 0.60021475E+00, 0.00000000E+00/ 3689 C DATA (GB( 9, 4,IC),IC=1,3) / 3690 C S 0.12550299E+01, 0.11116160E+01, 0.10000000E+01/ 3876 3691 C 3877 3692 C----- INTERVAL = 2 ----- T = 300.0 3878 3693 C 3879 3694 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3880 3881 3882 3883 3884 3885 3886 3887 3695 C DATA (GA(10, 3,IC),IC=1,3) / 3696 C S 0.12259226E+01, 0.54298448E+00, 0.00000000E+00/ 3697 C DATA (GB(10, 3,IC),IC=1,3) / 3698 C S 0.12259226E+01, 0.10882439E+01, 0.10000000E+01/ 3699 C DATA (GA(10, 4,IC),IC=1,3) / 3700 C S 0.12589256E+01, 0.60856112E+00, 0.00000000E+00/ 3701 C DATA (GB(10, 4,IC),IC=1,3) / 3702 C S 0.12589256E+01, 0.11151910E+01, 0.10000000E+01/ 3888 3703 C 3889 3704 C----- INTERVAL = 2 ----- T = 312.5 3890 3705 C 3891 3706 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3892 3893 3894 3895 3896 3897 3898 3899 3707 C DATA (GA(11, 3,IC),IC=1,3) / 3708 C S 0.12299344E+01, 0.55150227E+00, 0.00000000E+00/ 3709 C DATA (GB(11, 3,IC),IC=1,3) / 3710 C S 0.12299344E+01, 0.10918144E+01, 0.10000000E+01/ 3711 C DATA (GA(11, 4,IC),IC=1,3) / 3712 C S 0.12624402E+01, 0.61607594E+00, 0.00000000E+00/ 3713 C DATA (GB(11, 4,IC),IC=1,3) / 3714 C S 0.12624402E+01, 0.11184188E+01, 0.10000000E+01/ 3900 3715 C 3901 3716 C … … 3916 3731 C 3917 3732 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3918 3919 3920 3921 3922 3923 3924 3925 3733 C DATA (GA( 1, 7,IC),IC=1,3) / 3734 C S 0.10192131E+02, 0.80737799E+01, 0.00000000E+00/ 3735 C DATA (GB( 1, 7,IC),IC=1,3) / 3736 C S 0.10192131E+02, 0.82623280E+01, 0.10000000E+01/ 3737 C DATA (GA( 1, 8,IC),IC=1,3) / 3738 C S 0.92439050E+01, 0.77425778E+01, 0.00000000E+00/ 3739 C DATA (GB( 1, 8,IC),IC=1,3) / 3740 C S 0.92439050E+01, 0.79342219E+01, 0.10000000E+01/ 3926 3741 C 3927 3742 C----- INTERVAL = 3 ----- T = 200.0 3928 3743 C 3929 3744 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3930 3931 3932 3933 3934 3935 3936 3937 3745 C DATA (GA( 2, 7,IC),IC=1,3) / 3746 C S 0.97258602E+01, 0.79171158E+01, 0.00000000E+00/ 3747 C DATA (GB( 2, 7,IC),IC=1,3) / 3748 C S 0.97258602E+01, 0.81072291E+01, 0.10000000E+01/ 3749 C DATA (GA( 2, 8,IC),IC=1,3) / 3750 C S 0.87567422E+01, 0.75443460E+01, 0.00000000E+00/ 3751 C DATA (GB( 2, 8,IC),IC=1,3) / 3752 C S 0.87567422E+01, 0.77373458E+01, 0.10000000E+01/ 3938 3753 C 3939 3754 C----- INTERVAL = 3 ----- T = 212.5 3940 3755 C 3941 3756 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3942 3943 3944 3945 3946 3947 3948 3949 3757 C DATA (GA( 3, 7,IC),IC=1,3) / 3758 C S 0.92992890E+01, 0.77609605E+01, 0.00000000E+00/ 3759 C DATA (GB( 3, 7,IC),IC=1,3) / 3760 C S 0.92992890E+01, 0.79523834E+01, 0.10000000E+01/ 3761 C DATA (GA( 3, 8,IC),IC=1,3) / 3762 C S 0.83270144E+01, 0.73526151E+01, 0.00000000E+00/ 3763 C DATA (GB( 3, 8,IC),IC=1,3) / 3764 C S 0.83270144E+01, 0.75467334E+01, 0.10000000E+01/ 3950 3765 C 3951 3766 C----- INTERVAL = 3 ----- T = 225.0 3952 3767 C 3953 3768 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3954 3955 3956 3957 3958 3959 3960 3961 3769 C DATA (GA( 4, 7,IC),IC=1,3) / 3770 C S 0.89154021E+01, 0.76087371E+01, 0.00000000E+00/ 3771 C DATA (GB( 4, 7,IC),IC=1,3) / 3772 C S 0.89154021E+01, 0.78012527E+01, 0.10000000E+01/ 3773 C DATA (GA( 4, 8,IC),IC=1,3) / 3774 C S 0.79528337E+01, 0.71711188E+01, 0.00000000E+00/ 3775 C DATA (GB( 4, 8,IC),IC=1,3) / 3776 C S 0.79528337E+01, 0.73661786E+01, 0.10000000E+01/ 3962 3777 C 3963 3778 C----- INTERVAL = 3 ----- T = 237.5 3964 3779 C 3965 3780 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3966 3967 3968 3969 3970 3971 3972 3973 3781 C DATA (GA( 5, 7,IC),IC=1,3) / 3782 C S 0.85730084E+01, 0.74627112E+01, 0.00000000E+00/ 3783 C DATA (GB( 5, 7,IC),IC=1,3) / 3784 C S 0.85730084E+01, 0.76561458E+01, 0.10000000E+01/ 3785 C DATA (GA( 5, 8,IC),IC=1,3) / 3786 C S 0.76286839E+01, 0.70015571E+01, 0.00000000E+00/ 3787 C DATA (GB( 5, 8,IC),IC=1,3) / 3788 C S 0.76286839E+01, 0.71974319E+01, 0.10000000E+01/ 3974 3789 C 3975 3790 C----- INTERVAL = 3 ----- T = 250.0 3976 3791 C 3977 3792 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3978 3979 3980 3981 3982 3983 3984 3985 3793 C DATA (GA( 6, 7,IC),IC=1,3) / 3794 C S 0.82685838E+01, 0.73239981E+01, 0.00000000E+00/ 3795 C DATA (GB( 6, 7,IC),IC=1,3) / 3796 C S 0.82685838E+01, 0.75182174E+01, 0.10000000E+01/ 3797 C DATA (GA( 6, 8,IC),IC=1,3) / 3798 C S 0.73477879E+01, 0.68442532E+01, 0.00000000E+00/ 3799 C DATA (GB( 6, 8,IC),IC=1,3) / 3800 C S 0.73477879E+01, 0.70408543E+01, 0.10000000E+01/ 3986 3801 C 3987 3802 C----- INTERVAL = 3 ----- T = 262.5 3988 3803 C 3989 3804 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 3990 3991 3992 3993 3994 3995 3996 3997 3805 C DATA (GA( 7, 7,IC),IC=1,3) / 3806 C S 0.79978921E+01, 0.71929934E+01, 0.00000000E+00/ 3807 C DATA (GB( 7, 7,IC),IC=1,3) / 3808 C S 0.79978921E+01, 0.73878952E+01, 0.10000000E+01/ 3809 C DATA (GA( 7, 8,IC),IC=1,3) / 3810 C S 0.71035818E+01, 0.66987996E+01, 0.00000000E+00/ 3811 C DATA (GB( 7, 8,IC),IC=1,3) / 3812 C S 0.71035818E+01, 0.68960649E+01, 0.10000000E+01/ 3998 3813 C 3999 3814 C----- INTERVAL = 3 ----- T = 275.0 4000 3815 C 4001 3816 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4002 4003 4004 4005 4006 4007 4008 4009 3817 C DATA (GA( 8, 7,IC),IC=1,3) / 3818 C S 0.77568055E+01, 0.70697065E+01, 0.00000000E+00/ 3819 C DATA (GB( 8, 7,IC),IC=1,3) / 3820 C S 0.77568055E+01, 0.72652133E+01, 0.10000000E+01/ 3821 C DATA (GA( 8, 8,IC),IC=1,3) / 3822 C S 0.68903312E+01, 0.65644820E+01, 0.00000000E+00/ 3823 C DATA (GB( 8, 8,IC),IC=1,3) / 3824 C S 0.68903312E+01, 0.67623672E+01, 0.10000000E+01/ 4010 3825 C 4011 3826 C----- INTERVAL = 3 ----- T = 287.5 4012 3827 C 4013 3828 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4014 4015 4016 4017 4018 4019 4020 4021 3829 C DATA (GA( 9, 7,IC),IC=1,3) / 3830 C S 0.75416266E+01, 0.69539626E+01, 0.00000000E+00/ 3831 C DATA (GB( 9, 7,IC),IC=1,3) / 3832 C S 0.75416266E+01, 0.71500151E+01, 0.10000000E+01/ 3833 C DATA (GA( 9, 8,IC),IC=1,3) / 3834 C S 0.67032875E+01, 0.64405267E+01, 0.00000000E+00/ 3835 C DATA (GB( 9, 8,IC),IC=1,3) / 3836 C S 0.67032875E+01, 0.66389989E+01, 0.10000000E+01/ 4022 3837 C 4023 3838 C----- INTERVAL = 3 ----- T = 300.0 4024 3839 C 4025 3840 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4026 4027 4028 4029 4030 4031 4032 4033 3841 C DATA (GA(10, 7,IC),IC=1,3) / 3842 C S 0.73491694E+01, 0.68455144E+01, 0.00000000E+00/ 3843 C DATA (GB(10, 7,IC),IC=1,3) / 3844 C S 0.73491694E+01, 0.70420667E+01, 0.10000000E+01/ 3845 C DATA (GA(10, 8,IC),IC=1,3) / 3846 C S 0.65386461E+01, 0.63262376E+01, 0.00000000E+00/ 3847 C DATA (GB(10, 8,IC),IC=1,3) / 3848 C S 0.65386461E+01, 0.65252707E+01, 0.10000000E+01/ 4034 3849 C 4035 3850 C----- INTERVAL = 3 ----- T = 312.5 4036 3851 C 4037 3852 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4038 4039 4040 4041 4042 4043 4044 4045 3853 C DATA (GA(11, 7,IC),IC=1,3) / 3854 C S 0.71767400E+01, 0.67441020E+01, 0.00000000E+00/ 3855 C DATA (GB(11, 7,IC),IC=1,3) / 3856 C S 0.71767400E+01, 0.69411177E+01, 0.10000000E+01/ 3857 C DATA (GA(11, 8,IC),IC=1,3) / 3858 C S 0.63934377E+01, 0.62210701E+01, 0.00000000E+00/ 3859 C DATA (GB(11, 8,IC),IC=1,3) / 3860 C S 0.63934377E+01, 0.64206412E+01, 0.10000000E+01/ 4046 3861 C 4047 3862 C … … 4053 3868 C 4054 3869 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4055 4056 4057 4058 4059 4060 4061 4062 3870 C DATA (GA( 1, 9,IC),IC=1,3) / 3871 C S 0.24870635E+02, 0.10542131E+02, 0.00000000E+00/ 3872 C DATA (GB( 1, 9,IC),IC=1,3) / 3873 C S 0.24870635E+02, 0.10656640E+02, 0.10000000E+01/ 3874 C DATA (GA( 1,10,IC),IC=1,3) / 3875 C S 0.24586283E+02, 0.10490353E+02, 0.00000000E+00/ 3876 C DATA (GB( 1,10,IC),IC=1,3) / 3877 C S 0.24586283E+02, 0.10605856E+02, 0.10000000E+01/ 4063 3878 C 4064 3879 C----- INTERVAL = 4 ----- T = 200.0 4065 3880 C 4066 3881 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4067 4068 4069 4070 4071 4072 4073 4074 3882 C DATA (GA( 2, 9,IC),IC=1,3) / 3883 C S 0.24725591E+02, 0.10515895E+02, 0.00000000E+00/ 3884 C DATA (GB( 2, 9,IC),IC=1,3) / 3885 C S 0.24725591E+02, 0.10630910E+02, 0.10000000E+01/ 3886 C DATA (GA( 2,10,IC),IC=1,3) / 3887 C S 0.24441465E+02, 0.10463512E+02, 0.00000000E+00/ 3888 C DATA (GB( 2,10,IC),IC=1,3) / 3889 C S 0.24441465E+02, 0.10579514E+02, 0.10000000E+01/ 4075 3890 C 4076 3891 C----- INTERVAL = 4 ----- T = 212.5 4077 3892 C 4078 3893 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4079 4080 4081 4082 4083 4084 4085 4086 3894 C DATA (GA( 3, 9,IC),IC=1,3) / 3895 C S 0.24600320E+02, 0.10492949E+02, 0.00000000E+00/ 3896 C DATA (GB( 3, 9,IC),IC=1,3) / 3897 C S 0.24600320E+02, 0.10608399E+02, 0.10000000E+01/ 3898 C DATA (GA( 3,10,IC),IC=1,3) / 3899 C S 0.24311657E+02, 0.10439183E+02, 0.00000000E+00/ 3900 C DATA (GB( 3,10,IC),IC=1,3) / 3901 C S 0.24311657E+02, 0.10555632E+02, 0.10000000E+01/ 4087 3902 C 4088 3903 C----- INTERVAL = 4 ----- T = 225.0 4089 3904 C 4090 3905 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4091 4092 4093 4094 4095 4096 4097 4098 3906 C DATA (GA( 4, 9,IC),IC=1,3) / 3907 C S 0.24487300E+02, 0.10472049E+02, 0.00000000E+00/ 3908 C DATA (GB( 4, 9,IC),IC=1,3) / 3909 C S 0.24487300E+02, 0.10587891E+02, 0.10000000E+01/ 3910 C DATA (GA( 4,10,IC),IC=1,3) / 3911 C S 0.24196167E+02, 0.10417324E+02, 0.00000000E+00/ 3912 C DATA (GB( 4,10,IC),IC=1,3) / 3913 C S 0.24196167E+02, 0.10534169E+02, 0.10000000E+01/ 4099 3914 C 4100 3915 C----- INTERVAL = 4 ----- T = 237.5 4101 3916 C 4102 3917 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4103 4104 4105 4106 4107 4108 4109 4110 3918 C DATA (GA( 5, 9,IC),IC=1,3) / 3919 C S 0.24384935E+02, 0.10452961E+02, 0.00000000E+00/ 3920 C DATA (GB( 5, 9,IC),IC=1,3) / 3921 C S 0.24384935E+02, 0.10569156E+02, 0.10000000E+01/ 3922 C DATA (GA( 5,10,IC),IC=1,3) / 3923 C S 0.24093406E+02, 0.10397704E+02, 0.00000000E+00/ 3924 C DATA (GB( 5,10,IC),IC=1,3) / 3925 C S 0.24093406E+02, 0.10514900E+02, 0.10000000E+01/ 4111 3926 C 4112 3927 C----- INTERVAL = 4 ----- T = 250.0 4113 3928 C 4114 3929 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4115 4116 4117 4118 4119 4120 4121 4122 3930 C DATA (GA( 6, 9,IC),IC=1,3) / 3931 C S 0.24292341E+02, 0.10435562E+02, 0.00000000E+00/ 3932 C DATA (GB( 6, 9,IC),IC=1,3) / 3933 C S 0.24292341E+02, 0.10552075E+02, 0.10000000E+01/ 3934 C DATA (GA( 6,10,IC),IC=1,3) / 3935 C S 0.24001597E+02, 0.10380038E+02, 0.00000000E+00/ 3936 C DATA (GB( 6,10,IC),IC=1,3) / 3937 C S 0.24001597E+02, 0.10497547E+02, 0.10000000E+01/ 4123 3938 C 4124 3939 C----- INTERVAL = 4 ----- T = 262.5 4125 3940 C 4126 3941 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4127 4128 4129 4130 4131 4132 4133 4134 3942 C DATA (GA( 7, 9,IC),IC=1,3) / 3943 C S 0.24208572E+02, 0.10419710E+02, 0.00000000E+00/ 3944 C DATA (GB( 7, 9,IC),IC=1,3) / 3945 C S 0.24208572E+02, 0.10536510E+02, 0.10000000E+01/ 3946 C DATA (GA( 7,10,IC),IC=1,3) / 3947 C S 0.23919098E+02, 0.10364052E+02, 0.00000000E+00/ 3948 C DATA (GB( 7,10,IC),IC=1,3) / 3949 C S 0.23919098E+02, 0.10481842E+02, 0.10000000E+01/ 4135 3950 C 4136 3951 C----- INTERVAL = 4 ----- T = 275.0 4137 3952 C 4138 3953 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4139 4140 4141 4142 4143 4144 4145 4146 3954 C DATA (GA( 8, 9,IC),IC=1,3) / 3955 C S 0.24132642E+02, 0.10405247E+02, 0.00000000E+00/ 3956 C DATA (GB( 8, 9,IC),IC=1,3) / 3957 C S 0.24132642E+02, 0.10522307E+02, 0.10000000E+01/ 3958 C DATA (GA( 8,10,IC),IC=1,3) / 3959 C S 0.23844511E+02, 0.10349509E+02, 0.00000000E+00/ 3960 C DATA (GB( 8,10,IC),IC=1,3) / 3961 C S 0.23844511E+02, 0.10467553E+02, 0.10000000E+01/ 4147 3962 C 4148 3963 C----- INTERVAL = 4 ----- T = 287.5 4149 3964 C 4150 3965 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4151 4152 4153 4154 4155 4156 4157 4158 3966 C DATA (GA( 9, 9,IC),IC=1,3) / 3967 C S 0.24063614E+02, 0.10392022E+02, 0.00000000E+00/ 3968 C DATA (GB( 9, 9,IC),IC=1,3) / 3969 C S 0.24063614E+02, 0.10509317E+02, 0.10000000E+01/ 3970 C DATA (GA( 9,10,IC),IC=1,3) / 3971 C S 0.23776708E+02, 0.10336215E+02, 0.00000000E+00/ 3972 C DATA (GB( 9,10,IC),IC=1,3) / 3973 C S 0.23776708E+02, 0.10454488E+02, 0.10000000E+01/ 4159 3974 C 4160 3975 C----- INTERVAL = 4 ----- T = 300.0 4161 3976 C 4162 3977 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4163 4164 4165 4166 4167 4168 4169 4170 3978 C DATA (GA(10, 9,IC),IC=1,3) / 3979 C S 0.24000649E+02, 0.10379892E+02, 0.00000000E+00/ 3980 C DATA (GB(10, 9,IC),IC=1,3) / 3981 C S 0.24000649E+02, 0.10497402E+02, 0.10000000E+01/ 3982 C DATA (GA(10,10,IC),IC=1,3) / 3983 C S 0.23714816E+02, 0.10324018E+02, 0.00000000E+00/ 3984 C DATA (GB(10,10,IC),IC=1,3) / 3985 C S 0.23714816E+02, 0.10442501E+02, 0.10000000E+01/ 4171 3986 C 4172 3987 C----- INTERVAL = 4 ----- T = 312.5 4173 3988 C 4174 3989 C-- INDICES FOR PADE APPROXIMATION 1 28 37 45 4175 4176 4177 4178 4179 4180 4181 4182 3990 C DATA (GA(11, 9,IC),IC=1,3) / 3991 C S 0.23943021E+02, 0.10368736E+02, 0.00000000E+00/ 3992 C DATA (GB(11, 9,IC),IC=1,3) / 3993 C S 0.23943021E+02, 0.10486443E+02, 0.10000000E+01/ 3994 C DATA (GA(11,10,IC),IC=1,3) / 3995 C S 0.23658197E+02, 0.10312808E+02, 0.00000000E+00/ 3996 C DATA (GB(11,10,IC),IC=1,3) / 3997 C S 0.23658197E+02, 0.10431483E+02, 0.10000000E+01/ 4183 3998 C 4184 3999 C … … 4193 4008 C 4194 4009 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4195 4196 4197 4198 4199 4200 4201 4202 4010 C DATA (GA( 1, 5,IC),IC=1,3) / 4011 C S 0.15750172E+00,-0.22159303E-01, 0.00000000E+00/ 4012 C DATA (GB( 1, 5,IC),IC=1,3) / 4013 C S 0.15750172E+00, 0.38103212E+00, 0.10000000E+01/ 4014 C DATA (GA( 1, 6,IC),IC=1,3) / 4015 C S 0.17770551E+00,-0.24972399E-01, 0.00000000E+00/ 4016 C DATA (GB( 1, 6,IC),IC=1,3) / 4017 C S 0.17770551E+00, 0.41646579E+00, 0.10000000E+01/ 4203 4018 C 4204 4019 C----- INTERVAL = 5 ----- T = 200.0 4205 4020 C 4206 4021 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4207 4208 4209 4210 4211 4212 4213 4214 4022 C DATA (GA( 2, 5,IC),IC=1,3) / 4023 C S 0.16174076E+00,-0.22748917E-01, 0.00000000E+00/ 4024 C DATA (GB( 2, 5,IC),IC=1,3) / 4025 C S 0.16174076E+00, 0.38913800E+00, 0.10000000E+01/ 4026 C DATA (GA( 2, 6,IC),IC=1,3) / 4027 C S 0.18176757E+00,-0.25537247E-01, 0.00000000E+00/ 4028 C DATA (GB( 2, 6,IC),IC=1,3) / 4029 C S 0.18176757E+00, 0.42345095E+00, 0.10000000E+01/ 4215 4030 C 4216 4031 C----- INTERVAL = 5 ----- T = 212.5 4217 4032 C 4218 4033 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4219 4220 4221 4222 4223 4224 4225 4226 4034 C DATA (GA( 3, 5,IC),IC=1,3) / 4035 C S 0.16548628E+00,-0.23269898E-01, 0.00000000E+00/ 4036 C DATA (GB( 3, 5,IC),IC=1,3) / 4037 C S 0.16548628E+00, 0.39613651E+00, 0.10000000E+01/ 4038 C DATA (GA( 3, 6,IC),IC=1,3) / 4039 C S 0.18527967E+00,-0.26025624E-01, 0.00000000E+00/ 4040 C DATA (GB( 3, 6,IC),IC=1,3) / 4041 C S 0.18527967E+00, 0.42937476E+00, 0.10000000E+01/ 4227 4042 C 4228 4043 C----- INTERVAL = 5 ----- T = 225.0 4229 4044 C 4230 4045 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4231 4232 4233 4234 4235 4236 4237 4238 4046 C DATA (GA( 4, 5,IC),IC=1,3) / 4047 C S 0.16881124E+00,-0.23732392E-01, 0.00000000E+00/ 4048 C DATA (GB( 4, 5,IC),IC=1,3) / 4049 C S 0.16881124E+00, 0.40222421E+00, 0.10000000E+01/ 4050 C DATA (GA( 4, 6,IC),IC=1,3) / 4051 C S 0.18833348E+00,-0.26450280E-01, 0.00000000E+00/ 4052 C DATA (GB( 4, 6,IC),IC=1,3) / 4053 C S 0.18833348E+00, 0.43444062E+00, 0.10000000E+01/ 4239 4054 C 4240 4055 C----- INTERVAL = 5 ----- T = 237.5 4241 4056 C 4242 4057 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4243 4244 4245 4246 4247 4248 4249 4250 4058 C DATA (GA( 5, 5,IC),IC=1,3) / 4059 C S 0.17177839E+00,-0.24145123E-01, 0.00000000E+00/ 4060 C DATA (GB( 5, 5,IC),IC=1,3) / 4061 C S 0.17177839E+00, 0.40756010E+00, 0.10000000E+01/ 4062 C DATA (GA( 5, 6,IC),IC=1,3) / 4063 C S 0.19100108E+00,-0.26821236E-01, 0.00000000E+00/ 4064 C DATA (GB( 5, 6,IC),IC=1,3) / 4065 C S 0.19100108E+00, 0.43880316E+00, 0.10000000E+01/ 4251 4066 C 4252 4067 C----- INTERVAL = 5 ----- T = 250.0 4253 4068 C 4254 4069 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4255 4256 4257 4258 4259 4260 4261 4262 4070 C DATA (GA( 6, 5,IC),IC=1,3) / 4071 C S 0.17443933E+00,-0.24515269E-01, 0.00000000E+00/ 4072 C DATA (GB( 6, 5,IC),IC=1,3) / 4073 C S 0.17443933E+00, 0.41226954E+00, 0.10000000E+01/ 4074 C DATA (GA( 6, 6,IC),IC=1,3) / 4075 C S 0.19334122E+00,-0.27146657E-01, 0.00000000E+00/ 4076 C DATA (GB( 6, 6,IC),IC=1,3) / 4077 C S 0.19334122E+00, 0.44258354E+00, 0.10000000E+01/ 4263 4078 C 4264 4079 C----- INTERVAL = 5 ----- T = 262.5 4265 4080 C 4266 4081 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4267 4268 4269 4270 4271 4272 4273 4274 4082 C DATA (GA( 7, 5,IC),IC=1,3) / 4083 C S 0.17683622E+00,-0.24848690E-01, 0.00000000E+00/ 4084 C DATA (GB( 7, 5,IC),IC=1,3) / 4085 C S 0.17683622E+00, 0.41645142E+00, 0.10000000E+01/ 4086 C DATA (GA( 7, 6,IC),IC=1,3) / 4087 C S 0.19540288E+00,-0.27433354E-01, 0.00000000E+00/ 4088 C DATA (GB( 7, 6,IC),IC=1,3) / 4089 C S 0.19540288E+00, 0.44587882E+00, 0.10000000E+01/ 4275 4090 C 4276 4091 C----- INTERVAL = 5 ----- T = 275.0 4277 4092 C 4278 4093 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4279 4280 4281 4282 4283 4284 4285 4286 4094 C DATA (GA( 8, 5,IC),IC=1,3) / 4095 C S 0.17900375E+00,-0.25150210E-01, 0.00000000E+00/ 4096 C DATA (GB( 8, 5,IC),IC=1,3) / 4097 C S 0.17900375E+00, 0.42018474E+00, 0.10000000E+01/ 4098 C DATA (GA( 8, 6,IC),IC=1,3) / 4099 C S 0.19722732E+00,-0.27687065E-01, 0.00000000E+00/ 4100 C DATA (GB( 8, 6,IC),IC=1,3) / 4101 C S 0.19722732E+00, 0.44876776E+00, 0.10000000E+01/ 4287 4102 C 4288 4103 C----- INTERVAL = 5 ----- T = 287.5 4289 4104 C 4290 4105 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4291 4292 4293 4294 4295 4296 4297 4298 4106 C DATA (GA( 9, 5,IC),IC=1,3) / 4107 C S 0.18097099E+00,-0.25423873E-01, 0.00000000E+00/ 4108 C DATA (GB( 9, 5,IC),IC=1,3) / 4109 C S 0.18097099E+00, 0.42353379E+00, 0.10000000E+01/ 4110 C DATA (GA( 9, 6,IC),IC=1,3) / 4111 C S 0.19884918E+00,-0.27912608E-01, 0.00000000E+00/ 4112 C DATA (GB( 9, 6,IC),IC=1,3) / 4113 C S 0.19884918E+00, 0.45131451E+00, 0.10000000E+01/ 4299 4114 C 4300 4115 C----- INTERVAL = 5 ----- T = 300.0 4301 4116 C 4302 4117 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4303 4304 4305 4306 4307 4308 4309 4310 4118 C DATA (GA(10, 5,IC),IC=1,3) / 4119 C S 0.18276283E+00,-0.25673139E-01, 0.00000000E+00/ 4120 C DATA (GB(10, 5,IC),IC=1,3) / 4121 C S 0.18276283E+00, 0.42655211E+00, 0.10000000E+01/ 4122 C DATA (GA(10, 6,IC),IC=1,3) / 4123 C S 0.20029696E+00,-0.28113944E-01, 0.00000000E+00/ 4124 C DATA (GB(10, 6,IC),IC=1,3) / 4125 C S 0.20029696E+00, 0.45357095E+00, 0.10000000E+01/ 4311 4126 C 4312 4127 C----- INTERVAL = 5 ----- T = 312.5 4313 4128 C 4314 4129 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4315 4316 4317 4318 4319 4320 4321 4322 4130 C DATA (GA(11, 5,IC),IC=1,3) / 4131 C S 0.18440117E+00,-0.25901055E-01, 0.00000000E+00/ 4132 C DATA (GB(11, 5,IC),IC=1,3) / 4133 C S 0.18440117E+00, 0.42928533E+00, 0.10000000E+01/ 4134 C DATA (GA(11, 6,IC),IC=1,3) / 4135 C S 0.20159300E+00,-0.28294180E-01, 0.00000000E+00/ 4136 C DATA (GB(11, 6,IC),IC=1,3) / 4137 C S 0.20159300E+00, 0.45557797E+00, 0.10000000E+01/ 4323 4138 C 4324 4139 C … … 4332 4147 C 4333 4148 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4334 4335 4336 4337 4338 4339 4340 4341 4149 C DATA (GA( 1,11,IC),IC=1,3) / 4150 C S 0.11990218E+02,-0.12823142E+01, 0.00000000E+00/ 4151 C DATA (GB( 1,11,IC),IC=1,3) / 4152 C S 0.11990218E+02, 0.26681588E+02, 0.10000000E+01/ 4153 C DATA (GA( 1,12,IC),IC=1,3) / 4154 C S 0.79709806E+01,-0.74805226E+00, 0.00000000E+00/ 4155 C DATA (GB( 1,12,IC),IC=1,3) / 4156 C S 0.79709806E+01, 0.18377807E+02, 0.10000000E+01/ 4342 4157 C 4343 4158 C----- INTERVAL = 6 ----- T = 200.0 4344 4159 C 4345 4160 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4346 4347 4348 4349 4350 4351 4352 4353 4161 C DATA (GA( 2,11,IC),IC=1,3) / 4162 C S 0.10904073E+02,-0.10571588E+01, 0.00000000E+00/ 4163 C DATA (GB( 2,11,IC),IC=1,3) / 4164 C S 0.10904073E+02, 0.24728346E+02, 0.10000000E+01/ 4165 C DATA (GA( 2,12,IC),IC=1,3) / 4166 C S 0.75400737E+01,-0.56252739E+00, 0.00000000E+00/ 4167 C DATA (GB( 2,12,IC),IC=1,3) / 4168 C S 0.75400737E+01, 0.17643148E+02, 0.10000000E+01/ 4354 4169 C 4355 4170 C----- INTERVAL = 6 ----- T = 212.5 4356 4171 C 4357 4172 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4358 4359 4360 4361 4362 4363 4364 4365 4173 C DATA (GA( 3,11,IC),IC=1,3) / 4174 C S 0.89126838E+01,-0.74864953E+00, 0.00000000E+00/ 4175 C DATA (GB( 3,11,IC),IC=1,3) / 4176 C S 0.89126838E+01, 0.20551342E+02, 0.10000000E+01/ 4177 C DATA (GA( 3,12,IC),IC=1,3) / 4178 C S 0.81804377E+01,-0.46188072E+00, 0.00000000E+00/ 4179 C DATA (GB( 3,12,IC),IC=1,3) / 4180 C S 0.81804377E+01, 0.19296161E+02, 0.10000000E+01/ 4366 4181 C 4367 4182 C----- INTERVAL = 6 ----- T = 225.0 4368 4183 C 4369 4184 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4370 4371 4372 4373 4374 4375 4376 4377 4185 C DATA (GA( 4,11,IC),IC=1,3) / 4186 C S 0.85622405E+01,-0.58705980E+00, 0.00000000E+00/ 4187 C DATA (GB( 4,11,IC),IC=1,3) / 4188 C S 0.85622405E+01, 0.19955244E+02, 0.10000000E+01/ 4189 C DATA (GA( 4,12,IC),IC=1,3) / 4190 C S 0.10564339E+02,-0.40712065E+00, 0.00000000E+00/ 4191 C DATA (GB( 4,12,IC),IC=1,3) / 4192 C S 0.10564339E+02, 0.24951120E+02, 0.10000000E+01/ 4378 4193 C 4379 4194 C----- INTERVAL = 6 ----- T = 237.5 4380 4195 C 4381 4196 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4382 4383 4384 4385 4386 4387 4388 4389 4197 C DATA (GA( 5,11,IC),IC=1,3) / 4198 C S 0.94892164E+01,-0.49305772E+00, 0.00000000E+00/ 4199 C DATA (GB( 5,11,IC),IC=1,3) / 4200 C S 0.94892164E+01, 0.22227100E+02, 0.10000000E+01/ 4201 C DATA (GA( 5,12,IC),IC=1,3) / 4202 C S 0.46896789E+02,-0.15295996E+01, 0.00000000E+00/ 4203 C DATA (GB( 5,12,IC),IC=1,3) / 4204 C S 0.46896789E+02, 0.10957372E+03, 0.10000000E+01/ 4390 4205 C 4391 4206 C----- INTERVAL = 6 ----- T = 250.0 4392 4207 C 4393 4208 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4394 4395 4396 4397 4398 4399 4400 4401 4209 C DATA (GA( 6,11,IC),IC=1,3) / 4210 C S 0.13580937E+02,-0.51461431E+00, 0.00000000E+00/ 4211 C DATA (GB( 6,11,IC),IC=1,3) / 4212 C S 0.13580937E+02, 0.31770288E+02, 0.10000000E+01/ 4213 C DATA (GA( 6,12,IC),IC=1,3) / 4214 C S-0.30926524E+01, 0.43555255E+00, 0.00000000E+00/ 4215 C DATA (GB( 6,12,IC),IC=1,3) / 4216 C S-0.30926524E+01,-0.67432659E+01, 0.10000000E+01/ 4402 4217 C 4403 4218 C----- INTERVAL = 6 ----- T = 262.5 4404 4219 C 4405 4220 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4406 4407 4408 4409 4410 4411 4412 4413 4221 C DATA (GA( 7,11,IC),IC=1,3) / 4222 C S-0.32050918E+03, 0.12373350E+02, 0.00000000E+00/ 4223 C DATA (GB( 7,11,IC),IC=1,3) / 4224 C S-0.32050918E+03,-0.74061287E+03, 0.10000000E+01/ 4225 C DATA (GA( 7,12,IC),IC=1,3) / 4226 C S 0.85742941E+00, 0.50380874E+00, 0.00000000E+00/ 4227 C DATA (GB( 7,12,IC),IC=1,3) / 4228 C S 0.85742941E+00, 0.24550746E+01, 0.10000000E+01/ 4414 4229 C 4415 4230 C----- INTERVAL = 6 ----- T = 275.0 4416 4231 C 4417 4232 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4418 4419 4420 4421 4422 4423 4424 4425 4233 C DATA (GA( 8,11,IC),IC=1,3) / 4234 C S-0.37133165E+01, 0.44809588E+00, 0.00000000E+00/ 4235 C DATA (GB( 8,11,IC),IC=1,3) / 4236 C S-0.37133165E+01,-0.81329826E+01, 0.10000000E+01/ 4237 C DATA (GA( 8,12,IC),IC=1,3) / 4238 C S 0.19164038E+01, 0.68537352E+00, 0.00000000E+00/ 4239 C DATA (GB( 8,12,IC),IC=1,3) / 4240 C S 0.19164038E+01, 0.49089917E+01, 0.10000000E+01/ 4426 4241 C 4427 4242 C----- INTERVAL = 6 ----- T = 287.5 4428 4243 C 4429 4244 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4430 4431 4432 4433 4434 4435 4436 4437 4245 C DATA (GA( 9,11,IC),IC=1,3) / 4246 C S 0.18890836E+00, 0.46548918E+00, 0.00000000E+00/ 4247 C DATA (GB( 9,11,IC),IC=1,3) / 4248 C S 0.18890836E+00, 0.90279822E+00, 0.10000000E+01/ 4249 C DATA (GA( 9,12,IC),IC=1,3) / 4250 C S 0.23513199E+01, 0.89437630E+00, 0.00000000E+00/ 4251 C DATA (GB( 9,12,IC),IC=1,3) / 4252 C S 0.23513199E+01, 0.59008712E+01, 0.10000000E+01/ 4438 4253 C 4439 4254 C----- INTERVAL = 6 ----- T = 300.0 4440 4255 C 4441 4256 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4442 4443 4444 4445 4446 4447 4448 4449 4257 C DATA (GA(10,11,IC),IC=1,3) / 4258 C S 0.14209226E+01, 0.59121475E+00, 0.00000000E+00/ 4259 C DATA (GB(10,11,IC),IC=1,3) / 4260 C S 0.14209226E+01, 0.37532746E+01, 0.10000000E+01/ 4261 C DATA (GA(10,12,IC),IC=1,3) / 4262 C S 0.25566644E+01, 0.11127003E+01, 0.00000000E+00/ 4263 C DATA (GB(10,12,IC),IC=1,3) / 4264 C S 0.25566644E+01, 0.63532616E+01, 0.10000000E+01/ 4450 4265 C 4451 4266 C----- INTERVAL = 6 ----- T = 312.5 4452 4267 C 4453 4268 C-- INDICES FOR PADE APPROXIMATION 1 35 40 45 4454 4455 4456 4457 4458 4459 4460 4461 4269 C DATA (GA(11,11,IC),IC=1,3) / 4270 C S 0.19817679E+01, 0.74676119E+00, 0.00000000E+00/ 4271 C DATA (GB(11,11,IC),IC=1,3) / 4272 C S 0.19817679E+01, 0.50437916E+01, 0.10000000E+01/ 4273 C DATA (GA(11,12,IC),IC=1,3) / 4274 C S 0.26555181E+01, 0.13329782E+01, 0.00000000E+00/ 4275 C DATA (GB(11,12,IC),IC=1,3) / 4276 C S 0.26555181E+01, 0.65558627E+01, 0.10000000E+01/ 4462 4277 C 4463 4278 C … … 4477 4292 C 4478 4293 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4479 4480 4481 4482 4483 4484 4485 4486 4294 C DATA (GA( 1,13,IC),IC=1,3) / 4295 C S 0.87668459E-01, 0.13845511E+01, 0.00000000E+00/ 4296 C DATA (GB( 1,13,IC),IC=1,3) / 4297 C S 0.87668459E-01, 0.23203798E+01, 0.10000000E+01/ 4298 C DATA (GA( 1,14,IC),IC=1,3) / 4299 C S 0.74878820E-01, 0.11718758E+01, 0.00000000E+00/ 4300 C DATA (GB( 1,14,IC),IC=1,3) / 4301 C S 0.74878820E-01, 0.20206726E+01, 0.10000000E+01/ 4487 4302 C 4488 4303 C----- INTERVAL = 2 ----- T = 200.0 4489 4304 C 4490 4305 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4491 4492 4493 4494 4495 4496 4497 4498 4306 C DATA (GA( 2,13,IC),IC=1,3) / 4307 C S 0.83754276E-01, 0.13187042E+01, 0.00000000E+00/ 4308 C DATA (GB( 2,13,IC),IC=1,3) / 4309 C S 0.83754276E-01, 0.22288925E+01, 0.10000000E+01/ 4310 C DATA (GA( 2,14,IC),IC=1,3) / 4311 C S 0.71650966E-01, 0.11216131E+01, 0.00000000E+00/ 4312 C DATA (GB( 2,14,IC),IC=1,3) / 4313 C S 0.71650966E-01, 0.19441824E+01, 0.10000000E+01/ 4499 4314 C 4500 4315 C----- INTERVAL = 2 ----- T = 212.5 4501 4316 C 4502 4317 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4503 4504 4505 4506 4507 4508 4509 4510 4318 C DATA (GA( 3,13,IC),IC=1,3) / 4319 C S 0.80460283E-01, 0.12644396E+01, 0.00000000E+00/ 4320 C DATA (GB( 3,13,IC),IC=1,3) / 4321 C S 0.80460283E-01, 0.21515593E+01, 0.10000000E+01/ 4322 C DATA (GA( 3,14,IC),IC=1,3) / 4323 C S 0.68979615E-01, 0.10809473E+01, 0.00000000E+00/ 4324 C DATA (GB( 3,14,IC),IC=1,3) / 4325 C S 0.68979615E-01, 0.18807257E+01, 0.10000000E+01/ 4511 4326 C 4512 4327 C----- INTERVAL = 2 ----- T = 225.0 4513 4328 C 4514 4329 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4515 4516 4517 4518 4519 4520 4521 4522 4330 C DATA (GA( 4,13,IC),IC=1,3) / 4331 C S 0.77659686E-01, 0.12191543E+01, 0.00000000E+00/ 4332 C DATA (GB( 4,13,IC),IC=1,3) / 4333 C S 0.77659686E-01, 0.20855896E+01, 0.10000000E+01/ 4334 C DATA (GA( 4,14,IC),IC=1,3) / 4335 C S 0.66745345E-01, 0.10476396E+01, 0.00000000E+00/ 4336 C DATA (GB( 4,14,IC),IC=1,3) / 4337 C S 0.66745345E-01, 0.18275618E+01, 0.10000000E+01/ 4523 4338 C 4524 4339 C----- INTERVAL = 2 ----- T = 237.5 4525 4340 C 4526 4341 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4527 4528 4529 4530 4531 4532 4533 4534 4342 C DATA (GA( 5,13,IC),IC=1,3) / 4343 C S 0.75257056E-01, 0.11809511E+01, 0.00000000E+00/ 4344 C DATA (GB( 5,13,IC),IC=1,3) / 4345 C S 0.75257056E-01, 0.20288489E+01, 0.10000000E+01/ 4346 C DATA (GA( 5,14,IC),IC=1,3) / 4347 C S 0.64857571E-01, 0.10200373E+01, 0.00000000E+00/ 4348 C DATA (GB( 5,14,IC),IC=1,3) / 4349 C S 0.64857571E-01, 0.17825910E+01, 0.10000000E+01/ 4535 4350 C 4536 4351 C----- INTERVAL = 2 ----- T = 250.0 4537 4352 C 4538 4353 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4539 4540 4541 4542 4543 4544 4545 4546 4354 C DATA (GA( 6,13,IC),IC=1,3) / 4355 C S 0.73179175E-01, 0.11484154E+01, 0.00000000E+00/ 4356 C DATA (GB( 6,13,IC),IC=1,3) / 4357 C S 0.73179175E-01, 0.19796791E+01, 0.10000000E+01/ 4358 C DATA (GA( 6,14,IC),IC=1,3) / 4359 C S 0.63248495E-01, 0.99692726E+00, 0.00000000E+00/ 4360 C DATA (GB( 6,14,IC),IC=1,3) / 4361 C S 0.63248495E-01, 0.17442308E+01, 0.10000000E+01/ 4547 4362 C 4548 4363 C----- INTERVAL = 2 ----- T = 262.5 4549 4364 C 4550 4365 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4551 4552 4553 4554 4555 4556 4557 4558 4366 C DATA (GA( 7,13,IC),IC=1,3) / 4367 C S 0.71369063E-01, 0.11204723E+01, 0.00000000E+00/ 4368 C DATA (GB( 7,13,IC),IC=1,3) / 4369 C S 0.71369063E-01, 0.19367778E+01, 0.10000000E+01/ 4370 C DATA (GA( 7,14,IC),IC=1,3) / 4371 C S 0.61866970E-01, 0.97740923E+00, 0.00000000E+00/ 4372 C DATA (GB( 7,14,IC),IC=1,3) / 4373 C S 0.61866970E-01, 0.17112809E+01, 0.10000000E+01/ 4559 4374 C 4560 4375 C----- INTERVAL = 2 ----- T = 275.0 4561 4376 C 4562 4377 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4563 4564 4565 4566 4567 4568 4569 4570 4378 C DATA (GA( 8,13,IC),IC=1,3) / 4379 C S 0.69781812E-01, 0.10962918E+01, 0.00000000E+00/ 4380 C DATA (GB( 8,13,IC),IC=1,3) / 4381 C S 0.69781812E-01, 0.18991112E+01, 0.10000000E+01/ 4382 C DATA (GA( 8,14,IC),IC=1,3) / 4383 C S 0.60673632E-01, 0.96080188E+00, 0.00000000E+00/ 4384 C DATA (GB( 8,14,IC),IC=1,3) / 4385 C S 0.60673632E-01, 0.16828137E+01, 0.10000000E+01/ 4571 4386 C 4572 4387 C----- INTERVAL = 2 ----- T = 287.5 4573 4388 C 4574 4389 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4575 4576 4577 4578 4579 4580 4581 4582 4390 C DATA (GA( 9,13,IC),IC=1,3) / 4391 C S 0.68381606E-01, 0.10752229E+01, 0.00000000E+00/ 4392 C DATA (GB( 9,13,IC),IC=1,3) / 4393 C S 0.68381606E-01, 0.18658501E+01, 0.10000000E+01/ 4394 C DATA (GA( 9,14,IC),IC=1,3) / 4395 C S 0.59637277E-01, 0.94657562E+00, 0.00000000E+00/ 4396 C DATA (GB( 9,14,IC),IC=1,3) / 4397 C S 0.59637277E-01, 0.16580908E+01, 0.10000000E+01/ 4583 4398 C 4584 4399 C----- INTERVAL = 2 ----- T = 300.0 4585 4400 C 4586 4401 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4587 4588 4589 4590 4591 4592 4593 4594 4402 C DATA (GA(10,13,IC),IC=1,3) / 4403 C S 0.67139539E-01, 0.10567474E+01, 0.00000000E+00/ 4404 C DATA (GB(10,13,IC),IC=1,3) / 4405 C S 0.67139539E-01, 0.18363226E+01, 0.10000000E+01/ 4406 C DATA (GA(10,14,IC),IC=1,3) / 4407 C S 0.58732178E-01, 0.93430511E+00, 0.00000000E+00/ 4408 C DATA (GB(10,14,IC),IC=1,3) / 4409 C S 0.58732178E-01, 0.16365014E+01, 0.10000000E+01/ 4595 4410 C 4596 4411 C----- INTERVAL = 2 ----- T = 312.5 4597 4412 C 4598 4413 C-- INDICES FOR PADE APPROXIMATION 1 30 38 45 4599 4600 4601 4602 4603 4604 4605 4606 4414 C DATA (GA(11,13,IC),IC=1,3) / 4415 C S 0.66032012E-01, 0.10404465E+01, 0.00000000E+00/ 4416 C DATA (GB(11,13,IC),IC=1,3) / 4417 C S 0.66032012E-01, 0.18099779E+01, 0.10000000E+01/ 4418 C DATA (GA(11,14,IC),IC=1,3) / 4419 C S 0.57936092E-01, 0.92363528E+00, 0.00000000E+00/ 4420 C DATA (GB(11,14,IC),IC=1,3) / 4421 C S 0.57936092E-01, 0.16175164E+01, 0.10000000E+01/ 4607 4422 C 4608 4423 C … … 4624 4439 C 4625 4440 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4626 4627 4628 4629 4630 4631 4632 4633 4441 C DATA (GA( 1,15,IC),IC=1,3) / 4442 C S 0.13230067E+02, 0.22042132E+02, 0.00000000E+00/ 4443 C DATA (GB( 1,15,IC),IC=1,3) / 4444 C S 0.13230067E+02, 0.22051750E+02, 0.10000000E+01/ 4445 C DATA (GA( 1,16,IC),IC=1,3) / 4446 C S 0.13183816E+02, 0.22169501E+02, 0.00000000E+00/ 4447 C DATA (GB( 1,16,IC),IC=1,3) / 4448 C S 0.13183816E+02, 0.22178972E+02, 0.10000000E+01/ 4634 4449 C 4635 4450 C----- INTERVAL = 4 ----- T = 200.0 4636 4451 C 4637 4452 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4638 4639 4640 4641 4642 4643 4644 4645 4453 C DATA (GA( 2,15,IC),IC=1,3) / 4454 C S 0.13213564E+02, 0.22107298E+02, 0.00000000E+00/ 4455 C DATA (GB( 2,15,IC),IC=1,3) / 4456 C S 0.13213564E+02, 0.22116850E+02, 0.10000000E+01/ 4457 C DATA (GA( 2,16,IC),IC=1,3) / 4458 C S 0.13189991E+02, 0.22270075E+02, 0.00000000E+00/ 4459 C DATA (GB( 2,16,IC),IC=1,3) / 4460 C S 0.13189991E+02, 0.22279484E+02, 0.10000000E+01/ 4646 4461 C 4647 4462 C----- INTERVAL = 4 ----- T = 212.5 4648 4463 C 4649 4464 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4650 4651 4652 4653 4654 4655 4656 4657 4465 C DATA (GA( 3,15,IC),IC=1,3) / 4466 C S 0.13209140E+02, 0.22180915E+02, 0.00000000E+00/ 4467 C DATA (GB( 3,15,IC),IC=1,3) / 4468 C S 0.13209140E+02, 0.22190410E+02, 0.10000000E+01/ 4469 C DATA (GA( 3,16,IC),IC=1,3) / 4470 C S 0.13209485E+02, 0.22379193E+02, 0.00000000E+00/ 4471 C DATA (GB( 3,16,IC),IC=1,3) / 4472 C S 0.13209485E+02, 0.22388551E+02, 0.10000000E+01/ 4658 4473 C 4659 4474 C----- INTERVAL = 4 ----- T = 225.0 4660 4475 C 4661 4476 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4662 4663 4664 4665 4666 4667 4668 4669 4477 C DATA (GA( 4,15,IC),IC=1,3) / 4478 C S 0.13213894E+02, 0.22259478E+02, 0.00000000E+00/ 4479 C DATA (GB( 4,15,IC),IC=1,3) / 4480 C S 0.13213894E+02, 0.22268925E+02, 0.10000000E+01/ 4481 C DATA (GA( 4,16,IC),IC=1,3) / 4482 C S 0.13238789E+02, 0.22492992E+02, 0.00000000E+00/ 4483 C DATA (GB( 4,16,IC),IC=1,3) / 4484 C S 0.13238789E+02, 0.22502309E+02, 0.10000000E+01/ 4670 4485 C 4671 4486 C----- INTERVAL = 4 ----- T = 237.5 4672 4487 C 4673 4488 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4674 4675 4676 4677 4678 4679 4680 4681 4489 C DATA (GA( 5,15,IC),IC=1,3) / 4490 C S 0.13225963E+02, 0.22341039E+02, 0.00000000E+00/ 4491 C DATA (GB( 5,15,IC),IC=1,3) / 4492 C S 0.13225963E+02, 0.22350445E+02, 0.10000000E+01/ 4493 C DATA (GA( 5,16,IC),IC=1,3) / 4494 C S 0.13275017E+02, 0.22608508E+02, 0.00000000E+00/ 4495 C DATA (GB( 5,16,IC),IC=1,3) / 4496 C S 0.13275017E+02, 0.22617792E+02, 0.10000000E+01/ 4682 4497 C 4683 4498 C----- INTERVAL = 4 ----- T = 250.0 4684 4499 C 4685 4500 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4686 4687 4688 4689 4690 4691 4692 4693 4501 C DATA (GA( 6,15,IC),IC=1,3) / 4502 C S 0.13243806E+02, 0.22424247E+02, 0.00000000E+00/ 4503 C DATA (GB( 6,15,IC),IC=1,3) / 4504 C S 0.13243806E+02, 0.22433617E+02, 0.10000000E+01/ 4505 C DATA (GA( 6,16,IC),IC=1,3) / 4506 C S 0.13316096E+02, 0.22723843E+02, 0.00000000E+00/ 4507 C DATA (GB( 6,16,IC),IC=1,3) / 4508 C S 0.13316096E+02, 0.22733099E+02, 0.10000000E+01/ 4694 4509 C 4695 4510 C----- INTERVAL = 4 ----- T = 262.5 4696 4511 C 4697 4512 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4698 4699 4700 4701 4702 4703 4704 4705 4513 C DATA (GA( 7,15,IC),IC=1,3) / 4514 C S 0.13266104E+02, 0.22508089E+02, 0.00000000E+00/ 4515 C DATA (GB( 7,15,IC),IC=1,3) / 4516 C S 0.13266104E+02, 0.22517429E+02, 0.10000000E+01/ 4517 C DATA (GA( 7,16,IC),IC=1,3) / 4518 C S 0.13360555E+02, 0.22837837E+02, 0.00000000E+00/ 4519 C DATA (GB( 7,16,IC),IC=1,3) / 4520 C S 0.13360555E+02, 0.22847071E+02, 0.10000000E+01/ 4706 4521 C 4707 4522 C----- INTERVAL = 4 ----- T = 275.0 4708 4523 C 4709 4524 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4710 4711 4712 4713 4714 4715 4716 4717 4525 C DATA (GA( 8,15,IC),IC=1,3) / 4526 C S 0.13291782E+02, 0.22591771E+02, 0.00000000E+00/ 4527 C DATA (GB( 8,15,IC),IC=1,3) / 4528 C S 0.13291782E+02, 0.22601086E+02, 0.10000000E+01/ 4529 C DATA (GA( 8,16,IC),IC=1,3) / 4530 C S 0.13407324E+02, 0.22949751E+02, 0.00000000E+00/ 4531 C DATA (GB( 8,16,IC),IC=1,3) / 4532 C S 0.13407324E+02, 0.22958967E+02, 0.10000000E+01/ 4718 4533 C 4719 4534 C----- INTERVAL = 4 ----- T = 287.5 4720 4535 C 4721 4536 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4722 4723 4724 4725 4726 4727 4728 4729 4537 C DATA (GA( 9,15,IC),IC=1,3) / 4538 C S 0.13319961E+02, 0.22674661E+02, 0.00000000E+00/ 4539 C DATA (GB( 9,15,IC),IC=1,3) / 4540 C S 0.13319961E+02, 0.22683956E+02, 0.10000000E+01/ 4541 C DATA (GA( 9,16,IC),IC=1,3) / 4542 C S 0.13455544E+02, 0.23059032E+02, 0.00000000E+00/ 4543 C DATA (GB( 9,16,IC),IC=1,3) / 4544 C S 0.13455544E+02, 0.23068234E+02, 0.10000000E+01/ 4730 4545 C 4731 4546 C----- INTERVAL = 4 ----- T = 300.0 4732 4547 C 4733 4548 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4734 4735 4736 4737 4738 4739 4740 4741 4549 C DATA (GA(10,15,IC),IC=1,3) / 4550 C S 0.13349927E+02, 0.22756246E+02, 0.00000000E+00/ 4551 C DATA (GB(10,15,IC),IC=1,3) / 4552 C S 0.13349927E+02, 0.22765522E+02, 0.10000000E+01/ 4553 C DATA (GA(10,16,IC),IC=1,3) / 4554 C S 0.13504450E+02, 0.23165146E+02, 0.00000000E+00/ 4555 C DATA (GB(10,16,IC),IC=1,3) / 4556 C S 0.13504450E+02, 0.23174336E+02, 0.10000000E+01/ 4742 4557 C 4743 4558 C----- INTERVAL = 4 ----- T = 312.5 4744 4559 C 4745 4560 C-- INDICES FOR PADE APPROXIMATION 1 15 29 45 4746 DATA (GA(11,15,IC),IC=1,3) / 4747 S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/ 4748 DATA (GB(11,15,IC),IC=1,3) / 4749 S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/ 4750 DATA (GA(11,16,IC),IC=1,3) / 4751 S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/ 4752 DATA (GB(11,16,IC),IC=1,3) / 4753 S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/ 4561 C DATA (GA(11,15,IC),IC=1,3) / 4562 C S 0.13381108E+02, 0.22836093E+02, 0.00000000E+00/ 4563 C DATA (GB(11,15,IC),IC=1,3) / 4564 C S 0.13381108E+02, 0.22845354E+02, 0.10000000E+01/ 4565 C DATA (GA(11,16,IC),IC=1,3) / 4566 C S 0.13553282E+02, 0.23267456E+02, 0.00000000E+00/ 4567 C DATA (GB(11,16,IC),IC=1,3) / 4568 C S 0.13553282E+02, 0.23276638E+02, 0.10000000E+01/ 4569 C 4570 C ------------------------------------------------------------------ 4571 C DATA (( XP( J,K),J=1,6), K=1,6) / 4572 C S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03, 4573 C S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03, 4574 C S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03, 4575 C S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02, 4576 C S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03, 4577 C S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02, 4578 C S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03, 4579 C S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02, 4580 C S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02, 4581 C S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01, 4582 C S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03, 4583 C S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 / 4754 4584 4755 C ------------------------------------------------------------------4756 DATA (( XP( J,K),J=1,6), K=1,6) /4757 S 0.46430621E+02, 0.12928299E+03, 0.20732648E+03,4758 S 0.31398411E+03, 0.18373177E+03,-0.11412303E+03,4759 S 0.73604774E+02, 0.27887914E+03, 0.27076947E+03,4760 S-0.57322111E+02,-0.64742459E+02, 0.87238280E+02,4761 S 0.37050866E+02, 0.20498759E+03, 0.37558029E+03,4762 S 0.17401171E+03,-0.13350302E+03,-0.37651795E+02,4763 S 0.14930141E+02, 0.89161160E+02, 0.17793062E+03,4764 S 0.93433860E+02,-0.70646020E+02,-0.26373150E+02,4765 S 0.40386780E+02, 0.10855270E+03, 0.50755010E+02,4766 S-0.31496190E+02, 0.12791300E+00, 0.18017770E+01,4767 S 0.90811926E+01, 0.75073923E+02, 0.24654438E+03,4768 S 0.39332612E+03, 0.29385281E+03, 0.89107921E+02 /4769 4585 C 4770 4586 C … … 4774 4590 100 CONTINUE 4775 4591 C 4592 !cdir collapse 4776 4593 DO 102 JK = 1 , KFLEV+1 4777 4594 DO 101 JL = 1, KDLON … … 5625 5442 S , PADJD,PADJU,PCNTRB,PDBDT) 5626 5443 USE dimphy 5444 USE radiation_AR4_param, only : WG1 5627 5445 IMPLICIT none 5628 5446 cym#include "dimensions.h" … … 5685 5503 REAL(KIND=8) zwtr 5686 5504 c 5687 C* Data Block: 5688 c 5689 REAL(KIND=8) WG1(2) 5690 SAVE WG1 5691 c$OMP THREADPRIVATE(WG1) 5692 DATA (WG1(jk),jk=1,2) /1.0, 1.0/ 5505 5693 5506 C----------------------------------------------------------------------- 5694 5507 C … … 5903 5716 C 5904 5717 C 5718 !cdir collapse 5905 5719 DO 130 JA = 1 , 8 5906 5720 DO 120 JL = 1, KDLON … … 6080 5894 C 6081 5895 C 5896 5897 !CDIR ON_ADB(PUU1) 5898 !CDIR ON_ADB(PUU2) 5899 !CDIR COLLAPSE 6082 5900 DO 130 JA = 1 , 8 6083 5901 DO 120 JL = 1, KDLON -
LMDZ4/branches/LMDZ4-dev/makelmdz_fcm
r1249 r1250 162 162 # lecture des chemins propres à l'architecture de la machine # 163 163 ############################################################### 164 165 rm -f ./arch.path 166 ln -s ./arch/arch-${arch}.path ./arch.path 164 rm -f .void_file 165 echo > .void_file 166 rm -f arch.path 167 ln -s arch/arch-${arch}.path ./arch.path 167 168 source arch.path 168 169 … … 360 361 rm -f bin/${code}${SUFF_NAME}.e 361 362 rm -f arch.fcm 363 rm -f arch.opt 362 364 363 365 echo "%ARCH $arch" >> $config_fcm … … 380 382 381 383 ln -s arch/arch-${arch}.fcm arch.fcm 384 if test -f arch/arch-${arch}.opt && [ $compil_mod = "prod" ] 385 then 386 ln -s arch/arch-${arch}.opt arch.opt 387 else 388 ln -s .void_file arch.opt 389 fi 390 391 382 392 rm -f $LIBOGCM/${arch}${SUFF_NAME}/.config/fcm.bld.lock 383 393 ./build_gcm … … 387 397 ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config config 388 398 ln -s $LIBOGCM/${arch}${SUFF_NAME}/.config/tmp tmp_src 389 #\rm -f $libf/grid/dimensions.h
Note: See TracChangeset
for help on using the changeset viewer.
