Context Navigation

Changeset 927

Timestamp:

Apr 3, 2013, 2:20:46 PM (12 years ago)

Author:

slebonnois

Message:

SL: mise a jour de quelques docs + correction bug dans newstart (Venus)

Location:

trunk

Files:

-                      r124
+                      r927
 \end{itemize}
+The Exner function: definition.
+The Exner function:
+$pk = C_p \times (p/preff)^\kappa$.
 It corresponds to the pressure levels within the layers.
 Used for the computation of the potential temperature.
 …
 of the altitude (in km) of the atmospheric level (one per line of the file).
+For planetary applications, the usual way to go is to use \textsf{disvert\_noterre.F} together with \textsf{z2sig.def}.
+The sigma levels are computed as:
+\begin{itemize}
+\item $H$ is the first value in \textsf{z2sig.def}, $zsig$ are the following values in \textsf{z2sig.def}
+\item $\sigma(1)=1$
+\item l=2 to llm: $\sigma(l)= 0.5 \times (\exp(-zsig(l)/H)+\exp(-zsig(l-1)/H))$
+\item $\sigma(llm+1)=0$
+\end{itemize}
+The $ap$ and $bp$ values, defining the interface pressure levels, are then computed using these sigma values, the reference pressure $preff$, and the transition pressure $pa$. The values of $preff$ and $pa$ are read in the \textsf{start.nc} file, in the control array.
 \item Interface pressures:
 computed in \textsf{caldyn0.F, caldyn.F, integrd.F, leapfrog.F}

r819	r927
63	63	\end{itemize}
64	64
65		Currently, this tool is well adapted for Venus and Titan, so the corresponding routine is called \textsf{newstart-VT.F}, located in the \textsf{dyn3d} directory. As for the \textsf{start2archive} tool, it is compiled the same way as the GCM.
	65	Currently, this tool is well adapted for Venus and Titan, so the corresponding routine is called \textsf{newstart-VT.F}, located in the \textsf{dyn3d} directory. As for the \textsf{start2archive} tool, it is compiled the same way as the GCM (in sequential only).
66	66
67	67	As inputs, it takes the \textsf{start\_archive.nc} file with the initial state to be modified, and the \textsf{traceur.def} and \textsf{z2sig.def} files.

r892	r927
52	52	The goal of this tool is to initialize the model the same way \textsf{gcm.F} does on 3D, but only on a single column, so that the physics may be tested without any dynamics.
53	53
54		It can be compiled with a command like (e.g. for 50 layers):
	54	It can be compiled (sequential only) with a command like (e.g. for 50 layers):
55	55
56	56	\textsf{makelmdz -p venus -d 50 rcm1d}

-                      r849
+                      r927
       ip1jmp1old = (imold+1-1/imold)*(jmold+1-1/jmold)
+c dans run.def
+      CALL getin('planet_type',planet_type)
+      call ini_cpdet
 c=======================================================================
 c   CHANGEMENT DE CONSTANTES CONTENUES DANS tab_cntrl
 …
      &                   rlonuold,rlatvold,rlonu,rlatv)
       write (*,*) 'T ', T(1,jjp1,1)  ! INFO
+! pour info:
+! Si extension verticale, la T est extrapolee constante au-dessus de lmold
 c passage grille physique pour restartphy.nc
 …
           do i=1,iip1
             ucov( i,j,l ) = unat( i,j,l ) * cu(i,j)
+c           ucov( i,j,l ) = 0
+! pour info:
+! Si extension verticale, on impose u=0 au-dessus de lmold
+            if (l.gt.lmold) ucov( i,j,l ) = 0
           end do
         end do
 …
           do i=1,iim
             vcov( i,j,l ) = vnat( i,j,l ) * cv(i,j)
+c           vcov( i,j,l ) = 0
+! pour info:
+! Si extension verticale, on impose v=0 au-dessus de lmold
+            if (l.gt.lmold) vcov( i,j,l ) = 0
           end do
           vcov( iip1,j,l ) = vcov( 1,j,l )

Note: See TracChangeset for help on using the changeset viewer.