Changeset 465 for branches/GRISLIv3/SOURCES/initial-phy-2.f90

Timestamp:

01/26/24 10:24:00 (4 months ago)

Author:

aquiquet

Message:

Cleaning branch: start cleaning module3D

File:

• branches/GRISLIv3/SOURCES/initial-phy-2.f90 (modified) (7 diffs)

branches/GRISLIv3/SOURCES/initial-phy-2.f90

@@ -20 +20 @@
   use module3D_phy, only:icompteur,iout,reprcptr,igrdline,schoof,ibmelt_inv,dtmin,dtmax,&
        dtt,testdiag,num_param,num_rep_42,num_ritz,dt,ntmax,time,tafor,sealevel,&
-       sealevel_2d,secyear,nzz,nz,nzm,cl,da,nice,isynchro
+       sealevel_2d,secyear,isynchro
   use geography, only: geoplace
   use runparam, only :runname,itracebug,num_tracebug,tbegin,tend,dirsource,dirnameout,dttest,&
        nt
-  use param_phy_mod, only:rog,g,ro,row,rowg,rom,romg,rofresh,rofreshg,dice,coef_Sflot,&
-       coef_Bflot
+
 
   implicit none
@@ -40 +39 @@
   !--------------------------------------------------------------
   dirsource='../SOURCES'
-  ! dirsource='../../SOURCES'
-  ! filin='../'//trim(dirsource)//'/Fichiers-parametres/'//trim(geoplace)//'_param_list.dat'
   filin=trim(geoplace)//'_param_list.dat'
 
@@ -85 +82 @@
 
 
-
-  !     switch reprise compteur ou non: 1 => compteur, 0 => pas compteur
-  !     icompteur=2 reprise de tout sauf de la topo
-  !     icompteur=1 reprise de tout
-
-  !      reprcptr=TRIM('../CPTR-Heino/Heino096+k150.CPTR')
-
-
-
-
   write(6,*)'_________________________________________________'
   write(6,*)'runname=',runname
-  !      write(6,*) dirnameout
-  !      write(6,*) TRIM(DIRNAMEOUT)//'short'//runname//'.ritz'
 
 
@@ -104 +89 @@
   !------------------------------
 
-  ! filin='time-series/short'//runname//'.ritz'
+
   filin='short'//runname//'.ritz'
   filin=TRIM(DIRNAMEOUT)//TRIM(filin)
@@ -133 +118 @@
 
 
-  !     grounding line fixee
-  !--------------------------
-  !     Pour faire des expériences avec l'épaisseur des ice-shelves fixée
-  !     igrdline=1 !fixé, sinon (libre) igrdline=0
-  !     igrdline=0 
-
-
-
   ! lecture des parametres du run                          block timesteps
   !-----------------------------------------------------------------------
@@ -146 +123 @@
   read(num_param,timesteps)
 
-! dmr&aurel ##  namelist/timesteps/dtmin,dtmax,dtt,testdiag,tbegin,tend
   write(num_rep_42,'(A)')'!___________________________________________________________'  
   write(num_rep_42,'(A)') '&timesteps                 ! bloc time steps'
@@ -195 +171 @@
   secyear= 31556926 ! s /an
 
-!!$!     *** ABLATION
-!!$!     integrating step for positive degree days (degrees)
-!!$      DTP=2.0
-!!$!     number of months in 1 year, st. dev. for temp *)
-!!$      NYEAR=12
-!!$      SIGMA=5.0
-!!$!     proportion of melted water that can refreeze *)
-!!$      CSI=0.6
-!!$!     melting factors for snow and ice 
-!!$      Csnow=0.003
-!!$      Cice=0.008
-!!$!     ct for PDD calculation
-!!$      S22=0.5/SIGMA/SIGMA
-!!$      PY=2*PI/NYEAR
-!!$      PDDCT=DTP/SIGMA/sqrt(2.*PI)/NYEAR*365.
-
-
-! calcul de certain parametres de param_phy car ro et g ne sont plus parameter
-
-rog      = ro*g                !<  ro*g   (glace)
-rowg     = row*g               !<  row*g  (ocean) 
-romg     = rom*g               !<  rom*g  (asthenosphere)
-rofreshg = rofresh*g           !<  rofresh*g 
-dice     = ro/row              !<  ro/row
-
-coef_Sflot = (Row-Ro)/Row      !<   S = coef_Sflot * H + sealevel pour les shelves
-coef_Bflot = -Ro/Row           !<   B = coef_Bflot * H + sealevel pour les shelves
-
   !--------------------- 
 
-
-  !     *** TEMPERATURE IN ICE AND MANTLE
-  !     Total number of vertical grid points, grid step in mantle
-  NZZ=NZ+NZM
-
-  !     specific latent heat of fusion of ice J/Kg
-  CL=3.35E5
-
-  !     DA mantle diffusion
-  DA=4.E7
   !     switch 
-  NICE=1
   isynchro=0