Changeset 1211 for LMDZ4/branches/LMDZ4-dev/libf/phylmd

Timestamp:

Jul 17, 2009, 5:45:18 PM (15 years ago)

Author:

musat

Message:

Correction indices region;
type_ecri (ave,ave,ave,inst,inst);
freq. ecriture zoutm (secondes) en attendant les corrections de FH/AI
IM

File:

• LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90 (modified) (9 diffs)

LMDZ4/branches/LMDZ4-dev/libf/phylmd/phys_output_mod.F90

@@ -508 +508 @@
    type_ecri_files(1) = 'ave(X)'
    type_ecri_files(2) = 'ave(X)'
-   type_ecri_files(3) = 'inst(X)'
+   type_ecri_files(3) = 'ave(X)'
    type_ecri_files(4) = 'inst(X)'
-   type_ecri_files(5) = 'ave(X)'
+   type_ecri_files(5) = 'inst(X)'
 
    clef_files(1) = ok_mensuel
@@ -546 +546 @@
    type_ecri(:) = type_ecri_files(:)
 
+   print*,'phys_out_lonmin=',lonmin_reg
+   print*,'phys_out_lonmax=',lonmax_reg
+   print*,'phys_out_latmin=',latmin_reg
+   print*,'phys_out_latmax=',latmax_reg
+   print*,'phys_out_filenames=',name_files
+   print*,'phys_out_filetypes=',type_ecri_files
+   print*,'phys_out_filekeys=',clef_files
+   print*,'phys_out_filelevels=',lev_files
+   print*,'phys_out_filetimesteps=',ecrit_files
+
 !!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
 ! Appel de histbeg et histvert pour creer le fichier et les niveaux verticaux !!
@@ -568 +578 @@
     IF (clef_files(iff)) THEN
       
-      zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit
+      zoutm(iff) = ecrit_files(iff)*86400. ! Frequence ou l on ecrit
 
       idayref = day_ref
@@ -582 +592 @@
         jmax_ins=jjmp1
 
-! correction abderr        do i=1,iim-1
+! correction abderr        
         do i=1,iim
            print*,'io_lon(i)=',io_lon(i)
            if (io_lon(i).le.lonmin_reg(iff)) imin_ins=i
-           if (io_lon(i).le.lonmax_reg(iff)) imax_ins=i
+           if (io_lon(i).le.lonmax_reg(iff)) imax_ins=i+1
         enddo
 
-        do j=1,jjmp1-1
+        do j=1,jjmp1
             print*,'io_lat(j)=',io_lat(j)
             if (io_lat(j).ge.latmin_reg(iff)) jmax_ins=j+1
@@ -595 +605 @@
         enddo
 
-         print*,'On stoke le fichier hist sur, ', &
+        print*,'On stoke le fichier histoire numero ',iff,' sur ', &
          imin_ins,imax_ins,jmin_ins,jmax_ins
-         print*,'On stoke le fichier instantanne sur, ', &
+         print*,'longitudes : ', &
          io_lon(imin_ins),io_lon(imax_ins), &
-         io_lat(jmin_ins),io_lat(jmax_ins)
+         'latitudes : ', &
+         io_lat(jmax_ins),io_lat(jmin_ins)
 
  CALL histbeg(name_files(iff),iim,io_lon,jjmp1,io_lat, &
@@ -818 +829 @@
 !      on ecrit u v t q a 850 700 500 200 au niv 3
 
-   type_ecri(1) = 'inst(X)'
-   type_ecri(2) = 'inst(X)'
-   type_ecri(3) = 'inst(X)'
+   type_ecri(1) = 'ave(X)'
+   type_ecri(2) = 'ave(X)'
+   type_ecri(3) = 'ave(X)'
    type_ecri(4) = 'inst(X)'
    type_ecri(5) = 'inst(X)'
@@ -1010 +1021 @@
          ENDIF ! clef_files
 
-         ENDDO ! 
+         ENDDO !  iff
       end subroutine phys_output_open
 
@@ -1034 +1045 @@
        real zstophym
 
-       if (type_ecri(iff)=='inst(X)') then
+       if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
          zstophym=zoutm(iff)
        else
@@ -1074 +1085 @@
        call conf_physoutputs(nomvar,flag_var)
 
-       if (type_ecri(iff)=='inst(X)') then
+       if (type_ecri(iff)=='inst(X)'.OR.type_ecri(iff)=='once') then
          zstophym=zoutm(iff)
        else