! $Header$

IF (callinigrads) THEN
   string10='dyn'
   CALL inigrads(1,iip1 &
         ,rlonv,180./pi,-180.,180.,jjp1,rlatu,-90.,90.,180./pi &
         ,llm,presnivs,1. &
         ,dtvr*iperiod,string10,'dyn_zon ')

  callinigrads=.FALSE.


endif

string10='ps'
CALL wrgrads(1,1,ps,string10,string10)

string10='u'
CALL wrgrads(1,llm,unat,string10,string10)
string10='v'
CALL wrgrads(1,llm,vnat,string10,string10)
string10='teta'
CALL wrgrads(1,llm,teta,string10,string10)
DO iq=1,nqtot
   string10='q'
   WRITE(string10(2:2),'(i1)') iq
   CALL wrgrads(1,llm,q(:,:,iq),string10,string10)
enddo