Changeset 225 for trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

Timestamp:

12/21/18 10:01:52 (6 years ago)

Author:

dumas

Message:

format of the parameter file created by the code, possibility to use it in reading for Ant-40

File:

• trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90 (modified) (5 diffs)

trunk/SOURCES/GrIce2sea_files/climat_GrIce2sea_years_perturb_mod.f90

@@ -107 +107 @@
   namelist/clim_snap/nb_snap,time_depart_snaps,ecart_snap,file_smb_snap,massb_time
 
-428 format(A)
   rewind(num_param)                     ! pour revenir au debut du fichier param_list.dat
   read(num_param,clim_smb_T_gen)
 
-  write(num_rep_42,428)'!________________________________________________________________' 
-  write(num_rep_42,428)'!  module climat_Grice2sea_years_mod lecture climat ref          '
+  write(num_rep_42,'(A)')'!  module climat_Grice2sea_years_perturb_mod                     '
   write(num_rep_42,clim_smb_T_gen)
-  write(num_rep_42,428)'! smb_file          = fichier SMB (kg/m2/an)                     '
-  write(num_rep_42,428)'! coef_smb_unit     = coef passage m glace/an  (1/910 ou 1/918)  '
-  write(num_rep_42,428)'! temp_annual_file  = Temp moy annuelle  (°C)                    '
-  write(num_rep_42,428)'!________________________________________________________________'
+  write(num_rep_42,'(A)')'! smb_file          = fichier SMB (kg/m2/an)                     '
+  write(num_rep_42,'(A)')'! coef_smb_unit     = coef passage m glace/an  (1/910 ou 1/918)  '
+  write(num_rep_42,'(A)')'! temp_annual_file  = Temp moy annuelle  (°C)                    '
+  write(num_rep_42,'(A)')'!________________________________________________________________'
 
   
@@ -164 +162 @@
   Tjuly(:,:) = Tann(:,:)
 
-  rewind(num_param)        ! pour revenir au debut du fichier param_list.dat
-  read(num_param,clim_snap)
-
-! formats pour les ecritures dans 42
   rewind(num_param)                     ! pour revenir au debut du fichier param_list.dat
   read(num_param,clim_snap)
 
-  write(num_rep_42,428)'!_______________________________________________________________________' 
-  write(num_rep_42,428)'!  module climat_Grice2sea_years_mod                                    '
+  write(num_rep_42,'(A)')'!  module climat_Grice2sea_years_perturb_mod                            '
   write(num_rep_42,clim_snap)
-  write(num_rep_42,428)'! nb_snap           = nombre de snapshots                               '
-  write(num_rep_42,428)'! time_depart_snaps = debut du forçage                                  '
-  write(num_rep_42,428)'! ecart_snap        = ecart entre les snapshots                         '
-  write(num_rep_42,428)'! file_smb_snap     = fichier serie temp anomalie SMB de GCM            '
-  write(num_rep_42,428)'! massb_time        = 0:fixe, 1:interp snapshots, 2:snapsh+interp vert  '
-  write(num_rep_42,428)'!_______________________________________________________________________' 
+  write(num_rep_42,'(A)')'! nb_snap           = nombre de snapshots                               '
+  write(num_rep_42,'(A)')'! time_depart_snaps = debut du forçage                                  '
+  write(num_rep_42,'(A)')'! ecart_snap        = ecart entre les snapshots                         '
+  write(num_rep_42,'(A)')'! file_smb_snap     = fichier serie temp anomalie SMB de GCM            '
+  write(num_rep_42,'(A)')'! massb_time        = 0:fixe, 1:interp snapshots, 2:snapsh+interp vert  '
+  write(num_rep_42,'(A)')'!_______________________________________________________________________' 
 
   if (massb_time == 1) then ! lecture des snapshots
@@ -302 +295 @@
   bpert(:)=max(1.+bpert(:)*coefbmb,0.01) !freezing is not allowed
 
