Changeset 2000 for trunk/ICOSA_LMDZ/src/phyvenus

Timestamp:

Sep 19, 2018, 10:53:15 AM (6 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/phyvenus/interface_icosa_lmdz.f90 (modified) (9 diffs)

trunk/ICOSA_LMDZ/src/phyvenus/interface_icosa_lmdz.f90

@@ -38 +38 @@
   
   CHARACTER(len=30),SAVE,ALLOCATABLE :: tname(:) ! tracer names
+  CHARACTER(len=33),SAVE,ALLOCATABLE :: ttext(:) ! tracer long name for diagnostics
   REAL,SAVE :: pday ! number of ellapsed sols since Ls=0
   REAL,SAVE :: ptime ! "universal time" as fraction of sol (e.g. 0.5 for noon)
@@ -213 +214 @@
   CHARACTER(len=4)              :: type_trac
 !  CHARACTER(len=30),ALLOCATABLE :: tname(:)    ! tracer short name for restart and diagnostics
-  CHARACTER(len=33),ALLOCATABLE :: ttext(:)     ! tracer long name for diagnostics
+!  CHARACTER(len=33),ALLOCATABLE :: ttext(:)     ! tracer long name for diagnostics
   TYPE(t_field),POINTER,SAVE    :: f_ind_cell_glo(:)
   
@@ -259 +260 @@
     ENDDO
 
-!$OMP BARRIER
       
     CALL transfer_icosa_to_lmdz(f_ind_cell_glo,ind_cell_glo)
@@ -284 +284 @@
 
     ! Initialize tracer names, numbers, etc. for physics
+
+! init tracers model for standard lmdz case
 !$OMP MASTER
-
-! init tracers model for standard lmdz case
     ALLOCATE(tname(nqtot))
     ALLOCATE(ttext(nqtot))
+!$OMP END MASTER
+!$OMP BARRIER
 
 ! 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
@@ -306 +310 @@
 !$OMP END MASTER
 !$OMP BARRIER
+    ENDIF ! of (is_mpi_root)
+
+    DO i=1,nqtot
+      CALL bcast(tname(i))
+      CALL bcast(ttext(i))
+    ENDDO
 
    ! Get/set some constants for the physics
 
-!$OMP MASTER
 
     startphy_file=.true.
@@ -326 +335 @@
       ! Read in some information from the startphy.nc file
 
+      IF (is_mpi_root) THEN
+!$OMP MASTER      
       status=nf90_open('startphy.nc',NF90_NOWRITE,ncid)
       if (status.ne.nf90_noerr) then
@@ -347 +358 @@
 
       status=nf90_close(ncid)
+!$OMP END MASTER      
+!$OMP BARRIER
+      ENDIF ! of !IF (is_mpi_root)
+      
+      CALL bcast(day_ini)
+      CALL bcast(annee_ref)
+      CALL bcast(ptime)
 
     ELSE
@@ -365 +383 @@
     CALL getin('iflag_trac',iflag_trac)
 
-!$OMP END MASTER
-!$OMP BARRIER
 
     ! Initialize some physical constants
@@ -597 +613 @@
 
     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