source: codes/icosagcm/trunk/src/write_etat0.f90 @ 347

Last change on this file since 347 was 347, checked in by dubos, 9 years ago

Synced with aquaplanet branch HEAT@45 - tested with DCMIP4

File size: 946 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_mod
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.