source: LMDZ6/trunk/libf/dyn3d_common/ugeostr.f90 @ 5396

Last change on this file since 5396 was 5285, checked in by abarral, 7 weeks ago

As discussed internally, remove generic ONLY: ... for new _mod_h modules

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.4 KB
RevLine 
[524]1!
[1279]2! $Id: ugeostr.f90 5285 2024-10-28 13:33:29Z evignon $
[524]3!
[1474]4subroutine ugeostr(phi,ucov)
[524]5
[1474]6  ! Calcul du vent covariant geostrophique a partir du champ de
7  ! geopotentiel.
8  ! We actually compute: (1 - cos^8 \phi) u_g
9  ! to have a wind going smoothly to 0 at the equator.
10  ! We assume that the surface pressure is uniform so that model
11  ! levels are pressure levels.
[524]12
[5281]13  USE comgeom2_mod_h
[2597]14  use comconst_mod, only: omeg, rad
[5281]15
[5271]16  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
[5285]17USE paramet_mod_h
[5271]18implicit none
[524]19
[5271]20
[5272]21
[524]22
[1474]23  real ucov(iip1,jjp1,llm),phi(iip1,jjp1,llm)
24  real um(jjm,llm),fact,u(iip1,jjm,llm)
25  integer i,j,l
[524]26
[1474]27  real zlat
[524]28
[1474]29  um(:,:)=0 ! initialize um()
[524]30
[1474]31  DO j=1,jjm
[1279]32
[1474]33     if (abs(sin(rlatv(j))).lt.1.e-4) then
34        zlat=1.e-4
35     else
36        zlat=rlatv(j)
37     endif
38     fact=cos(zlat)
39     fact=fact*fact
40     fact=fact*fact
41     fact=fact*fact
42     fact=(1.-fact)/ &
43          (2.*omeg*sin(zlat)*(rlatu(j+1)-rlatu(j)))
44     fact=-fact/rad
45     DO l=1,llm
46        DO i=1,iim
47           u(i,j,l)=fact*(phi(i,j+1,l)-phi(i,j,l))
48           um(j,l)=um(j,l)+u(i,j,l)/REAL(iim)
49        ENDDO
50     ENDDO
51  ENDDO
52  call dump2d(jjm,llm,um,'Vent-u geostrophique')
[524]53
[1474]54  !   calcul des champ de vent:
[524]55
[1474]56  DO l=1,llm
57     DO i=1,iip1
58        ucov(i,1,l)=0.
59        ucov(i,jjp1,l)=0.
60     end DO
61     DO  j=2,jjm
62        DO  i=1,iim
63           ucov(i,j,l) = 0.5*(u(i,j,l)+u(i,j-1,l))*cu(i,j)
64        end DO
65        ucov(iip1,j,l)=ucov(1,j,l)
66     end DO
67  end DO
[524]68
[1474]69  print *, 301
[524]70
[1474]71end subroutine ugeostr
Note: See TracBrowser for help on using the repository browser.