Changeset 3535

Timestamp:

Jun 13, 2019, 10:33:32 AM (5 years ago)

Author:

yann meurdesoif

Message:

For unstructured grid (dynamico).
A file can be used to impose land fraction without using Relief.nc.
The land franction is available at 0.05°, 0.25° and 0.5° degree resolution.
If land_water_0.05.nc or land_water_0.25.nc or land_water_0.50.nc is present in the runing directory, it will be use to impose ter/oce fraction.

YM

File:

• LMDZ6/trunk/libf/phylmd/create_etat0_unstruct.F90 (modified) (1 diff)

LMDZ6/trunk/libf/phylmd/create_etat0_unstruct.F90

@@ -18 +18 @@
    ! for coupling activate ocean fraction reading from file "ocean_fraction.nc"
     IF (is_omp_master) THEN
-      IF(NF90_OPEN("ocean_fraction.nc", NF90_NOWRITE, file_id)==NF90_NOERR) THEN
+
+      IF (NF90_OPEN("ocean_fraction.nc", NF90_NOWRITE, file_id)==NF90_NOERR) THEN
         CALL xios_set_file_attr("frac_ocean",enabled=.TRUE.)
         CALL xios_set_field_attr("mask",field_ref="frac_ocean_read")
         iret=NF90_CLOSE(file_id)
+      ELSE IF (NF90_OPEN("land_water_0.05.nc", NF90_NOWRITE, file_id)==NF90_NOERR) THEN
+        CALL xios_set_file_attr("land_water",name="land_water_0.05",enabled=.TRUE.)
+        CALL xios_set_field_attr("mask",field_ref="land_water")
+        iret=NF90_CLOSE(file_id)
+      ELSE IF (NF90_OPEN("land_water_0.25.nc", NF90_NOWRITE, file_id)==NF90_NOERR) THEN
+        CALL xios_set_file_attr("land_water",name="land_water_0.25",enabled=.TRUE.)
+        CALL xios_set_field_attr("mask",field_ref="land_water")
+        iret=NF90_CLOSE(file_id)
+      ELSE IF (NF90_OPEN("land_water_0.50.nc", NF90_NOWRITE, file_id)==NF90_NOERR) THEN
+        CALL xios_set_file_attr("land_water",name="land_water_0.50",enabled=.TRUE.)
+        CALL xios_set_field_attr("mask",field_ref="land_water")
+        iret=NF90_CLOSE(file_id)
       ENDIF
+
     ENDIF