Changeset 916


Ignore:
Timestamp:
Feb 27, 2008, 4:19:37 PM (17 years ago)
Author:
Laurent Fairhead
Message:

travail de Yann sur dynzon
"débranchement" de dyn_hist, dyn_histv, dyn_hist_ave
SD
LF

Location:
LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar
Files:
5 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/bilan_dyn_p.F

    r764 r916  
    1515      use misc_mod
    1616      use write_field
     17      USE Write_Field_p
    1718      IMPLICIT NONE
    1819
     
    164165      type(Request) :: Req
    165166
     167! definition du domaine d'ecriture pour le rebuild
     168
     169      INTEGER,DIMENSION(1) :: ddid
     170      INTEGER,DIMENSION(1) :: dsg
     171      INTEGER,DIMENSION(1) :: dsl
     172      INTEGER,DIMENSION(1) :: dpf
     173      INTEGER,DIMENSION(1) :: dpl
     174      INTEGER,DIMENSION(1) :: dhs
     175      INTEGER,DIMENSION(1) :: dhe
     176     
     177      INTEGER :: bilan_dyn_domain_id
     178
    166179
    167180c=====================================================================
     
    233246      jje=jj_end
    234247      jjn=jj_nb
    235       if (pole_sud) jjn=jj_nb-1
     248      IF (pole_sud) THEN
     249        jjn=jj_nb-1
     250        jje=jj_end-1
     251      ENDIF
     252
     253      ddid=(/ 2 /)
     254      dsg=(/ jjm /)
     255      dsl=(/ jjn /)
     256      dpf=(/ jjb /)
     257      dpl=(/ jje /)
     258      dhs=(/ 0 /)
     259      dhe=(/ 0 /)
     260
     261      call flio_dom_set(mpi_size,mpi_rank,ddid,dsg,dsl,dpf,dpl,dhs,dhe,
     262     .                 'box',bilan_dyn_domain_id)
    236263       
    237       call histbeg(trim(infile)//'_'//trim(int2str(mpi_rank)),
     264      call histbeg(trim(infile),
    238265     .             1, rlong(jjb:jje), jjn, rlatg(jjb:jje),
    239266     .             1, 1, 1, jjn,
    240      .             tau0, zjulian, dt_cum, thoriid, fileid)
     267     .             tau0, zjulian, dt_cum, thoriid, fileid,
     268     .             bilan_dyn_domain_id)
    241269
    242270C
     
    352380c
    353381      if(icum.EQ.0) then
    354          ps_cum=0.
    355          masse_cum=0.
    356          flux_u_cum=0.
    357          flux_v_cum=0.
    358          Q_cum=0.
    359          flux_vQ_cum=0.
    360          flux_uQ_cum=0.
     382         jjb=jj_begin
     383         jje=jj_end
     384
     385         ps_cum(:,jjb:jje)=0.
     386         masse_cum(:,jjb:jje,:)=0.
     387         flux_u_cum(:,jjb:jje,:)=0.
     388         Q_cum(:,jjb:jje,:,:)=0.
     389         flux_uQ_cum(:,jjb:jje,:,:)=0.
     390         flux_v_cum(:,jjb:jje,:)=0.
     391         if (pole_sud) jje=jj_end-1
     392         flux_v_cum(:,jjb:jje,:)=0.
     393         flux_vQ_cum(:,jjb:jje,:,:)=0.
    361394      endif
    362395
     
    366399
    367400c   accumulation des flux de masse horizontaux
    368       ps_cum=ps_cum+ps
    369       masse_cum=masse_cum+masse
    370       flux_u_cum=flux_u_cum+flux_u
    371       flux_v_cum=flux_v_cum+flux_v
     401      jjb=jj_begin
     402      jje=jj_end
     403
     404      ps_cum(:,jjb:jje)=ps_cum(:,jjb:jje)+ps(:,jjb:jje)
     405      masse_cum(:,jjb:jje,:)=masse_cum(:,jjb:jje,:)+masse(:,jjb:jje,:)
     406      flux_u_cum(:,jjb:jje,:)=flux_u_cum(:,jjb:jje,:)
     407     .                       +flux_u(:,jjb:jje,:)
     408      if (pole_sud) jje=jj_end-1
     409      flux_v_cum(:,jjb:jje,:)=flux_v_cum(:,jjb:jje,:)
     410     .                         +flux_v(:,jjb:jje,:)
     411
     412      jjb=jj_begin
     413      jje=jj_end
     414
    372415      do iQ=1,nQ
    373416        Q_cum(:,jjb:jje,:,iQ)=Q_cum(:,jjb:jje,:,iQ)
     
    396439c    -------------
    397440      do iQ=1,nQ
     441        call Register_Hallo(Q(1,1,1,iQ),ip1jmp1,llm,0,1,1,0,Req)
     442      enddo
     443      call SendRequest(Req)
     444      call WaitRequest(Req)
     445     
     446      jjb=jj_begin
     447      jje=jj_end
     448      if (pole_sud) jje=jj_end-1
     449     
     450      do iQ=1,nQ
    398451         do l=1,llm
    399452            do j=jjb,jje
     
    426479      call convmas_p(flux_u_cum,flux_v_cum,convm)
    427480      CALL vitvert_p(convm,w)
     481
     482      jjb=jj_begin
     483      jje=jj_end
    428484
    429485      do iQ=1,nQ
     
    455511      enddo
    456512      zz=1./float(ncum)
    457       ps_cum=ps_cum*zz
    458       masse_cum=masse_cum*zz
    459       flux_u_cum=flux_u_cum*zz
    460       flux_v_cum=flux_v_cum*zz
    461       flux_uQ_cum=flux_uQ_cum*zz
    462       flux_vQ_cum=flux_vQ_cum*zz
    463       dQ=dQ*zz
     513
     514      jjb=jj_begin
     515      jje=jj_end
     516
     517      ps_cum(:,jjb:jje)=ps_cum(:,jjb:jje)*zz
     518      masse_cum(:,jjb:jje,:)=masse_cum(:,jjb:jje,:)*zz
     519      flux_u_cum(:,jjb:jje,:)=flux_u_cum(:,jjb:jje,:)*zz
     520      flux_uQ_cum(:,jjb:jje,:,:)=flux_uQ_cum(:,jjb:jje,:,:)*zz
     521      dQ(:,jjb:jje,:,:)=dQ(:,jjb:jje,:,:)*zz
     522     
     523      IF (pole_sud) jje=jj_end-1
     524      flux_v_cum(:,jjb:jje,:)=flux_v_cum(:,jjb:jje,:)*zz
     525      flux_vQ_cum(:,jjb:jje,:,:)=flux_vQ_cum(:,jjb:jje,:,:)*zz
     526
     527      jjb=jj_begin
     528      jje=jj_end
    464529
    465530
     
    476541c   cumul zonal des masses des mailles
    477542c   ----------------------------------
    478       zv=0.
    479       zmasse=0.
    480       call massbar(masse_cum,massebx,masseby)
     543      jjb=jj_begin
     544      jje=jj_end
     545      if (pole_sud) jje=jj_end-1
     546
     547      zv(jjb:jje,:)=0.
     548      zmasse(jjb:jje,:)=0.
     549
     550      call Register_Hallo(masse_cum,ip1jmp1,llm,1,1,1,1,Req)
     551      do iQ=1,nQ
     552        call Register_Hallo(Q_cum(1,1,1,iQ),ip1jmp1,llm,0,1,1,0,Req)
     553      enddo
     554
     555      call SendRequest(Req)
     556      call WaitRequest(Req)
     557
     558      call massbar_p(masse_cum,massebx,masseby)
    481559     
    482560      jjb=jj_begin
     
    524602c   ----------------------------------------
    525603
     604      jjb=jj_begin
     605      jje=jj_end
     606      if (pole_sud) jje=jj_end-1
     607     
    526608      zvQ=0.
    527609      psiQ=0.
     
    560642
    561643c   fonction de courant pour la circulation meridienne moyenne
    562       psi=0.
     644      psi(jjb:jje,:)=0.
    563645      do l=llm,1,-1
    564646         do j=jjb,jje
     
    588670      enddo
    589671
    590       call histwrite(fileid,'masse',itau,zmasse(jjb:jje,1)
     672      call histwrite(fileid,'masse',itau,zmasse(jjb:jje,1:llm)
    591673     s   ,jjn*llm,ndex3d)
    592       call histwrite(fileid,'v',itau,zv(jjb:jje,1)
     674      call histwrite(fileid,'v',itau,zv(jjb:jje,1:llm)
    593675     s   ,jjn*llm,ndex3d)
    594       psi=psi*1.e-9
     676      psi(jjb:jje,:)=psi(jjb:jje,:)*1.e-9
    595677      call histwrite(fileid,'psi',itau,psi(jjb:jje,1:llm),
    596678     s               jjn*llm,ndex3d)
     
    603685c   -----------------
    604686
    605       zamasse=0.
     687      zamasse(jjb:jje)=0.
    606688      do l=1,llm
    607689         zamasse(jjb:jje)=zamasse(jjb:jje)+zmasse(jjb:jje,l)
    608690      enddo
    609       zavQ=0.
     691     
     692      zavQ(jjb:jje,:,:)=0.
    610693      do iQ=1,nQ
    611694         do itr=2,ntr
  • LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/gcm.F

    r853 r916  
    404404      t_ops = iecri * daysec
    405405      t_wrt = iecri * daysec
    406       CALL inithist_p(dynhist_file,day_ref,annee_ref,time_step,
    407      .              t_ops, t_wrt, nqmx, histid, histvid)
     406c      CALL inithist_p(dynhist_file,day_ref,annee_ref,time_step,
     407c     .              t_ops, t_wrt, nqmx, histid, histvid)
    408408
    409409      t_ops = iperiod * time_step
    410410      t_wrt = periodav * daysec
    411       CALL initdynav_p(dynhistave_file,day_ref,annee_ref,time_step,
    412      .              t_ops, t_wrt, nqmx, histaveid)
     411c      CALL initdynav_p(dynhistave_file,day_ref,annee_ref,time_step,
     412c     .              t_ops, t_wrt, nqmx, histaveid)
    413413
    414414      dtav = iperiod*dtvr/daysec
  • LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/gr_u_scal_p.F

    r790 r916  
    5050      ijb=ij_begin
    5151      ije=ij_end
    52       if (pole_nord) ijb=ij_begin+1
     52      if (pole_nord) ijb=ij_begin+iip1
    5353     
    5454      DO l=1,nx
  • LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/gr_v_scal_p.F

    r774 r916  
    4949      ijb=ij_begin
    5050      ije=ij_end
    51       if (pole_nord) ijb=ij_begin+1
    52       if (pole_sud)  ije=ij_end-1
     51      if (pole_nord) ijb=ij_begin+iip1
     52      if (pole_sud)  ije=ij_end-iip1
    5353     
    5454      DO l=1,nx
  • LMDZ4/branches/LMDZ4_V3_patches/libf/dyn3dpar/leapfrog_p.F

    r774 r916  
    66#define IO_DEBUG
    77
    8 !#undef CPP_IOIPSL
     8c#undef CPP_IOIPSL
    99
    1010      SUBROUTINE leapfrog_p(ucov,vcov,teta,ps,masse,phis,nq,q,clesphy0,
     
    481481c$OMP MASTER
    482482       call VTb(VThallo)
     483c$OMP END MASTER
     484
    483485       call Register_Hallo(ucov,ip1jmp1,llm,1,1,1,1,TestRequest)
    484486       call Register_Hallo(vcov,ip1jm,llm,1,1,1,1,TestRequest)
     
    496498
    497499       call SendRequest(TestRequest)
     500c$OMP BARRIER
    498501       call WaitRequest(TestRequest)
     502
     503c$OMP MASTER
    499504       call VTe(VThallo)
    500 
     505c$OMP END MASTER
     506c$OMP BARRIER
    501507     
    502508      if (debug) then   
     
    514520     .                reshape(q(:,:,j),(/iip1,jmp1,llm/)))
    515521        enddo       
     522c$OMP BARRIER
    516523      endif
    517 c$OMP END MASTER
    518 c$OMP BARRIER
     524
    519525     
    520526      True_itau=True_itau+1
     
    531537      call VTb(VTcaldyn)
    532538c$OMP END MASTER
    533 c$OMP BARRIER
    534539      var_time=time+iday-day_ini
    535 cc$OMP PARALLEL DEFAULT(SHARED)
    536 cc$OMP+         PRIVATE(rdaym_ini,rdayvrai,ijb,ije,
    537 cc$OMP+                 tppn,tpn,tpps,tps)
    538 
    539 cc$OMP+         SHARED(itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
    540 cc$OMP+                phi,conser,du,dv,dteta,dp,w, pbaru,pbarv,
    541 cc$OMP+                time, iday,day_ini,forward,leapf, iapptrac,
    542 cc$OMP+                q,dq,p,VTcaldyn,offline,dtvr,itau)     
    543 
     540
     541c$OMP BARRIER
     542!      CALL FTRACE_REGION_BEGIN("caldyn")
    544543      CALL caldyn_p
    545544     $  ( itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
    546545     $    phi,conser,du,dv,dteta,dp,w, pbaru,pbarv, time+iday-day_ini )
    547546
    548 ccc$OMP END PARALLEL     
     547!      CALL FTRACE_REGION_END("caldyn")
    549548c$OMP MASTER
    550549      call VTe(VTcaldyn)
    551550c$OMP END MASTER     
     551
     552cc$OMP BARRIER
     553cc$OMP MASTER
    552554c      call WriteField_p('du',reshape(du,(/iip1,jmp1,llm/)))
    553555c      call WriteField_p('dv',reshape(dv,(/iip1,jjm,llm/)))
     
    557559c      call WriteField_p('pbaru',reshape(pbaru,(/iip1,jmp1,llm/)))
    558560c      call WriteField_p('pbarv',reshape(pbarv,(/iip1,jjm,llm/)))
     561c      call WriteField_p('p',reshape(p,(/iip1,jmp1,llmp1/)))
     562c      call WriteField_p('masse',reshape(masse,(/iip1,jmp1,llm/)))
     563c      call WriteField_p('pk',reshape(pk,(/iip1,jmp1,llm/)))
     564cc$OMP END MASTER
    559565
    560566c-----------------------------------------------------------------------
     
    587593       IF (offline) THEN
    588594Cmaf stokage du flux de masse pour traceurs OFF-LINE
    589 #undef CPP_IOIPSL
     595
    590596#ifdef CPP_IOIPSL
    591            CALL fluxstokenc(pbaru,pbarv,masse,teta,phi,phis,
     597           CALL fluxstokenc_p(pbaru,pbarv,masse,teta,phi,phis,
    592598     .   dtvr, itau)
    593599#endif
     
    662668c$OMP END MASTER
    663669
    664 c$OMP BARRIER
    665670         CALL pression_p (  ip1jmp1, ap, bp, ps,  p      )
    666671
     672c$OMP BARRIER
     673
    667674         CALL exner_hyb_p(  ip1jmp1, ps, p,alpha,beta,pks, pk, pkf )
     675c$OMP BARRIER
    668676           rdaym_ini  = itau * dtvr / daysec
    669677           rdayvrai   = rdaym_ini  + day_ini
     
    690698c$OMP MASTER
    691699        call VTb(VThallo)
     700c$OMP END MASTER
     701
    692702        call SetTag(Request_physic,800)
    693703       
     
    726736     *                               jj_Nb_physic,2,2,Request_physic)
    727737#endif
     738       
     739        call SendRequest(Request_Physic)
     740c$OMP BARRIER
     741        call WaitRequest(Request_Physic)       
     742
     743c$OMP BARRIER
     744c$OMP MASTER
    728745        call SetDistrib(jj_nb_Physic)
    729        
    730         call SendRequest(Request_Physic)
    731         call WaitRequest(Request_Physic)       
    732        
    733746        call VTe(VThallo)
    734        
    735747        call VTb(VTphysiq)
    736748c$OMP END MASTER
     
    773785c$OMP MASTER
    774786        call SetDistrib(jj_nb_Physic_bis)
    775 
    776787        call VTb(VThallo)
     788c$OMP END MASTER
     789c$OMP BARRIER
    777790 
    778791        call Register_Hallo(dufi,ip1jmp1,llm,
     
    794807       
    795808        call SendRequest(Request_Physic)
     809c$OMP BARRIER
    796810        call WaitRequest(Request_Physic)
    797811             
     812c$OMP BARRIER
     813c$OMP MASTER
    798814        call VTe(VThallo)
    799  
    800815        call SetDistrib(jj_nb_Physic)
    801816c$OMP END MASTER
    802817c$OMP BARRIER   
    803                 ijb=ij_begin
     818       
     819        ijb=ij_begin
    804820        if (.not. pole_nord) then
    805821       
     
    842858c$OMP MASTER
    843859        call VTe(VTphysiq)
    844 
    845860        call VTb(VThallo)
     861c$OMP END MASTER
    846862
    847863        call SetTag(Request_physic,800)
     
    878894
    879895        call SendRequest(Request_Physic)
     896c$OMP BARRIER
    880897        call WaitRequest(Request_Physic)     
    881898
    882         call VTe(VThallo)
    883 
    884         call SetDistrib(jj_Nb_caldyn)
     899c$OMP BARRIER
     900c$OMP MASTER
     901       call VTe(VThallo)
     902       call SetDistrib(jj_Nb_caldyn)
    885903c$OMP END MASTER
    886904c$OMP BARRIER
     
    931949       ENDIF
    932950
    933 c$OMP BARRIER
    934951        CALL pression_p ( ip1jmp1, ap, bp, ps, p                  )
     952c$OMP BARRIER
     953
     954
    935955        CALL exner_hyb_p( ip1jmp1, ps, p,alpha,beta, pks, pk, pkf )
    936956c$OMP BARRIER
     
    956976
    957977c$OMP BARRIER
    958 c$OMP MASTER
     978
    959979        call Register_SwapFieldHallo(ucov,ucov,ip1jmp1,llm,
    960980     *                          jj_Nb_dissip,1,1,Request_dissip)
     
    973993
    974994        call SendRequest(Request_dissip)       
     995c$OMP BARRIER
    975996        call WaitRequest(Request_dissip)       
     997
     998c$OMP BARRIER
     999c$OMP MASTER
    9761000        call SetDistrib(jj_Nb_dissip)
    977        
    978         call VTe(VThallo)
    979 
     1001        call VTe(VThallo)
    9801002        call VTb(VTdissipation)
    981        
    9821003        call start_timer(timer_dissip)
    9831004c$OMP END MASTER
     
    10161037            call suspend_timer(timer_dissip)
    10171038            call VTb(VThallo)
    1018 
     1039c$OMP END MASTER
    10191040            call Register_Hallo(ucov,ip1jmp1,llm,1,1,1,1,Request_Dissip)
    10201041            call Register_Hallo(vcov,ip1jm,llm,1,1,1,1,Request_Dissip)
    10211042            call SendRequest(Request_Dissip)
     1043c$OMP BARRIER
    10221044            call WaitRequest(Request_Dissip)
     1045c$OMP MASTER
    10231046            call VTe(VThallo)
    10241047            call resume_timer(timer_dissip)
     
    11151138c$OMP MASTER
    11161139        call VTe(VTdissipation)
    1117 
    11181140        call stop_timer(timer_dissip)
    1119        
    11201141        call VTb(VThallo)
    1121 
     1142c$OMP END MASTER
    11221143        call Register_SwapField(ucov,ucov,ip1jmp1,llm,
    11231144     *                          jj_Nb_caldyn,Request_dissip)
     
    11361157
    11371158        call SendRequest(Request_dissip)       
     1159c$OMP BARRIER
    11381160        call WaitRequest(Request_dissip)       
     1161
     1162c$OMP BARRIER
     1163c$OMP MASTER
    11391164        call SetDistrib(jj_Nb_caldyn)
    11401165        call VTe(VThallo)
     
    11421167        print *,'fin dissipation'
    11431168c$OMP END MASTER
     1169c$OMP BARRIER
    11441170      END IF
    11451171
     
    12111237        call finalize_parallel
    12121238c$OMP END MASTER
     1239c$OMP BARRIER
    12131240        RETURN
    12141241      ENDIF
     
    12441271            IF(MOD(itau,iperiod).EQ.0 .OR. itau.EQ.itaufin) THEN
    12451272c$OMP BARRIER
    1246 c$OMP MASTER
    12471273               IF(itau.EQ.itaufin) THEN
    12481274                  iav=1
     
    12501276                  iav=0
    12511277               ENDIF
    1252 #ifdef CPP_IOIPSL
    1253              call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
    1254               call SendRequest(TestRequest)
    1255               call WaitRequest(TestRequest)
    1256 
    1257               CALL writedynav_p(histaveid, nqmx, itau,vcov ,
    1258      ,                          ucov,teta,pk,phi,q,masse,ps,phis)
     1278c#ifdef CPP_IOIPSL
     1279c             call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
     1280c             call SendRequest(TestRequest)
     1281cc$OMP BARRIER
     1282c              call WaitRequest(TestRequest)
     1283cc$OMP MASTER
     1284c              CALL writedynav_p(histaveid, nqmx, itau,vcov ,
     1285c     ,                          ucov,teta,pk,phi,q,masse,ps,phis)
     1286cc$OMP END MASTER
    12591287c               call bilan_dyn_p(2,dtvr*iperiod,dtvr*day_step*periodav,
    12601288c     ,           ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,q)
    1261 #endif
    1262 c$OMP END MASTER
     1289c#endif
    12631290            ENDIF
    12641291
     
    13021329        enddo
    13031330       
    1304 #ifdef CPP_IOIPSL
     1331c#ifdef CPP_IOIPSL
    13051332 
    1306         CALL writehist_p(histid,histvid, nqmx,itau,vcov,
    1307      s                       ucov,teta,phi,q,masse,ps,phis)
     1333c        CALL writehist_p(histid,histvid, nqmx,itau,vcov,
     1334c     s                       ucov,teta,phi,q,masse,ps,phis)
    13081335c#else
    13091336c       call Gather_Field(unat,ip1jmp1,llm,0)
     
    13181345c#include "write_grads_dyn.h"
    13191346c       endif
    1320 #endif
     1347c#endif
    13211348c$OMP END MASTER
    13221349           ENDIF
     
    14021429                  iav=0
    14031430               ENDIF
    1404 #ifdef CPP_IOIPSL
    1405 c$OMP BARRIER
    1406 c$OMP MASTER
    1407 
    1408               call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
    1409               call SendRequest(TestRequest)
    1410               call WaitRequest(TestRequest)
    1411 
    1412               CALL writedynav_p(histaveid, nqmx, itau,vcov ,
    1413      ,                          ucov,teta,pk,phi,q,masse,ps,phis)
     1431c#ifdef CPP_IOIPSL
     1432cc$OMP BARRIER
     1433
     1434c              call Register_Hallo(vcov,ip1jm,llm,1,0,0,1,TestRequest)
     1435c              call SendRequest(TestRequest)
     1436cc$OMP BARRIER
     1437c              call WaitRequest(TestRequest)
     1438
     1439cc$OMP MASTER
     1440c              CALL writedynav_p(histaveid, nqmx, itau,vcov ,
     1441c     ,                          ucov,teta,pk,phi,q,masse,ps,phis)
    14141442c               call bilan_dyn_p (2,dtvr*iperiod,dtvr*day_step*periodav,
    14151443c     ,           ps,masse,pk,pbaru,pbarv,teta,phi,ucov,vcov,q)
    1416 c$OMP END MASTER
    1417 #endif
     1444cc$OMP END MASTER
     1445c#endif
    14181446            ENDIF
    14191447
     
    14511479        enddo
    14521480
    1453 #ifdef CPP_IOIPSL
    1454 
    1455        CALL writehist_p( histid, histvid, nqmx, itau,vcov ,
    1456      ,                           ucov,teta,phi,q,masse,ps,phis)
     1481c#ifdef CPP_IOIPSL
     1482
     1483c       CALL writehist_p( histid, histvid, nqmx, itau,vcov ,
     1484c     ,                           ucov,teta,phi,q,masse,ps,phis)
    14571485c#else
    14581486c      call Gather_Field(unat,ip1jmp1,llm,0)
     
    14671495c#include "write_grads_dyn.h"
    14681496c      endif
    1469 #endif
     1497c#endif
    14701498
    14711499c$OMP END MASTER
Note: See TracChangeset for help on using the changeset viewer.