-  
-! ecriture de verification
-!file_smb_snap = 'test_sortie_smb.nc'
-!file_smb_snap  = trim(dirnameout)//trim(file_smb_snap)
-
-!ncidloc = 501
-
-!call  lect_netcdf_type
-!write(6,*) 'ncdf_type'
-
-!if (ncdf_type.eq.32) then
-!   status  = nf90_create(TRIM(file_smb_snap),NF90_WRITE,ncidloc)    ! ouverture du fichier
-!else if (ncdf_type.eq.64) then
-!   status  = nf90_create(trim(file_smb_snap),and(nf90_write,nf90_64bit_offset),ncidloc) ! r2d2
-!end if
-
-!   status  = nf90_close(ncidloc)                                 ! fermeture
-
-! call write_ncdf_dim('x',trim(file_smb_snap),nx)               ! dimensions des tableaux
-! call write_ncdf_dim('y',trim(file_smb_snap),ny)
-! call write_ncdf_dim('time',trim(file_smb_snap),nb_snap)
-
-!Tab3D(:,:,:) = smb_snap (:,:,:)
-!dimtestname(1) = 'x'
-!dimtestname(2) = 'y'
-!dimtestname(3) = 'time'
-
-!call write_ncdf_var(trim('smb'),dimtestname,trim(file_smb_snap),tab3D,'double')
+  if (massb_time == 1) then ! lecture gradients smb
+     call init_grad_smb
+  endif
    
 end subroutine input_clim
@@ -346 +314 @@
   read(num_param,lapse_rates)
 
-! formats pour les ecritures dans 42
-428 format(A)
-
-  rewind(num_param)                     ! pour revenir au debut du fichier param_list.dat
-  read(num_param,lapse_rates)
-
-  write(num_rep_42,428)'!________________________________________________________________' 
-  write(num_rep_42,428)'!  module climat_Grice2sea_years_mod                             '
+  write(num_rep_42,'(A)')'!  module climat_Grice2sea_years_perturb_mod                     '
   write(num_rep_42,lapse_rates)
-  write(num_rep_42,428)'!T_lapse_rate       = lapse rate temp annuelle                   '
-  write(num_rep_42,428)'!________________________________________________________________' 
+  write(num_rep_42,'(A)')'!T_lapse_rate       = lapse rate temp annuelle                   '
+  write(num_rep_42,'(A)')'!________________________________________________________________' 
 
   
@@ -362 +323 @@
   read(num_param,clim_pert_massb)
   
-  write(num_rep_42,428)'!___________________________________________________________'
-  write(num_rep_42,428) '&clim_pert                      ! module climat_perturb_mod '
-  write(num_rep_42,*)
-  write(num_rep_42,*) 'coefT        = ', coefT
-  write(num_rep_42,'(A,A)') ' filforc      = ', filforc
-  write(num_rep_42,*) 'pertsmb        = ', pertsmb
-  write(num_rep_42,*) 'rapsmb         = ', rapsmb
-  write(num_rep_42,*) 'pertbmb        = ', pertbmb
-  write(num_rep_42,*) 'coefbmb         = ', coefbmb
-  write(num_rep_42,*)'/'
-  write(num_rep_42,*)
-
-
-! appelle la routine de lecture des smb annuels
-  call input_clim
-  if (massb_time == 1) then ! lecture gradients smb
-     call init_grad_smb
-  endif
+  write(num_rep_42,'(A)') '! module climat_Grice2sea_years_perturb_mod                      '
+  write(num_rep_42,clim_pert_massb)
+  write(num_rep_42,'(A)') '! coefT : coef amplitude perturbation T'
+  write(num_rep_42,'(A)') '! filforc : fichier de forcage temporel'
+  write(num_rep_42,'(A)') '! pertsmb : 0 SMB fixe, 1 SMB exp temperature'
+  write(num_rep_42,'(A)') '! rapsmb : coef variation SMB'
+  write(num_rep_42,'(A)') '! pertbmb : 0 bmelt fixe, 1 bmelt variable'
+  write(num_rep_42,'(A)') '! coefbmb : coef variation bmelt'
+  write(num_rep_42,'(A)') '!________________________________________________________________'
 
   return