source: codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/write_etat0.f90 @ 316

Last change on this file since 316 was 262, checked in by ymipsl, 10 years ago

Save starting iteration itau0 in dynamico restart file.
Timeloop now begin from itau0 which can be read from start file or set to 0 by default.

YM
File size: 942 bytes
Line 
1MODULE write_etat0_mod
2
3
4
5CONTAINS
6
7  SUBROUTINE write_etat0(it,f_ps,f_phis,f_theta_rhodz,f_u, f_q)
8  USE icosa
9  USE restart_mod
10  USE wind_mod
11  USE write_field
12  IMPLICIT NONE
13    INTEGER,INTENT(IN)    :: it
14    TYPE(t_field),POINTER :: f_ps(:)
15    TYPE(t_field),POINTER :: f_phis(:)
16    TYPE(t_field),POINTER :: f_theta_rhodz(:)
17    TYPE(t_field),POINTER :: f_u(:)
18    TYPE(t_field),POINTER :: f_q(:)
19 
20    TYPE(t_field),POINTER,SAVE :: f_ulon(:)
21    TYPE(t_field),POINTER,SAVE :: f_ulat(:)
22   
23   
24    CALL allocate_field(f_ulon,field_t,type_real,llm,name='ulon')
25    CALL allocate_field(f_ulat,field_t,type_real,llm,name='ulat')
26
27    CALL transfert_request(f_u,req_e1_vect)
28    CALL un2ulonlat(f_u, f_ulon, f_ulat)
29   
30    CALL write_restart(it,f_ps,f_phis,f_theta_rhodz,f_u, f_ulon, f_ulat, f_q)
31
32    CALL deallocate_field(f_ulon)
33    CALL deallocate_field(f_ulat)
34   
35  END SUBROUTINE write_etat0
36   
37END MODULE write_etat0_mod
Note: See TracBrowser for help on using the repository browser.