1 |
|
|
! |
2 |
|
|
! $Header$ |
3 |
|
|
! |
4 |
|
|
if (callinigrads) then |
5 |
|
|
|
6 |
|
|
string10='dyn' |
7 |
|
|
call inigrads(1,iip1 |
8 |
|
|
s ,rlonv,180./pi,-180.,180.,jjp1,rlatu,-90.,90.,180./pi |
9 |
|
|
s ,llm,presnivs,1. |
10 |
|
|
s ,dtvr*iperiod,string10,'dyn_zon ') |
11 |
|
|
|
12 |
|
|
callinigrads=.false. |
13 |
|
|
|
14 |
|
|
|
15 |
|
|
endif |
16 |
|
|
|
17 |
|
|
string10='ps' |
18 |
|
|
CALL wrgrads(1,1,ps,string10,string10) |
19 |
|
|
|
20 |
|
|
string10='u' |
21 |
|
|
CALL wrgrads(1,llm,unat,string10,string10) |
22 |
|
|
string10='v' |
23 |
|
|
CALL wrgrads(1,llm,vnat,string10,string10) |
24 |
|
|
string10='teta' |
25 |
|
|
CALL wrgrads(1,llm,teta,string10,string10) |
26 |
|
|
do iq=1,nqtot |
27 |
|
|
string10='q' |
28 |
|
|
write(string10(2:2),'(i1)') iq |
29 |
|
|
CALL wrgrads(1,llm,q(:,:,iq),string10,string10) |
30 |
|
|
enddo |
31 |
|
|
|