Changeset 108 for palm/trunk/SOURCE/init_pegrid.f90

Timestamp:

Aug 24, 2007 3:10:38 PM (17 years ago)

Author:

letzel

Message:

• Improved coupler: evaporation - salinity-flux coupling for humidity = .T.,

avoid MPI hangs when coupled runs terminate, add DOC/app/chapter_3.8;

• Optional calculation of km and kh from initial TKE e_init;
• Default initialization of km,kh = 0.00001 for ocean = .T.;
• Allow data_output_pr= q, wq, w"q", w*q* for humidity = .T.;
• Bugfix: Rayleigh damping for ocean fixed.

File:

• palm/trunk/SOURCE/init_pegrid.f90 (modified) (6 diffs)

palm/trunk/SOURCE/init_pegrid.f90

@@ -4 +4 @@
 ! Actual revisions:
 ! -----------------
+! TEST OUTPUT (TO BE REMOVED) logging mpi2 ierr values
 ! Intercommunicator (comm_inter) and derived data type (type_xy) for
-! coupled model runs created,
+! coupled model runs created, assign coupling_mode_remote
 ! indices nxlu and nysv are calculated (needed for non-cyclic boundary
 ! conditions)
@@ -467 +468 @@
 
 !
-!--    Recieve data from the other PEs
+!--    Receive data from the other PEs
        DO  i = 1,numprocs-1
           CALL MPI_RECV( ibuf, 12, MPI_INTEGER, i, MPI_ANY_TAG, comm2d, status, &
@@ -494 +495 @@
 
           CALL MPI_OPEN_PORT( MPI_INFO_NULL, port_name, ierr )
+!
+!--       TEST OUTPUT (TO BE REMOVED)
+          WRITE(9,*)  TRIM( coupling_mode ),  &
+               ', ierr after MPI_OPEN_PORT: ', ierr
+          CALL LOCAL_FLUSH( 9 )
+
           CALL MPI_PUBLISH_NAME( 'palm_coupler', MPI_INFO_NULL, port_name, &
                                  ierr )
+!
+!--       TEST OUTPUT (TO BE REMOVED)
+          WRITE(9,*)  TRIM( coupling_mode ),  &
+               ', ierr after MPI_PUBLISH_NAME: ', ierr
+          CALL LOCAL_FLUSH( 9 )
+
 !
 !--       Write a flag file for the ocean model and the other atmosphere
@@ -519 +532 @@
 
           CALL MPI_LOOKUP_NAME( 'palm_coupler', MPI_INFO_NULL, port_name, ierr )
+!
+!--       TEST OUTPUT (TO BE REMOVED)
+          WRITE(9,*)  TRIM( coupling_mode ),  &
+               ', ierr after MPI_LOOKUP_NAME: ', ierr
+          CALL LOCAL_FLUSH( 9 )
+
 
        ENDIF
@@ -535 +554 @@
        print*, '--- ierr = ', ierr
        print*, '--- comm_inter atmosphere = ', comm_inter
+
+       coupling_mode_remote = 'ocean_to_atmosphere'
 
     ELSEIF ( coupling_mode == 'ocean_to_atmosphere' )  THEN
@@ -544 +565 @@
        print*, '--- ierr = ', ierr
        print*, '--- comm_inter ocean      = ', comm_inter
+
+       coupling_mode_remote = 'atmosphere_to_ocean'
 
     ENDIF