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	module writehist_m
2
3	implicit none
4
5	contains
6
7	subroutine writehist(time, vcov, ucov, teta, phi, masse, ps)
8
9	! From writehist.F , revision 1403, 2010-07-01 09:02:53Z
10	! Écriture du fichier histoire
11	! L. Fairhead, LMD, 03/99
12
13	use covnat_m, only: covnat
14	use dimens_m, only: llm
15	use histsync_m, only: histsync
16	use histwrite_m, only: histwrite
17	use inithist_m, only: histid, histvid, histuid
18	use paramet_m, only: ip1jm, ip1jmp1
19	use temps, only: itau_dyn
20
21	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
28	! Local:
29	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