Changeset 1395 for trunk/LMDZ.MARS/libf/dyn3d/gcm.F
- Timestamp:
- Mar 12, 2015, 12:45:17 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/dyn3d/gcm.F
r1130 r1395 5 5 & nday_r, idissip, iconser, ecritstart, 6 6 & ecritphy 7 use comgeomphy, only: initcomgeomphy7 ! use comgeomphy, only: initcomgeomphy 8 8 IMPLICIT NONE 9 9 … … 143 143 c variables pour l'initialisation de la physique : 144 144 c ------------------------------------------------ 145 INTEGER ngridmx146 PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm )147 REAL zcufi(ngridmx),zcvfi(ngridmx)148 REAL latfi(ngridmx),lonfi(ngridmx)149 REAL airefi(ngridmx)150 SAVE latfi, lonfi, airefi151 INTEGER i,j145 ! INTEGER ngridmx 146 ! PARAMETER( ngridmx = 2+(jjm-1)*iim - 1/jjm ) 147 ! REAL zcufi(ngridmx),zcvfi(ngridmx) 148 ! REAL latfi(ngridmx),lonfi(ngridmx) 149 ! REAL airefi(ngridmx) 150 ! SAVE latfi, lonfi, airefi 151 ! INTEGER i,j 152 152 153 153 c----------------------------------------------------------------------- … … 168 168 !#ifdef CPP_PHYS 169 169 CALL init_phys_lmdz(iim,jjp1,llm,1,(/(jjm-1)*iim+2/)) 170 call initcomgeomphy 170 ! call initcomgeomphy ! now done in iniphysiq 171 171 !#endif 172 172 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! … … 227 227 228 228 ! IF (call_iniphys.and.(iflag_phys==1.or.iflag_phys>=100)) THEN 229 latfi(1)=rlatu(1)230 lonfi(1)=0.231 zcufi(1) = cu(1)232 zcvfi(1) = cv(1)233 DO j=2,jjm234 DO i=1,iim235 latfi((j-2)*iim+1+i)= rlatu(j)236 lonfi((j-2)*iim+1+i)= rlonv(i)237 zcufi((j-2)*iim+1+i) = cu((j-1)*iip1+i)238 zcvfi((j-2)*iim+1+i) = cv((j-1)*iip1+i)239 ENDDO240 ENDDO241 latfi(ngridmx)= rlatu(jjp1)242 lonfi(ngridmx)= 0.243 zcufi(ngridmx) = cu(ip1jm+1)244 zcvfi(ngridmx) = cv(ip1jm-iim)245 246 ! build airefi(), mesh area on physics grid247 CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi)248 ! Poles are single points on physics grid249 airefi(1)=airefi(1)*iim250 airefi(ngridmx)=airefi(ngridmx)*iim229 ! latfi(1)=rlatu(1) 230 ! lonfi(1)=0. 231 ! zcufi(1) = cu(1) 232 ! zcvfi(1) = cv(1) 233 ! DO j=2,jjm 234 ! DO i=1,iim 235 ! latfi((j-2)*iim+1+i)= rlatu(j) 236 ! lonfi((j-2)*iim+1+i)= rlonv(i) 237 ! zcufi((j-2)*iim+1+i) = cu((j-1)*iip1+i) 238 ! zcvfi((j-2)*iim+1+i) = cv((j-1)*iip1+i) 239 ! ENDDO 240 ! ENDDO 241 ! latfi(ngridmx)= rlatu(jjp1) 242 ! lonfi(ngridmx)= 0. 243 ! zcufi(ngridmx) = cu(ip1jm+1) 244 ! zcvfi(ngridmx) = cv(ip1jm-iim) 245 ! 246 ! ! build airefi(), mesh area on physics grid 247 ! CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,aire,airefi) 248 ! ! Poles are single points on physics grid 249 ! airefi(1)=airefi(1)*iim 250 ! airefi(ngridmx)=airefi(ngridmx)*iim 251 251 252 252 ! Initialisation de la physique: pose probleme quand on tourne … … 254 254 ! Il faut une cle CPP_PHYS 255 255 !#ifdef CPP_PHYS 256 ! CALL iniphysiq( ngridmx,llm,daysec,day_ini,dtphys/nsplit_phys,257 CALL iniphysiq( ngridmx,llm,daysec,day_ini,dtphys,258 & latfi,lonfi,airefi,zcufi,zcvfi,rad,g,r,cpp,256 ! CALL iniphysiq(iim,jjm,llm,daysec,day_ini,dtphys/nsplit_phys, 257 CALL iniphysiq(iim,jjm,llm,daysec,day_ini,dtphys, 258 & rlatu,rlonv,aire,cu,cv,rad,g,r,cpp, 259 259 & 1) 260 260 ! & iflag_phys)
Note: See TracChangeset
for help on using the changeset viewer.