trunk/dyn3d/writehist.f

module writehist_m

  implicit none

contains

  subroutine writehist(time, vcov, ucov, teta, phi, masse, ps)

    ! From writehist.F , revision 1403, 2010-07-01 09:02:53Z
    ! Écriture du fichier histoire
    ! L. Fairhead, LMD, 03/99

    use covnat_m, only: covnat
    use dimens_m, only: llm
    use histsync_m, only: histsync
    use histwrite_m, only: histwrite
    use inithist_m, only: histid, histvid, histuid
    use paramet_m, only: ip1jm, ip1jmp1
    use temps, only: itau_dyn

    integer, intent(in):: time ! temps de l'ecriture
    REAL, intent(in):: vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vent covariant
    REAL, intent(in):: teta(ip1jmp1, llm) ! temperature potentielle
    REAL, intent(in):: phi(ip1jmp1, llm) ! geopotentiel instantane
    REAL, intent(in):: masse(ip1jmp1, llm)
    REAL, intent(in):: ps(ip1jmp1) ! pression au sol

    ! Local:
    logical ok_sync
    integer itau_w
    REAL vnat(ip1jm, llm), unat(ip1jmp1, llm)

    !---------------------------------------------------------------------

    ! Initialisations

    ok_sync =.TRUE.
    itau_w = itau_dyn + time
    ! Passage aux composantes naturelles du vent
    call covnat(llm, ucov, vcov, unat, vnat)

    ! Appels a histwrite pour l'ecriture des variables a sauvegarder

    ! Vents U

    call histwrite(histuid, 'u', itau_w, unat)

    ! Vents V

    call histwrite(histvid, 'v', itau_w, vnat)

    ! Temperature potentielle

    call histwrite(histid, 'teta', itau_w, teta)

    ! Geopotentiel

    call histwrite(histid, 'phi', itau_w, phi)

    ! Masse

    call histwrite(histid, 'masse', itau_w, masse)

    ! Pression au sol
    call histwrite(histid, 'ps', itau_w, ps)

    if (ok_sync) then
       call histsync(histid)
       call histsync(histvid)
       call histsync(histuid)
    endif

  end subroutine writehist

end module writehist_m
1	guez	69	module writehist_m
2	guez	56
3	guez	69	implicit none
4	guez	56
5	guez	69	contains
6	guez	56
7	guez	138	subroutine writehist(time, vcov, ucov, teta, phi, masse, ps)
8	guez	56
9	guez	212	! From writehist.F , revision 1403, 2010-07-01 09:02:53Z
10			! Écriture du fichier histoire
11	guez	69	! L. Fairhead, LMD, 03/99
12	guez	56
13	guez	212	use covnat_m, only: covnat
14	guez	178	use dimens_m, only: llm
15	guez	212	use histsync_m, only: histsync
16			use histwrite_m, only: histwrite
17			use inithist_m, only: histid, histvid, histuid
18	guez	178	use paramet_m, only: ip1jm, ip1jmp1
19	guez	69	use temps, only: itau_dyn
20	guez	56
21	guez	212	integer, intent(in):: time ! temps de l'ecriture
22			REAL, intent(in):: vcov(ip1jm, llm), ucov(ip1jmp1, llm) ! vent covariant
23			REAL, intent(in):: teta(ip1jmp1, llm) ! temperature potentielle
24			REAL, intent(in):: phi(ip1jmp1, llm) ! geopotentiel instantane
25			REAL, intent(in):: masse(ip1jmp1, llm)
26			REAL, intent(in):: ps(ip1jmp1) ! pression au sol
27	guez	56
28	guez	212	! Local:
29	guez	69	logical ok_sync
30			integer itau_w
31			REAL vnat(ip1jm, llm), unat(ip1jmp1, llm)
32
33			!---------------------------------------------------------------------
34
35			! Initialisations
36
37			ok_sync =.TRUE.
38			itau_w = itau_dyn + time
39			! Passage aux composantes naturelles du vent
40			call covnat(llm, ucov, vcov, unat, vnat)
41
42			! Appels a histwrite pour l'ecriture des variables a sauvegarder
43
44			! Vents U
45
46			call histwrite(histuid, 'u', itau_w, unat)
47
48			! Vents V
49
50			call histwrite(histvid, 'v', itau_w, vnat)
51
52			! Temperature potentielle
53
54			call histwrite(histid, 'teta', itau_w, teta)
55
56			! Geopotentiel
57
58			call histwrite(histid, 'phi', itau_w, phi)
59
60			! Masse
61
62			call histwrite(histid, 'masse', itau_w, masse)
63
64			! Pression au sol
65			call histwrite(histid, 'ps', itau_w, ps)
66
67			if (ok_sync) then
68			call histsync(histid)
69			call histsync(histvid)
70			call histsync(histuid)
71			endif
72
73			end subroutine writehist
74
75			end module writehist_m