/[lmdze]/trunk/libf/dyn3d/inigeom.f90
ViewVC logotype

Diff of /trunk/libf/dyn3d/inigeom.f90

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 57 by guez, Mon Jan 30 12:54:02 2012 UTC revision 60 by guez, Mon Jan 30 14:37:26 2012 UTC
# Line 76  contains Line 76  contains
76      REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm      REAL eps, x1, xo1, f, df, xdm, y1, yo1, ydm
77      REAL coslatm, coslatp, radclatm, radclatp      REAL coslatm, coslatp, radclatm, radclatp
78      REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m      REAL, dimension(iip1, jjp1):: cuij1, cuij2, cuij3, cuij4 ! in m
79      REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4      REAL, dimension(iip1, jjp1):: cvij1, cvij2, cvij3, cvij4 ! in m
80      REAL rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)      REAL rlatu1(jjm), yprimu1(jjm), rlatu2(jjm), yprimu2(jjm)
81      real yprimv(jjm), yprimu(jjp1)      real yprimv(jjm), yprimu(jjp1)
82      REAL gamdi_gdiv, gamdi_grot, gamdi_h      REAL gamdi_gdiv, gamdi_grot, gamdi_h
# Line 332  contains Line 332  contains
332                 + cvij3(i, j))                 + cvij3(i, j))
333            cuv(i, j) = 0.5 * (cuij2(i, j) + cuij3(i, j) + cuij1(i, j + 1) &            cuv(i, j) = 0.5 * (cuij2(i, j) + cuij3(i, j) + cuij1(i, j + 1) &
334                 + cuij4(i, j + 1))                 + cuij4(i, j + 1))
335            unscv2_2d(i, j) = 1. / (cv_2d(i, j) * cv_2d(i, j))            unscv2_2d(i, j) = 1. / cv_2d(i, j)**2
336         END DO         END DO
337         DO i = 1, iim         DO i = 1, iim
338            cuvsurcv_2d(i, j) = airev_2d(i, j) * unscv2_2d(i, j)            cuvsurcv_2d(i, j) = airev_2d(i, j) * unscv2_2d(i, j)
# Line 356  contains Line 356  contains
356         DO i = 1, iim         DO i = 1, iim
357            cu_2d(i, j) = 0.5 * (cuij1(i, j) + cuij4(i + 1, j) + cuij2(i, j) &            cu_2d(i, j) = 0.5 * (cuij1(i, j) + cuij4(i + 1, j) + cuij2(i, j) &
358                 + cuij3(i + 1, j))                 + cuij3(i + 1, j))
359            unscu2_2d(i, j) = 1. / (cu_2d(i, j) * cu_2d(i, j))            unscu2_2d(i, j) = 1. / cu_2d(i, j)**2
360            cvusurcu_2d(i, j) = aireu_2d(i, j) * unscu2_2d(i, j)            cvusurcu_2d(i, j) = aireu_2d(i, j) * unscu2_2d(i, j)
361            cusurcvu_2d(i, j) = 1. / cvusurcu_2d(i, j)            cusurcvu_2d(i, j) = 1. / cvusurcu_2d(i, j)
362            cvuscugam1_2d(i, j) = cvusurcu_2d(i, j)**(-gamdi_gdiv)            cvuscugam1_2d(i, j) = cvusurcu_2d(i, j)**(-gamdi_gdiv)
Legend:
Removed from v.57  
changed lines
  Added in v.60
  ViewVC Help
Powered by ViewVC 1.1.21