/[lmdze]/trunk/Sources/dyn3d/Guide/guide.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/Guide/guide.f

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

revision 20 by guez, Wed Oct 15 16:19:57 2008 UTC revision 22 by guez, Fri Jul 31 15:18:47 2009 UTC
# Line 54  SUBROUTINE guide(itau, ucov, vcov, teta, Line 54  SUBROUTINE guide(itau, ucov, vcov, teta,
54    
55    !=======================================================================    !=======================================================================
56    
57    !   ...  Dans inigeom , nouveaux calculs pour les elongations  cu , cv    ! Dans inigeom , nouveaux calculs pour les elongations  cu , cv
58    !        et possibilite d'appeler une fonction f(y)  a derivee tangente    ! et possibilite d'appeler une fonction f(y)  a derivee tangente
59    !        hyperbolique a la  place de la fonction a derivee sinusoidale.            ! hyperbolique a la  place de la fonction a derivee sinusoidale.        
60    
61    !   ...  Possibilite de choisir le shema de Van-leer pour l'advection de    !   ...  Possibilite de choisir le shema de Van-leer pour l'advection de
62    !         q  , en faisant iadv = 10  dans   traceur  (29/04/97) .    !         q  , en faisant iadv = 10  dans   traceur  (29/04/97) .
# Line 249  SUBROUTINE guide(itau, ucov, vcov, teta, Line 249  SUBROUTINE guide(itau, ucov, vcov, teta,
249       PRINT *, 'nlev', nlev       PRINT *, 'nlev', nlev
250       rcod = nf90_close(ncidpl)       rcod = nf90_close(ncidpl)
251       !   Lecture du premier etat des reanalyses.       !   Lecture du premier etat des reanalyses.
252       CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, masserea2, &       CALL read_reanalyse(1, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &
253            psrea2, 1, nlev)            masserea2, psrea2, 1, nlev)
254       qrea2(:, :) = max(qrea2(:, :), 0.1)       qrea2(:, :) = max(qrea2(:, :), 0.1)
255    
256    
# Line 286  SUBROUTINE guide(itau, ucov, vcov, teta, Line 286  SUBROUTINE guide(itau, ucov, vcov, teta,
286               count_no_rea, ' non lectures'               count_no_rea, ' non lectures'
287          step_rea = step_rea + 1          step_rea = step_rea + 1
288          itau_test = itau          itau_test = itau
289          CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, qrea2, &          CALL read_reanalyse(step_rea, ps, ucovrea2, vcovrea2, tetarea2, &
290               masserea2, psrea2, 1, nlev)               qrea2, masserea2, psrea2, 1, nlev)
291          qrea2(:, :) = max(qrea2(:, :), 0.1)          qrea2(:, :) = max(qrea2(:, :), 0.1)
292          factt = dtvr*iperiod/daysec          factt = dtvr*iperiod/daysec
293          ztau(:) = factt/max(alpha_t(:), 1.E-10)          ztau(:) = factt/max(alpha_t(:), 1.E-10)
Legend:
Removed from v.20  
changed lines
  Added in v.22
  ViewVC Help
Powered by ViewVC 1.1.21