/[lmdze]/trunk/dyn3d/tourpot.f
ViewVC logotype

Annotation of /trunk/dyn3d/tourpot.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 104 - (hide annotations)
Thu Sep 4 10:05:52 2014 UTC (9 years, 8 months ago) by guez
File size: 1348 byte(s) 
Removed procedure sortvarc0. Called sortvarc with an additional
argument resetvarc instead. (Following LMDZ.) Moved current time
computations and some printing statements from sortvarc to
caldyn. Could then remove arguments itau and time_0 of sortvarc, and
could remove "use dynetat0". Better to keep "dynetat0.f" as a gcm-only
file.

Moved some variables from module ener to module sortvarc.

Split file "mathelp.f" into single-procedure files.

Removed unused argument nadv of adaptdt. Removed dimension arguments
of bernoui.

Removed unused argument nisurf of interfoce_lim. Changed the size of
argument lmt_sst of interfoce_lim from klon to knon. Removed case when
newlmt is false.

dynredem1 is called only once in each run, either ce0l or gcm. So
variable nb in call to nf95_put_var was always 1. Removed variable nb.

Removed dimension arguments of calcul_fluxs. Removed unused arguments
precip_rain, precip_snow, snow of calcul_fluxs. Changed the size of
all the arrays in calcul_fluxs from klon to knon.

Removed dimension arguments of fonte_neige. Changed the size of all
the arrays in fonte_neige from klon to knon.

Changed the size of arguments tsurf and tsurf_new of interfsurf_hq
from klon to knon. Changed the size of argument ptsrf of soil from
klon to knon.
1 guez 78 module tourpot_m
2 guez 3
3 guez 78 IMPLICIT NONE
4 guez 3
5 guez 78 contains
6 guez 27
7 guez 78 SUBROUTINE tourpot(vcov, ucov, massebxy, vorpot)
8 guez 3
9 guez 78 ! From LMDZ4/libf/dyn3d/tourpot.F, version 1.1.1.1 2004/05/19 12:53:06
10 guez 3
11 guez 79 ! Author: P. Le Van
12 guez 78 ! Objet : calcul du tourbillon potentiel
13 guez 3
14 guez 104 USE comgeom, ONLY: fext_2d
15 guez 78 USE dimens_m, ONLY: iim, jjm, llm
16     use filtreg_m, only: filtreg
17 guez 3
18 guez 79 REAL, intent(in):: vcov(:, :, :) ! (iim + 1, jjm, llm)
19     REAL, intent(in):: ucov(:, :, :) ! (iim + 1, jjm + 1, llm)
20     REAL, intent(in):: massebxy(:, :, :) ! (iim + 1, jjm, llm) mass of grid cell
21 guez 3
22 guez 79 real, intent(out):: vorpot(:, :, :) ! (iim + 1, jjm, llm)
23 guez 78 ! = (Filtre(d(vcov)/dx - d(ucov)/dy) + fext) / massebxy
24 guez 3
25 guez 78 ! Local:
26 guez 104
27 guez 79 REAL rot(iim + 1, jjm, llm)
28     ! relative vorticity multiplied by cell area, in m2 s-1
29 guez 3
30 guez 79 INTEGER l, i, j
31    
32 guez 78 !---------------------------------------------------------------
33 guez 3
34 guez 78 ! Calcul du rotationnel du vent puis filtrage
35 guez 3
36 guez 79 forall (i = 1: iim, j = 1: jjm) rot(i, j, :) &
37     = vcov(i + 1, j, :) - vcov(i, j, :) + ucov(i, j + 1, :) - ucov(i, j, :)
38     rot(iim + 1, :, :) = rot(1, :, :)
39 guez 3
40 guez 78 CALL filtreg(rot, jjm, llm, 2, 1, .FALSE.)
41 guez 3
42 guez 79 forall (l = 1: llm) vorpot(:iim, :, l) &
43     = (rot(:iim, :, l) + fext_2d(:iim, :)) / massebxy(:iim, :, l)
44     vorpot(iim + 1, :, :)= vorpot(1, :, :)
45 guez 3
46 guez 78 END SUBROUTINE tourpot
47 guez 3
48 guez 78 end module tourpot_m
  ViewVC Help
Powered by ViewVC 1.1.21