Changeset 5086 for LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm
- Timestamp:
- Jul 19, 2024, 7:54:50 PM (5 months ago)
- Location:
- LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/dump2ds.F
r5082 r5086 27 27 REAL zmin,zmax,zllu,zllm 28 28 write(F1000,'(''(4x,'',I3,''(1H-))'')')im+3 29 DO 10001i=1,20029 DO i=1,200 30 30 jline(1+(i-1)*5:5*i)='. ' 31 10001 CONTINUE 31 END DO 32 32 10002 zmin=z(1,1) 33 33 imin=1 … … 37 37 jmax=1 38 38 kzero=0 39 DO 10003j=1,jm40 DO 10005i=1,im39 DO j=1,jm 40 DO i=1,im 41 41 IF(.NOT.( z(i,j)>zmax))GOTO 10007 42 42 zmax=z(i,j) … … 50 50 kzero=kzero+1 51 51 10011 CONTINUE 52 10005 CONTINUE 52 END DO 53 53 10006 CONTINUE 54 10003 CONTINUE 54 END DO 55 55 10004 zsign=(sign(1.,zmin)*sign(1.,zmax)>0.) 56 56 WRITE(*,*)'>>> dump2ds: ',trim(nom_z) … … 73 73 zinf=.false. 74 74 znan=.false. 75 DO 10017j=1,jm76 DO 10019i=1,im75 DO j=1,jm 76 DO i=1,im 77 77 az=abs(z(i,j)) 78 78 IF(.NOT.( az==0.))GOTO 10021 … … 97 97 kchar(i)=32-kchar(i) 98 98 10027 CONTINUE 99 10019 CONTINUE 99 END DO 100 100 10020 WRITE(*,'(1x,i3,''|'',1000a)')j+1,(iform(kchar(i)),i=1,im),'|','| 101 101 *' 102 10017 CONTINUE 102 END DO 103 103 10018 write(*,F1000) 104 104 WRITE(*,'(5x,1000i1)')(mod(i,10),i=1,im) -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/eq_regions_mod.F90
r5082 r5086 353 353 do k = 2, n 354 354 w = w * (x - k) 355 end do355 END DO 356 356 else 357 357 w = 1 358 358 do k = 0, -n - 1 359 359 y = y * (x + k) 360 end do360 END DO 361 361 end if 362 362 gamma_res = w / y -
LMDZ6/branches/Amaury_dev/libf/phylmd/rrtm/set99.F
r5082 r5086 15 15 NIL=0 16 16 NHL=(N/2)-1 17 DO 10K=NIL,NHL17 DO K=NIL,NHL 18 18 ANGLE=FLOAT(K)*DEL 19 19 TRIGS(2*K+1)=COS(ANGLE) 20 20 TRIGS(2*K+2)=SIN(ANGLE) 21 10 CONTINUE21 END DO 22 22 C 23 23 C FIND FACTORS OF N (8,6,5,4,3,2; ONLY ONE 8 ALLOWED)
Note: See TracChangeset
for help on using the changeset viewer.