Changeset 2000 for trunk/ICOSA_LMDZ/src/phystd

Timestamp:

Sep 19, 2018, 10:53:15 AM (7 years ago)

Author:

emillour

Message:

Dynamico-physics interfaces:

• some OpenMP fixes.
• small bug fix: in Dynamico vertical flux is positive when upwards, but in the physics positive is downwards, so change the sign when sending values to the physics.

EM

File:

• trunk/ICOSA_LMDZ/src/phystd/interface_icosa_lmdz.f90 (modified) (10 diffs)

trunk/ICOSA_LMDZ/src/phystd/interface_icosa_lmdz.f90 ¶

@@ -255 +255 @@
     ENDDO
 
-!$OMP BARRIER
       
     CALL transfer_icosa_to_lmdz(f_ind_cell_glo,ind_cell_glo)
@@ -276 +275 @@
 
     ! Initialize tracer names, numbers, etc. for physics
-!$OMP MASTER
     !Config  Key  = type_trac
     !Config  Desc = Choix de couplage avec model de chimie INCA ou REPROBUS
@@ -290 +288 @@
 !    nqo=2
 !    nbtr=2
+!$OMP MASTER
     ALLOCATE(tname(nqtot))
+!$OMP END MASTER
+!$OMP BARRIER
 !    ALLOCATE(ttext(nqtot))
 !    ALLOCATE(niadv(nqtot))
@@ -299 +300 @@
 
 ! read tname() from traceur.def file
+    IF (is_mpi_root) THEN
+!$OMP MASTER
     OPEN(unit=42,file="traceur.def",form="formatted",status="old",iostat=ierr)
     IF (ierr==0) THEN
@@ -313 +316 @@
 !$OMP END MASTER
 !$OMP BARRIER
+    ENDIF ! of (is_mpi_root)
+
+    DO i=1,nqtot
+      CALL bcast(tname(i))
+    ENDDO
 
 !    CALL init_infotrac_phy(nqtot,nqo,nbtr,tname,ttext,type_trac,&
@@ -325 +333 @@
 !    CALL inifis(klon_omp,nlayer,nqtot,pdayref,punjours,nday,ptimestep, &
 !            latitude,longitude,cell_area,prad,pg,pr,pcpp)
-!$OMP MASTER
 
     startphy_file=.true.
@@ -332 +339 @@
     IF (startphy_file) THEN
 
+      IF (is_mpi_root) THEN
+!$OMP MASTER      
       status=nf90_open('startfi.nc',NF90_NOWRITE,ncid)
       if (status.ne.nf90_noerr) then
@@ -348 +357 @@
       status=nf90_get_var(ncid,varid,tab_cntrl)
       start_day=tab_cntrl(3)
-      print*,"start_day YORGL",start_day
+      print*,"initialize_physics_omp: start_day",start_day
       status=nf90_close(ncid)
+!$OMP END MASTER      
+!$OMP BARRIER
+      ENDIF ! of !IF (is_mpi_root)
+      
+      CALL bcast(start_day)
 
     ELSE
@@ -364 +378 @@
     ndays=nint(itaumax*(dt/day_length))! number of days to run
     physics_timestep=dt*itau_physics
-!$OMP END MASTER
-!$OMP BARRIER
+
     CALL inifis(klon_omp,llm,nqtot,start_day,day_length,ndays,physics_timestep, &
             latfi,lonfi,airefi,radius,g,kappa*cpp,cpp)
@@ -596 +609 @@
 
     DO l=1,llm
-      wflux_phy(:,l)=wflux_phy(:,l)*cell_area(:)
+      ! Warning: In the physics, vertical flux convention is positive if downwards!
+      wflux_phy(:,l)= - wflux_phy(:,l)*cell_area(:)
+      ! Compute relative geopotential
       phi_phy(:,l)=phi_phy(:,l)-phis_phy(:)
     ENDDO