Changeset 346

Timestamp:

04/30/21 19:25:36 (3 years ago)

Author:

dumas

Message:

nzoc, number of vertical levels in T/S files is read in netcdf file

File:

• trunk/SOURCES/bmelt-beckmann-gcm_mod.f90 (modified) (4 diffs)

trunk/SOURCES/bmelt-beckmann-gcm_mod.f90

@@ -29 +29 @@
   implicit none
 
-  integer, parameter :: nzoc = 40 !< number of vertical levels in the ocean
-  integer, parameter :: nbassins = 10 !< number of vertical levels in the ocean
+  integer, parameter :: nbassins = 10 !< number of sectors in the ocean
+  integer :: nzoc !< number of vertical levels in the ocean (read in netcdf T/S file)
   real*8 :: coef_OM !combined coefficient by DeConto et Pollard (m/yr/C2)
   real*8 :: K_t
   real*8 :: n_tour
-  real*8, dimension (nx,ny,nzoc) :: temp_ocean , mask_oce2, temp_ocean_2 !< thermal forcing            , input
-  real*8, dimension (nx,ny,nzoc) :: salinity_ocean !< thermal forcing            , input
+  real*8, dimension (:,:,:), allocatable :: temp_ocean , mask_oce2, temp_ocean_2 !< thermal forcing            , input
+  real*8, dimension (:,:,:), allocatable :: salinity_ocean !< thermal forcing            , input
+  
   real :: bmelt_empty ! bmelt value for bassins without any GCM data
       
 !  integer, dimension (nx,ny) :: profondeur
   integer, dimension(nx,ny) :: bassin
-  real*8, dimension (nzoc)       :: zoc             !< depth of oceanic levels    , input
+  real*8, dimension (:), allocatable :: zoc             !< depth of oceanic levels    , input
   real*8, dimension (nx,ny)      :: mesh_area       !< grid cell area
   real*8, dimension (nx,ny)      :: ice_draft       !< ice draft (S-H)
   character(len=200) :: TO_file, SO_file, Bassin_file   ! fichiers de forcage
-  logical, dimension(nzoc) :: mask_bassin_vide  ! mask true if bassin without any GCM data
+  logical, dimension(nbassins) :: mask_bassin_vide  ! mask true if bassin without any GCM data
   
 contains
@@ -67 +68 @@
     integer :: depth_max_grid ! max depth level were T S are defined (all grid)
     integer, dimension(nx,ny) :: depth_max ! max depth were T S are defined for every point
-    integer, dimension(nx,ny,nzoc) :: mask_gcm  ! mask were T S are defined
-    integer, dimension(nx,ny,nzoc) :: mask_gcm_init ! mask were T S are defined in GCM before interpolation
+    integer, dimension(:,:,:), allocatable :: mask_gcm  ! mask were T S are defined
+    integer, dimension(:,:,:), allocatable :: mask_gcm_init ! mask were T S are defined in GCM before interpolation
     integer, dimension(nbassins) :: max_depth_bassin ! max level were T S are defined for every bassins
     integer :: nmax ! to avoid inifite loop 
+    integer :: err
     namelist/bmelt_beckmann_gcm_mod/TO_file,SO_file,Bassin_file,bmelt_empty
     
@@ -82 +84 @@
 
     file_inputs=TRIM(DIRNAMEINP)//'Snapshots-GCM/'//TRIM(TO_file)
+    
+    ! read the number of levels in T/S files :
+    nzoc = read_nzoc(file_inputs)
+    ! print*, 'nzoc = ',nzoc
+    
+    ! allocate tabs :
+    if (.not.allocated(zoc)) then
+      allocate(zoc(nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau zoc",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(temp_ocean)) then
+      allocate(temp_ocean(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau temp_ocean",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(mask_oce2)) then
+      allocate(mask_oce2(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau mask_oce2",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(temp_ocean_2)) then
+      allocate(temp_ocean_2(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau temp_ocean_2",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(salinity_ocean)) then
+      allocate(salinity_ocean(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau salinity_ocean",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(mask_gcm)) then
+      allocate(mask_gcm(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau mask_gcm",err
+          stop 4
+       end if
+    end if
+    if (.not.allocated(mask_gcm_init)) then
+      allocate(mask_gcm_init(nx,ny,nzoc),stat=err)
+      if (err/=0) then
+          print *,"erreur a l'allocation du tableau mask_gcm_init",err
+          stop 4
+       end if
+    end if
+     
     call Read_Ncdf_var('lev',file_inputs,tab1d)
     zoc(:)  = -tab1d(:)
@@ -378 +436 @@
 
   end subroutine bmeltshelf
+
+  function read_nzoc(filename)
+    ! Reads number of vertical levels in ocean forcing NetCDF File.
+    ! Returns the levels number 
+
+    ! Reads
+    use netcdf
+
+    implicit none
+    ! input
+    character (len = *), intent(in) :: filename
+    integer :: read_nzoc
+
+    ! This will be the netCDF ID for the file and data variable.
+    integer :: ncid, varid, status
+    integer, dimension(3) :: count
+    integer, dimension(3) :: start
+
+    ! Open the file. 
+    status = nf90_open(filename, nf90_nowrite, ncid)
+    if (status/=nf90_noerr) then
+       write(*,*)"unable to open netcdf file : ",filename
+       stop
+    endif
+
+    ! Get the varids of the netCDF variable.
+    status = nf90_inq_dimid(ncid, "lev", varid)
+    status = nf90_inquire_dimension(ncid, varid, len = read_nzoc)
+
+    ! Close the file. This frees up any internal netCDF resources
+    ! associated with the file.
+    status = nf90_close(ncid)
+
+  end function read_nzoc
+
 end module bmelt_beckmann_gcm_mod