- Timestamp:
- Jul 24, 2024, 4:23:34 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/read_surface.F90
r5116 r5117 18 18 character*10 varname 19 19 20 realtmp_dyn(iip1,jjp1)21 realtmp_dyn_glo(nbp_lon+1,nbp_lat)20 REAL tmp_dyn(iip1,jjp1) 21 REAL tmp_dyn_glo(nbp_lon+1,nbp_lat) 22 22 REAL tmp_dyn_invers(iip1,jjp1) 23 realtmp_dyn_invers_glo(nbp_lon+1,nbp_lat)24 realtmp_fi(klon)25 realtmp_fi_glo(klon_glo)26 realsurfa(klon,5)27 realsurfa_glo(klon_glo,5)23 REAL tmp_dyn_invers_glo(nbp_lon+1,nbp_lat) 24 REAL tmp_fi(klon) 25 REAL tmp_fi_glo(klon_glo) 26 REAL surfa(klon,5) 27 REAL surfa_glo(klon_glo,5) 28 28 29 integerncid30 integervarid31 integerrcode32 integerstart(2),count(2),status33 integeri,j,l,ig29 INTEGER ncid 30 INTEGER varid 31 INTEGER rcode 32 INTEGER start(2),count(2),status 33 INTEGER i,j,l,ig 34 34 character*1 str1 35 35 36 36 !JE20140526<< 37 37 character*4 :: latstr,aux4s 38 logical:: outcycle, isinversed39 real, dimension(jjp1) :: lats40 real, dimension(nbp_lat) :: lats_glo38 LOGICAL :: outcycle, isinversed 39 REAL, DIMENSION(jjp1) :: lats 40 REAL, DIMENSION(nbp_lat) :: lats_glo 41 41 REAL :: rcode2 42 integer, dimension(1) :: startj,endj42 INTEGER, DIMENSION(1) :: startj,endj 43 43 !JE20140526>> 44 44 !$OMP MASTER … … 53 53 isinversed=.FALSE. 54 54 do i=1,5 55 if(i==1) aux4s='latu'56 if(i==2) aux4s='LATU'57 if(i==3) aux4s='LatU'58 if(i==4) aux4s='Latu'59 if(i==5) aux4s='latU'55 IF (i==1) aux4s='latu' 56 IF (i==2) aux4s='LATU' 57 IF (i==3) aux4s='LatU' 58 IF (i==4) aux4s='Latu' 59 IF (i==5) aux4s='latU' 60 60 status = nf90_inq_varid (ncid, aux4s, rcode) 61 61 ! print *,'stat,i',status,i,outcycle,aux4s 62 62 ! print *,'ifclause',status.NE. nf90_noerr ,outcycle == .FALSE. 63 IF ((. not.(status/= nf90_noerr) ).and.( .not. outcycle )) THEN63 IF ((.NOT.(status/= nf90_noerr) ).AND.( .NOT. outcycle )) THEN 64 64 outcycle=.TRUE. 65 65 latstr=aux4s … … 78 78 79 79 ! check if netcdf is latitude inversed or not. 80 if(lats_glo(1)<lats_glo(2)) isinversed=.TRUE.80 IF (lats_glo(1)<lats_glo(2)) isinversed=.TRUE. 81 81 ! JE20140526>> 82 82 … … 118 118 !JE20140526<< 119 119 ! CALL gr_dyn_fi(1, iip1, jjp1, klon, tmp_dyn_invers, tmp_fi) 120 if(isinversed) THEN120 IF (isinversed) THEN 121 121 CALL gr_dyn_fi(1, nbp_lon+1, nbp_lat, klon_glo, & 122 122 tmp_dyn_invers_glo, tmp_fi_glo)
Note: See TracChangeset
for help on using the changeset viewer.