Changeset 1318 for LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/phys_cosp.F90

Timestamp:

Feb 23, 2010, 9:12:08 PM (14 years ago)

Author:

yann meurdesoif

Message:

YM : Parallelisation COSP MPI+OpenMP

File:

• LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/phys_cosp.F90 (modified) (10 diffs)

LMDZ4/branches/LMDZ4V5.0-dev/libf/cosp/phys_cosp.F90

@@ -70 +70 @@
   USE MOD_COSP
   USE mod_phys_lmdz_para
+  USE mod_grid_phy_lmdz
   use ioipsl
   use iophy
@@ -76 +77 @@
 
   ! Local variables
-  character(len=64)  :: cosp_input_nl='cosp_input_nl.txt'
-  character(len=64)  :: cosp_output_nl='cosp_output_nl.txt'
+  character(len=64),PARAMETER  :: cosp_input_nl='cosp_input_nl.txt'
+  character(len=64),PARAMETER  :: cosp_output_nl='cosp_output_nl.txt'
   character(len=512), save :: finput ! Input file name
   character(len=512), save :: cmor_nl
@@ -84 +85 @@
   integer,parameter :: Ncollmdz=20
   integer, save :: Npoints      ! Number of gridpoints
+!$OMP THREADPRIVATE(Npoints)
   integer, save :: Nlevels      ! Number of levels
   Integer :: Nptslmdz,Nlevlmdz ! Nb de points issus de physiq.F
@@ -90 +92 @@
   integer :: i
   type(cosp_config),save :: cfg   ! Configuration options
+!$OMP THREADPRIVATE(cfg)
   type(cosp_gridbox) :: gbx ! Gridbox information. Input for COSP
   type(cosp_subgrid) :: sgx     ! Subgrid outputs
@@ -103 +106 @@
   integer :: Nlon,Nlat,geomode
   real,save :: radar_freq,k2,ZenAng,co2,ch4,n2o,co,emsfc_lw
+!$OMP THREADPRIVATE(emsfc_lw)
   integer,dimension(RTTOV_MAX_CHANNELS),save :: Channels
   real,dimension(RTTOV_MAX_CHANNELS),save :: Surfem
@@ -115 +119 @@
   integer :: nhori,nvert,nvertp,nvertisccp,nvertm,nvertcol
   integer, save :: nid_day_cosp,nid_mth_cosp,nid_hf_cosp
+!$OMP THREADPRIVATE(nid_day_cosp,nid_mth_cosp,nid_hf_cosp)
   logical, save :: debut_cosp=.true.
+!$OMP THREADPRIVATE(debut_cosp)
   integer :: itau_wcosp
-  character(len=10),dimension(Ncollmdz) :: chcol=(/'c01','c02','c03','c04','c05','c06','c07','c08','c09','c10', &
+  character(len=10),dimension(Ncollmdz),parameter :: chcol=(/'c01','c02','c03','c04','c05','c06','c07','c08','c09','c10', &
                                                    'c11','c12','c13','c14','c15','c16','c17','c18','c19','c20'/)
   real,dimension(Ncollmdz) :: column_ax
   integer, save :: Nlevout
+!$OMP THREADPRIVATE(Nlevout)
 
   include "dimensions.h"
@@ -139 +146 @@
 !
    namelist/COSP_INPUT/cmor_nl,overlap,isccp_topheight,isccp_topheight_direction, &
-              npoints,npoints_it,ncolumns,nlevels,use_vgrid,nlr,csat_vgrid,finput, &
+              npoints_it,ncolumns,nlevels,use_vgrid,nlr,csat_vgrid,finput, &
               radar_freq,surface_radar,use_mie_tables, &
               use_gas_abs,do_ray,melt_lay,k2,Nprmts_max_hydro,Naero,Nprmts_max_aero, &
@@ -154 +161 @@
 
  if (debut_cosp) then
+  NPoints=Nptslmdz
 ! Lecture du namelist input 
-  open(10,file=cosp_input_nl,status='old')
-  read(10,nml=cosp_input)
-  close(10)
+  CALL read_cosp_input
+
 ! Clefs Outputs 
   call read_cosp_output_nl(cosp_output_nl,cfg)
 
-    if ( (Ncollmdz.ne.Ncolumns).or.(Nptslmdz.ne.Npoints).or.(Nlevlmdz.ne.Nlevels) ) then
+    if ( (Ncollmdz.ne.Ncolumns).or. (Nlevlmdz.ne.Nlevels) ) then
        print*,'Nb points Horiz, Vert, Sub-col passes par physiq.F = ', &
                Nptslmdz, Nlevlmdz, Ncollmdz
@@ -169 +176 @@
        call abort
     endif
-
+    
     if (overlaplmdz.ne.overlap) then
        print*,'Attention overlaplmdz different de overlap lu dans namelist '
@@ -462 +469 @@
 !  call system_clock(t1,count_rate,count_max)
 !  print *,(t1-t0)*1.0/count_rate
+ 
+  CONTAINS 
+  
+  SUBROUTINE read_cosp_input
     
+    IF (is_master) THEN
+      OPEN(10,file=cosp_input_nl,status='old')
+      READ(10,nml=cosp_input)
+      CLOSE(10)
+    ENDIF
+    CALL bcast(cmor_nl)
+    CALL bcast(overlap)
+    CALL bcast(isccp_topheight)
+    CALL bcast(isccp_topheight_direction)
+    CALL bcast(npoints_it)
+    CALL bcast(ncolumns)
+    CALL bcast(nlevels)
+    CALL bcast(use_vgrid)
+    CALL bcast(nlr)
+    CALL bcast(csat_vgrid)
+    CALL bcast(finput)
+    CALL bcast(radar_freq)
+    CALL bcast(surface_radar)
+    CALL bcast(use_mie_tables)
+    CALL bcast(use_gas_abs)
+    CALL bcast(do_ray)
+    CALL bcast(melt_lay)
+    CALL bcast(k2)
+    CALL bcast(Nprmts_max_hydro)
+    CALL bcast(Naero)
+    CALL bcast(Nprmts_max_aero)
+    CALL bcast(lidar_ice_type)
+    CALL bcast(use_precipitation_fluxes)
+    CALL bcast(use_reff)
+    CALL bcast(platform)
+    CALL bcast(satellite)
+    CALL bcast(Instrument)
+    CALL bcast(Nchannels)
+    CALL bcast(Channels)
+    CALL bcast(Surfem)
+    CALL bcast(ZenAng)
+    CALL bcast(co2)
+    CALL bcast(ch4)
+    CALL bcast(n2o)
+    CALL bcast(co)
+!$OMP BARRIER  
+  END SUBROUTINE read_cosp_input 
+
 end subroutine phys_cosp