Changeset 86 for trunk/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90

Timestamp:

09/15/16 17:52:05 (8 years ago)

Author:

dumas

Message:

Bug correction in flottab : call to determin_tache suppressed | sealevel added in Netcdf output class 1

File:

• trunk/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90 (modified) (7 diffs)

trunk/SOURCES/Netcdf-routines/sortie_netcdf_GRISLI_mod.0.2-hassine.f90

@@ -111 +111 @@
                                                              !< des variables dans chaque fichier 
   integer,dimension(:),allocatable:: idef_time               !< 1 ou 0 si le temp est defini ou non
+  integer,dimension(:),allocatable:: idef_sealevel           !< 1 ou 0 si sealevel est defini ou non
   integer,dimension(:,:),allocatable:: idef                  !< 1 ou 0 si la varaible est defini ou non
   integer,dimension(:),allocatable :: num_ncdf_file          !< compteur des fichiers netcdf par class  
@@ -504 +505 @@
     if (.not.allocated(nrecs) .and. .not.allocated(nbsnap) .and. .not.allocated(num_ncdf_file) &
          .and. .not. allocated(idef) .and. .not. allocated(idef_time)) then
-       allocate(nrecs(nclassout),nbsnap(nclassout),num_ncdf_file(nclassout),idef(nclassout,ntab),idef_time(nclassout))
+       allocate(nrecs(nclassout),nbsnap(nclassout),num_ncdf_file(nclassout),idef(nclassout,ntab),idef_time(nclassout),idef_sealevel(nclassout))
        nrecs=1
        idef=0
        idef_time=0
+       idef_sealevel=0
        num_ncdf_file=0
     end if
@@ -516 +518 @@
           idef(j,:)=0
           idef_time(j)=0
+          idef_sealevel(j)=0
           nbsnap(j)=0
           ! numerote le fichier sortie
@@ -558 +561 @@
        idef(posis,:)=0
        idef_time(posis)=0
+       idef_sealevel(posis)=0
        nbsnap(posis)=0
        ! numerote le fichier sortie
@@ -654 +658 @@
     character(len=20) :: nametmp                         !< nom intermediaire
     real*8,pointer,dimension(:) :: liste_time => null()  !< liste des snapshot des variables ecrites en netcdf
+    real*8,pointer,dimension(:) :: sealevel_p => null()  !< pointeur vers sealevel pour ecriture netcdf
     real*8,pointer,dimension(:) :: x,y,x1,y1,z,nzzm   
     real*8,pointer,dimension(:,:):: lat,lon => null()
@@ -683 +688 @@
        liste_time(1)=-1
     end if
+    if (.not.associated(sealevel_p)) then
+       allocate(sealevel_p(1))
+       sealevel_p(1)=-1
+    end if
 
 liste_times:  if ((liste_time(1) .ne.timetmp) .or.(liste_time(1) .eq. -1) ) then
        liste_time(1)= timetmp
+       sealevel_p(1)= sealevel
        ! print*,"time outncdf=",liste_time(1)
        write(charint,'(i0)') floor(timetmp)
@@ -720 +730 @@
 ! ecrit le temps
           call write_ncdf_var('time','time',trim(fil_sortie(k)),liste_time,nbsnap(k)+1,idef_time(k),'double')
+          
+          sealevel_p=sealevel
+          call write_ncdf_var('sealevel','time',trim(fil_sortie(k)),sealevel_p,nrecs(k),idef_sealevel(k),'double')
 
           fait = .FALSE.