Changeset 360 for LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90
- Timestamp:
- May 24, 2002, 3:54:35 PM (22 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ.3.3/branches/rel-LF/libf/phylmd/interface_surf.F90
r353 r360 354 354 & tq_cdrag, petAcoef, peqAcoef, petBcoef, peqBcoef, & 355 355 & precip_rain, precip_snow, sollwdown, swnet, swdown, & 356 & tsurf, p1lay/100., ps , radsol, &356 & tsurf, p1lay/100., ps/100., radsol, & 357 357 & evap, fluxsens, fluxlat, & 358 358 & tsol_rad, tsurf_new, alb_new, alblw, & … … 844 844 ! Attention aux poles 845 845 ! 846 !!$ do igrid = 1, knon847 !!$ index = ktindex(igrid)848 !!$ ij = index - int((index-1)/iim)*iim - 1849 !!$ jj = 2 + int((index-1)/iim)850 !!$ if (mod(index,iim) == 1 ) then851 !!$ jj = 1 + int((index-1)/iim)852 !!$ ij = iim853 !!$ endif854 !!$ correspond(ij,jj) = igrid855 !!$ write(50,*)'igrid, i, j =',igrid,ij,jj856 !!$ enddo857 ! Pb de correspondances de grilles!858 !!$ do igrid = 1, knon859 !!$ index = ktindex(igrid)860 !!$ ij = ig(index)861 !!$ jj = jg(index)862 !!$ correspond(ij,jj) = igrid863 !!$ write(51,*)'igrid, i, j =',igrid,ij,jj864 !!$ enddo865 846 do igrid = 1, knon 866 847 index = ktindex(igrid) … … 869 850 correspond(ij,jj) = igrid 870 851 enddo 871 872 !!$ index = 0873 !!$ do jj = 1, jjm+1874 !!$ do ij = 1, iim875 !!$ index = index + 1876 !!$ correspond(ij,jj) = index877 !!$ enddo878 !!$ enddo879 852 880 853 ! Allouer et initialiser le tableau de coordonnees du sol … … 973 946 index = iglob + offset(i) 974 947 ireal = (min(max(1, index - iim + 1), klon)) 975 ! if (index <= 1) index = 1976 ! if (index >= klon) index = klon977 948 if (pctsrf(ireal, is_ter) > EPSFRA) then 978 949 jj = int((index - 1)/iim) + 1 979 950 ij = index - (jj - 1) * iim 980 !!$ ij = index - int((index-1)/iim)*iim - 1981 !!$ jj = 2 + int((index-1)/iim)982 !!$ if (mod(index,iim) == 1 ) then983 !!$ jj = 1 + int((index-1)/iim)984 !!$ ij = iim985 !!$ endif986 !!$! Pb de correspondances de grilles!987 !!$ ij = ig(index)988 !!$ jj = jg(index)989 !!$! write(*,*)'correspond',igrid, ireal,index,ij,jj990 ! if ( ij >= 1 .and. ij <= iim .and. jj >= 1 .and. jj <= jjm) then991 ! write(*,*)'correspond',igrid, ireal,index,ij,jj992 951 neighbours(igrid, i) = correspond(ij, jj) 993 ! endif994 952 endif 995 953 enddo 996 954 enddo 997 998 write(*,*)'Neighbours = '999 1000 write(*,*)neighbours(1,8), neighbours(1,1),neighbours(1,2)1001 write(*,*)neighbours(1,7), ktindex(1), neighbours(1,3)1002 write(*,*)neighbours(1,6), neighbours(1,5),neighbours(1,4)1003 1004 write(*,*)neighbours(250,8), neighbours(250,1),neighbours(250,2)1005 write(*,*)neighbours(250,7), ktindex(250), neighbours(250,3)1006 write(*,*)neighbours(250,6), neighbours(250,5),neighbours(250,4)1007 1008 OPEN (unit=12, file="neighbours.9671")1009 DO i=1,knon1010 WRITE(12,*) '-----------------------------'1011 WRITE(12,'(I7,f8.5, " ",3I6)') knon, contfrac(i), &1012 & neighbours(i,8), neighbours(i,1), neighbours(i,2)1013 WRITE(12,'(f10.5," ",3I6)') lalo(i,2), neighbours(i,7), &1014 & ktindex(i), neighbours(i,3)1015 WRITE(12,'(f10.5," ",3I6)') lalo(i,1), neighbours(i,6), &1016 & neighbours(i,5), neighbours(i,3)1017 ENDDO1018 CLOSE(12)1019 1020 955 1021 956 ! … … 1064 999 & zlev, u1_lay, v1_lay, spechum, temp_air, epot_air, ccanopy, & 1065 1000 & cdrag, petA_orc, peqA_orc, petB_orc, peqB_orc, & 1066 & precip_rain, precip_snow, lwdown, swnet, swdown, p 1lay, &1001 & precip_rain, precip_snow, lwdown, swnet, swdown, ps, & 1067 1002 & evap, fluxsens, fluxlat, coastalflow, riverflow, & 1068 1003 & tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, & … … 1075 1010 & zlev, u1_lay, v1_lay, spechum, temp_air, epot_air, ccanopy, & 1076 1011 & cdrag, petA_orc, peqA_orc, petB_orc, peqB_orc, & 1077 & precip_rain, precip_snow, lwdown, swnet, swdown, p 1lay, &1012 & precip_rain, precip_snow, lwdown, swnet, swdown, ps, & 1078 1013 & evap, fluxsens, fluxlat, coastalflow, riverflow, & 1079 1014 & tsol_rad, tsurf_new, qsurf, albedo_out, emis_new, z0_new, & … … 1082 1017 bidule=0. 1083 1018 bidule(1:knon)=riverflow(1:knon) 1084 call gath2cpl(bidule, tmp_rriv, klon, knon,iim,jjm, 1019 call gath2cpl(bidule, tmp_rriv, klon, knon,iim,jjm,ktindex) 1085 1020 bidule=0. 1086 1021 bidule(1:knon)=coastalflow(1:knon) 1087 call gath2cpl(bidule, tmp_rcoa, klon, knon,iim,jjm, 1022 call gath2cpl(bidule, tmp_rcoa, klon, knon,iim,jjm,ktindex) 1088 1023 alb_new(1:knon) = albedo_out(1:knon,1) 1089 1024 alblw(1:knon) = albedo_out(1:knon,2)
Note: See TracChangeset
for help on using the changeset viewer.