Changeset 1575 for LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

Timestamp:

Sep 21, 2011, 3:57:48 PM (13 years ago)

Author:

jghattas

Message:

• Added suffix _mpi_rank to lmdz.out text file. Each processus now write into seperate file but only if lunout/=6 (lunout is set in run.def).

• Change some print* into write(lunout,*)

• phytrac.F90 : always include ini_histrac and write_histrac. The file histrac.nc is written if ecrit_tra> 0 (set in physiq.def). Change default value of ecrit_tra into 0.

File:

• LMDZ5/trunk/libf/phylmd/phys_output_mod.F90 (modified) (12 diffs)

LMDZ5/trunk/libf/phylmd/phys_output_mod.F90

@@ -631 +631 @@
     include "thermcell.h"
     include "comvert.h"
+    include "iniprint.h"
 
     real,dimension(klon),intent(in) :: rlon
@@ -687 +688 @@
     real, dimension(nfiles), save     :: phys_out_latmax        = (/ 90., 90., 90., 90., 90., 90. /)
 
-    print*,'Debut phys_output_mod.F90'
+    write(lunout,*) 'Debut phys_output_mod.F90'
     ! Initialisations (Valeurs par defaut
 
@@ -761 +762 @@
     type_ecri_files(:)=type_ecri(:)
 
-    print*,'phys_out_lonmin=',phys_out_lonmin
-    print*,'phys_out_lonmax=',phys_out_lonmax
-    print*,'phys_out_latmin=',phys_out_latmin
-    print*,'phys_out_latmax=',phys_out_latmax
-    print*,'phys_out_filenames=',phys_out_filenames
-    print*,'phys_out_filetypes=',type_ecri
-    print*,'phys_out_filekeys=',clef_files
-    print*,'phys_out_filestations=',clef_stations
-    print*,'phys_out_filelevels=',lev_files
+    write(lunout,*)'phys_out_lonmin=',phys_out_lonmin
+    write(lunout,*)'phys_out_lonmax=',phys_out_lonmax
+    write(lunout,*)'phys_out_latmin=',phys_out_latmin
+    write(lunout,*)'phys_out_latmax=',phys_out_latmax
+    write(lunout,*)'phys_out_filenames=',phys_out_filenames
+    write(lunout,*)'phys_out_filetypes=',type_ecri
+    write(lunout,*)'phys_out_filekeys=',clef_files
+    write(lunout,*)'phys_out_filestations=',clef_stations
+    write(lunout,*)'phys_out_filelevels=',lev_files
 
 !!!!!!!!!!!!!!!!!!!!!!! Boucle sur les fichiers !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
@@ -785 +786 @@
     enddo
     !          if(prt_level.ge.1) then
-    print*,'Ap Hybrid = ',Ahyb(1:klev)
-    print*,'Bp Hybrid = ',Bhyb(1:klev)
-    print*,'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
+    write(lunout,*)'Ap Hybrid = ',Ahyb(1:klev)
+    write(lunout,*)'Bp Hybrid = ',Bhyb(1:klev)
+    write(lunout,*)'Alt approx des couches pour une haut d echelle de 8km = ',Alt(1:klev)
     !          endif
     DO iff=1,nfiles
@@ -799 +800 @@
              call convers_timesteps(chtimestep(iff),dtime,ecrit_files(iff)) 
           endif
-          print*,'ecrit_files(',iff,')= ',ecrit_files(iff)
+          write(lunout,*)'ecrit_files(',iff,')= ',ecrit_files(iff)
 
           zoutm(iff) = ecrit_files(iff) ! Frequence ou l on ecrit en seconde
@@ -817 +818 @@
              ! correction abderr        
              do i=1,iim
-                print*,'io_lon(i)=',io_lon(i)
+                write(lunout,*)'io_lon(i)=',io_lon(i)
                 if (io_lon(i).le.phys_out_lonmin(iff)) imin_ins=i
                 if (io_lon(i).le.phys_out_lonmax(iff)) imax_ins=i+1
@@ -823 +824 @@
 
              do j=1,jjmp1
-                print*,'io_lat(j)=',io_lat(j)
+                write(lunout,*)'io_lat(j)=',io_lat(j)
                 if (io_lat(j).ge.phys_out_latmin(iff)) jmax_ins=j+1
                 if (io_lat(j).ge.phys_out_latmax(iff)) jmin_ins=j
              enddo
 
-             print*,'On stoke le fichier histoire numero ',iff,' sur ', &
+             write(lunout,*)'On stoke le fichier histoire numero ',iff,' sur ', &
                   imin_ins,imax_ins,jmin_ins,jmax_ins
-             print*,'longitudes : ', &
+             write(lunout,*)'longitudes : ', &
                   io_lon(imin_ins),io_lon(imax_ins), &
                   'latitudes : ', &
@@ -843 +844 @@
           else if (clef_stations(iff)) THEN
 
-             print*,'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)
+             write(lunout,*)'phys_output_mod phys_out_filenames=',phys_out_filenames(iff)
 
              call histbeg_phy_all(rlon,rlat,pim,tabij,ipt,jpt,plon,plat,plon_bounds,plat_bounds, &
@@ -1709 +1710 @@
 
     ENDDO !  iff
-    print*,'Fin phys_output_mod.F90'
+    write(lunout,*)'Fin phys_output_mod.F90'
   end subroutine phys_output_open
 
@@ -1841 +1842 @@
     include "temps.h"
     include "comconst.h"
+    include "iniprint.h"
 
     ipos=scan(str,'0123456789.',.true.)
     !  
     il=len_trim(str)
-    print*,ipos,il
+    write(lunout,*)ipos,il
     read(str(1:ipos),*) ttt
-    print*,ttt
+    write(lunout,*)ttt
     type=str(ipos+1:il)
 
@@ -1857 +1859 @@
     if ( type == 'day'.or.type == 'days'.or.type == 'jours'.or.type == 'jour' ) timestep = ttt * dayseconde
     if ( type == 'mounths'.or.type == 'mth'.or.type == 'mois' ) then
-       print*,'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref)
+       write(lunout,*)'annee_ref,day_ref mon_len',annee_ref,day_ref,ioget_mon_len(annee_ref,day_ref)
        timestep = ttt * dayseconde * mth_len
     endif
@@ -1865 +1867 @@
     if ( type == 'TS' ) timestep = ttt * dtime
 
-    print*,'type =      ',type
-    print*,'nb j/h/m =  ',ttt
-    print*,'timestep(s)=',timestep
+    write(lunout,*)'type =      ',type
+    write(lunout,*)'nb j/h/m =  ',ttt
+    write(lunout,*)'timestep(s)=',timestep
 
   END SUBROUTINE convers_timesteps