[3331] | 1 | MODULE cfmip_point_locations |
---|
| 2 | IMPLICIT NONE |
---|
| 3 | |
---|
| 4 | CONTAINS |
---|
| 5 | |
---|
| 6 | SUBROUTINE read_CFMIP_point_locations(npCFMIP, tab, lonCFMIP, latCFMIP) |
---|
| 7 | IMPLICIT none |
---|
| 8 | INTEGER :: npCFMIP |
---|
| 9 | REAL, DIMENSION(npCFMIP) :: lonCFMIP, latCFMIP |
---|
| 10 | INTEGER :: tab(npCFMIP), np |
---|
| 11 | |
---|
| 12 | WRITE(*,*) 'npCFMIP=',npCFMIP |
---|
| 13 | ! OPEN(20, file="pointlocations.txt",status='old') |
---|
| 14 | OPEN(20, file="pointlocations.txt",status='old',err=999) |
---|
| 15 | OPEN(21, file="pointlocations_lon180.txt") |
---|
| 16 | np=1 |
---|
| 17 | 10 READ(20,*) tab(np), lonCFMIP(np), latCFMIP(np) |
---|
| 18 | !!! passage de 0-360 a -180/180 |
---|
| 19 | IF (lonCFMIP(np).GT.180.) THEN |
---|
| 20 | lonCFMIP(np)=lonCFMIP(np)-360. |
---|
| 21 | ENDIF |
---|
| 22 | WRITE(21,*) np, lonCFMIP(np), latCFMIP(np) |
---|
| 23 | np=np+1 |
---|
| 24 | IF(np.LE.npCFMIP) THEN |
---|
| 25 | GOTO 10 |
---|
| 26 | ENDIF |
---|
| 27 | CLOSE(20) |
---|
| 28 | CLOSE(21) |
---|
| 29 | 999 RETURN |
---|
| 30 | END SUBROUTINE read_CFMIP_point_locations |
---|
| 31 | |
---|
| 32 | SUBROUTINE LMDZ_CFMIP_point_locations(npCFMIP, lonCFMIP, latCFMIP, & |
---|
| 33 | tabijGCM, lonGCM, latGCM, ipt, jpt) |
---|
| 34 | USE dimphy |
---|
| 35 | USE iophy |
---|
| 36 | USE mod_grid_phy_lmdz, ONLY: nbp_lon, nbp_lat, klon_glo |
---|
| 37 | |
---|
| 38 | IMPLICIT none |
---|
| 39 | INTEGER :: npCFMIP |
---|
| 40 | REAL, DIMENSION(npCFMIP) :: lonCFMIP, latCFMIP |
---|
| 41 | INTEGER :: i, j, np, ip |
---|
| 42 | INTEGER, DIMENSION(npCFMIP) :: ipt, jpt |
---|
| 43 | REAL :: dlon1, dlon2 |
---|
| 44 | REAL :: dlat1, dlat2 |
---|
| 45 | REAL, DIMENSION(nbp_lon+1) :: lon |
---|
| 46 | INTEGER, DIMENSION(npCFMIP) :: tabijGCM |
---|
| 47 | REAL, DIMENSION(npCFMIP) :: lonGCM, latGCM |
---|
| 48 | |
---|
| 49 | lon(1:nbp_lon)=io_lon(:) |
---|
| 50 | lon(nbp_lon+1)=-1*lon(1) |
---|
| 51 | OPEN(22, file="LMDZ_pointsCFMIP.txt") |
---|
| 52 | DO np=1, npCFMIP |
---|
| 53 | DO i=1, nbp_lon |
---|
| 54 | ! |
---|
| 55 | ! PRINT*,'IM np i lonCF lonGCM lonGCM+1',np,i,lonCFMIP(np),lon(i), & |
---|
| 56 | ! lon(i+1) |
---|
| 57 | ! |
---|
| 58 | IF(lonCFMIP(np).GE.lon(i).AND.lonCFMIP(np).LT.lon(i+1)) THEN |
---|
| 59 | dlon1 = abs (lonCFMIP(np) - lon(i)) |
---|
| 60 | dlon2 = abs (lonCFMIP(np) - lon(i+1)) |
---|
| 61 | IF (dlon1.LE.dlon2) THEN |
---|
| 62 | ipt(np)=i |
---|
| 63 | ELSE |
---|
| 64 | ipt(np)=i+1 |
---|
| 65 | ENDIF |
---|
| 66 | ENDIF |
---|
| 67 | END DO |
---|
| 68 | END DO |
---|
| 69 | ! |
---|
| 70 | np=1 |
---|
| 71 | 30 j=1 |
---|
| 72 | 40 IF(latCFMIP(np).LE.io_lat(j).AND.latCFMIP(np).GE.io_lat(j+1)) THEN |
---|
| 73 | dlat1 = abs (latCFMIP(np) - io_lat(j)) |
---|
| 74 | dlat2 = abs (latCFMIP(np) - io_lat(j+1)) |
---|
| 75 | IF (dlat1.LE.dlat2) THEN |
---|
| 76 | jpt(np)=j |
---|
| 77 | ELSE |
---|
| 78 | jpt(np)=j+1 |
---|
| 79 | ENDIF |
---|
| 80 | np=np+1 |
---|
| 81 | IF(np.LE.npCFMIP) THEN |
---|
| 82 | GOTO 30 |
---|
| 83 | ENDIF |
---|
| 84 | ELSE |
---|
| 85 | j=j+1 |
---|
| 86 | IF(j.LE.nbp_lat-1) THEN |
---|
| 87 | GOTO 40 |
---|
| 88 | ENDIF |
---|
| 89 | ENDIF |
---|
| 90 | |
---|
| 91 | DO np=1, npCFMIP |
---|
| 92 | WRITE(22,*) lon(ipt(np)), io_lat(jpt(np)) |
---|
| 93 | ENDDO |
---|
| 94 | CLOSE(22) |
---|
| 95 | |
---|
| 96 | OPEN(23, file="pointsCFMIPvsLMDZ.txt") |
---|
| 97 | DO ip=1, npCFMIP |
---|
| 98 | lonGCM(ip)=lon(ipt(ip)) |
---|
| 99 | latGCM(ip)=io_lat(jpt(ip)) |
---|
| 100 | if(jpt(ip).GE.2.AND.jpt(ip).LE.nbp_lat-1) THEN |
---|
| 101 | tabijGCM(ip)=1+(jpt(ip)-2)*nbp_lon+ipt(ip) |
---|
| 102 | else if(jpt(ip).EQ.1) THEN |
---|
| 103 | tabijGCM(ip)=1 |
---|
| 104 | else if(jpt(ip).EQ.nbp_lat) THEN |
---|
| 105 | tabijGCM(ip)=klon_glo |
---|
| 106 | else |
---|
| 107 | print*,'ip jpt tabijGCM',ip,jpt(ip),tabijGCM(ip) |
---|
| 108 | endif |
---|
| 109 | ! PRINT*,'CFMIP ip lon lat tabijGCM',ip,lonGCM(ip),latGCM(ip),tabijGCM(ip) |
---|
| 110 | ENDDO |
---|
| 111 | DO ip=1, npCFMIP |
---|
| 112 | if(lonGCM(ip).EQ.io_lon(1)) lonGCM(ip)=360.+lonGCM(ip) |
---|
| 113 | ENDDO |
---|
| 114 | DO i=1, npCFMIP |
---|
| 115 | WRITE(23,*) i, lonCFMIP(i), latCFMIP(i), lonGCM(i), latGCM(i), tabijGCM(i) |
---|
| 116 | ENDDO |
---|
| 117 | CLOSE(23) |
---|
| 118 | END SUBROUTINE LMDZ_CFMIP_point_locations |
---|
| 119 | |
---|
| 120 | END MODULE CFMIP_point_locations |
---|