4 |
|
|
5 |
contains |
contains |
6 |
|
|
7 |
subroutine writehist(time, vcov, ucov, teta, phi, q, masse, ps) |
subroutine writehist(time, vcov, ucov, teta, phi, masse, ps) |
8 |
|
|
9 |
! From writehist.F 1403 2010-07-01 09:02:53Z |
! From writehist.F 1403 2010-07-01 09:02:53Z |
10 |
! Écriture du fichier histoire au format IOIPSL |
! Écriture du fichier histoire au format IOIPSL |
12 |
! L. Fairhead, LMD, 03/99 |
! L. Fairhead, LMD, 03/99 |
13 |
|
|
14 |
use dimens_m, only: nqmx, llm, jjm |
use dimens_m, only: nqmx, llm, jjm |
|
USE iniadvtrac_m, ONLY: ttext |
|
15 |
use com_io_dyn, only: histid, histvid, histuid |
use com_io_dyn, only: histid, histvid, histuid |
16 |
use paramet_m, only: ip1jm, ip1jmp1, iip1, jjp1 |
use paramet_m, only: ip1jm, ip1jmp1, iip1, jjp1 |
17 |
use temps, only: itau_dyn |
use temps, only: itau_dyn |
25 |
! ucov: vents u covariants |
! ucov: vents u covariants |
26 |
! teta: temperature potentielle |
! teta: temperature potentielle |
27 |
! phi : geopotentiel instantane |
! phi : geopotentiel instantane |
|
! q : traceurs |
|
28 |
! masse: masse |
! masse: masse |
29 |
! ps :pression au sol |
! ps :pression au sol |
30 |
|
|
32 |
|
|
33 |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) |
REAL vcov(ip1jm, llm), ucov(ip1jmp1, llm) |
34 |
REAL teta(ip1jmp1, llm), phi(ip1jmp1, llm) |
REAL teta(ip1jmp1, llm), phi(ip1jmp1, llm) |
35 |
REAL ps(ip1jmp1), masse(ip1jmp1, llm) |
REAL, intent(in):: ps(ip1jmp1), masse(ip1jmp1, llm) |
|
REAL q(ip1jmp1, llm, nqmx) |
|
36 |
integer time |
integer time |
37 |
|
|
38 |
! This routine needs IOIPSL to work |
! This routine needs IOIPSL to work |
39 |
! Variables locales |
! Variables locales |
40 |
|
|
|
integer iq, ii, ll |
|
41 |
integer ndexu(ip1jmp1*llm), ndexv(ip1jm*llm), ndex2d(ip1jmp1) |
integer ndexu(ip1jmp1*llm), ndexv(ip1jm*llm), ndex2d(ip1jmp1) |
42 |
logical ok_sync |
logical ok_sync |
43 |
integer itau_w |
integer itau_w |
73 |
|
|
74 |
call histwrite(histid, 'phi', itau_w, phi) |
call histwrite(histid, 'phi', itau_w, phi) |
75 |
|
|
|
! Traceurs |
|
|
|
|
|
! DO iq=1, nqmx |
|
|
! call histwrite(histid, ttext(iq), itau_w, q(:, :, iq), |
|
|
! . iip1*jjp1*llm, ndexu) |
|
|
! enddo |
|
|
!C |
|
76 |
! Masse |
! Masse |
77 |
|
|
78 |
call histwrite(histid, 'masse', itau_w, masse) |
call histwrite(histid, 'masse', itau_w, masse) |